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ADVANCED  CONCEPTS  THEORY  ANNUAL  REPORT  1983 


INTRODUCTION 

The  objectives  of  the  Advanced  Concepts  Theory  Program  at  the  Naval 
Research  Laboratory  are  to  develop  a  variety  of  theoretical  models,  tools, 
and  computational  strategies  to  understand,  guide,  and  predict  the  behavior 
of  high  brightness,  laboratory  x-ray  plasma  sources  designed  to  simulate 
nuclear  weapons  effects.  Specifically,  to  develop  the  capability  to 
determine  self-consistently  the  radiative  properties  of  different  materials 
In  a  variety  of  geometrical  configurations  that  occur  In  exploding  wire, 
gas  puff,  and  pellet  Implosion  experiments.  These  theoretical  models  must 
treat  hard  and  soft  x-ray  emission  from  different  electron  configurations 
with  K,  L,  and  M  shells,  and  they  must  include  accurate  treatments  of  the 
different  broadening  mechanisms  that  contribute  to  the  emission  and 
absorption  profiles.  The  self-consistency  of  these  nonlinear  models  is 
critical  if  they  are  to  be  used  to  (1)  accurately  diagnose  both  ongoing 
experiments  and  calculations  and  (2)  properly  understand  the  energy  flow 
and  partitioning  that  occurs  during  energy  absorption,  implosion,  and  x-ray 
emission.  Such  models  can  then  be  employed  to  interpret,  explain,  and 
guide  present  day  experiments,  on  the  one  hand,  and  to  predict  how  these 
experiments  should  be  scaled  with  improved  x-ray  conversion  efficiency  to 
higher  energy  and  power  output  levels  on  the  other 

In  this  report  we  will  focus  primarily  on  those  activities  not  dealing 
with  load  dynamics  which  will  appear  in  a  separate  report,  but  instead  will 
focus  on  ion  beam  target  interaction,  ionization  dynamics  and  radiation 
physics.  Since  a  number  of  papers  are  published  each  year  in  the  open 
literature  emphasing  the  eraission/absorption  physics  and  radiation 
transport  aspects  of  the  program,  it  was  decided  to  omit  lengthy  discussion 
of  these  activities  here  and  refer  the  reader  to  the  open  literature  for 
the  details. 

Finally,  this  report  represents  the  combined  efforts  of  a  number  of 
people,  namely,  J.  Davis,  D.  Duston,  P.  Kepple,  J.  Apruzese,  R.  Clark,  J. 
Rogerson,  V.  Jacobs  and  F.  Cochran  (of  Berkeley  Research  Associates).  We 
would  also  like  to  acknowledge  the  assistance  of  Maxwell  Laboratory  for 
providing  the  data  analyzed  in  this  report;  in  particular  to  M.  Gersten  for 
making  the  reduced  spectral  data  available. 

Manuicript  approved  March  2,  1984. 


.  Detailed  Modeling  and  Approximate  Scaling  Laws  for  K-Shell  Line 
Power  in  Z-Plnch  Plasmas 

I.  INTRODUCTION 

Substantial  progress  was  registered  in  two  important  areas  directly 
relating  to  the  x-ray  yield  obtained  from  high-energy  Z  pinch  plasmas.  For 
the  first  time  the  effects  of  fine  structure  opacity  profiles  on  emission 
in  hydrogen-like  lines  have  been  calculated.  These  effects  were  found  to 
be  of  considerable  importance  for  any  future  machines  in  the  50  TW  or 
higher  power  range,  although  the  Voigt  profile  approximation  we  have 
employed  all  along  was  found  to  be  valid  for  present  and  earlier 
generations  of  pulsed-power  devices. 

Additionally,  the  analytic  scaling  law  for  K-shell  yield  which  was 
developed  previously  for  argon  has  been  successfully  extended  to  all  low- 
to-medium  Z  elements.  This  scaling  law  gives  the  K-shell  yield  in 
kilojoules  as  a  function  of  such  plasma  parameters  as  temperature  and 
density.  The  next  step  is  to  couple  such  analytic  expressions  for  accurate 
radiation  yields  to  simplified  hydrodynamics  models  in  the  effort  to 
achieve  a  precise  overall  yield  scaling  law  in  terms  of  machine  parameters. 

II.  DETAILED  OPACITY  MODELING 

The  effects  of  hydrogen-like  line  fine  structure  profiles  may  at  first 
appear  to  be  a  minor,  irrelevant  physical  detail.  However,  it  turns  out  to 
be  of  substantial  importance  for  contemplated  future  machines  where  the 
plasma  optical  depths  would  be  substantially  higher  than  those  achieved  to 
date.  Each  Lyman  line  has  normally  been  modeled  as  a  single  line  whose 
opacity  is  characterized  by  one  Voigt  profile.  However,  the  3ingle  "line" 
is  in  reality  two  lines  separated  by  as  much  as  several  thermal  Doppler 
widths.  The  more  realistic  calculation  of  photon  absorption  therefore 
involves  less  opacity  spread  out  over  a  greater  frequency  range.  The  true 
absorption  is  the  superposition  of  two  Voigt  profiles  whose  separation 
varies  both  as  a  function  of  atomic  number  and  the  exact  Lyman  line  under 
consideration.  Obviously,  the  higher  the  optical  depth,  the  greater  is  the 


photon  absorption  and  the  more  critical  a  precise  description  of  the 
opacity  becomes. 

Complete  details  of  our  investigation  are  presented  in  the  Appendix  A, 
"Influence  of  Lyman  Series  Fine  Structure  Opacity  on  the  K-Shell  Spectrum 
and  Level  Populations  of  Low-to-Medium  Z  Plasmas.”  The  finding  most 
relevant  to  DNA  applications  is  that  the  effect  of  the  reduced  opacity 
would  produce  power  output  enhancements  of  30-50%  in  the  principal 
hydrogen-like  lines  for  plasmas  of  optical  depths  5-10  times  greater  than 
those  produced  with  present  machines.  Thus  the  construction  of  more 
powerful  machines  capable  of  producing  these  larger  plasma  masses  could 
well  result  in  even  -greater  K-shell  yields  than  indicated  by  extrapolations 
from  opacity  modeling  which  does  not  include  this  effect.  The  Lyman-line 
power  enhancements  are  quite  clearly  seen  in  the  line  ratio  curves 
presented  in  Figs.  3  and  5  of  the  Appendix  a. 

An  additional  result  of  this  investigation  is  the  discovery  of  a 
potentially  valuable,  independent  plasma  density  diagnostic.  In  general 
for  typical  Z-pinch  conditions,  spectroscopic  diagnosis  of  plasma  density 
is  far  more  difficult  than  that  of  temperature.  In  many  cases,  Z  pinch 
plasmas  produced  by  DNA  machines  are  subjected  to  a  density  diagnosis  based 
on  matching  the  measured  yield  and  size  with  the  required  density.  As  is 
illustrated  in  Fig.  6  of  the  Appendix  A,  the  emitted  Lyman-n  line  profile 
varies  quite  sensitively  with  density.  To  exploit  this  potential,  however, 
development  of  spectrometers  with  spectral  resolution  of  1-2  m.1  will  be 
required,  regardless  of  the  plasma  element  excited  in  the  device. 

III.  ANALYTIC  SCALING  OF  K-SHELL  YIELDS 

In  the  FY82  final  report  an  analytic  expression  for  predicting  the  K- 
shell  yield  of  argon  gas  puff  implosions  was  presented.  This  expression 
agrees  within  10%  with  both  detailed  code  predictions  and  experimental 
measurements  for  each  of  the  30  spatial  elements  in  the  most  elaborately 
and  completely  diagnosed  argon  shot.  In  terms  of  the  plasma  radiating  time 
t  (in  ns),  the  axial  length  1  (in  cm),  the  diameter  d  (in  mm),  the  ion 
density  N  (cm“^)  and  temperature  T  (in  KeV) ,  the  argon  K-shell  yield  is 
given  as 


YOU)  -  142  (ij)  <})  (--S-I?)1'9  <TT3J 

3.7x10 


.2  .-2.0.  ,  ..  2.0. 
)  exp  (-^ — )  min  ( 1 , — ^— ) 


We  have  generalized  Eq.  (1)  to  apply  to  arbitrary  atomic  number  Z 
through  use  of  basic  atomic  data  and  its  scalings.  Favorable  comparison  of 
yields  for  neon,  aluminum,  and  titanium  as  well  as  argon  has  been  achieved. 

The  temperature  at  which  K-shell  ions  dominate  the  species  of  a  plasma 
of  atomic  number  Z  scales  as  Z^.  Therefore,  the  temperature  functions  of 
Eq.  (1)  -  developed  only  for  argon  with  Z  ■  18  -  may  be  generalized  to 

exp  (-2.0  (jg)2  /T)  min  (1,  2.0  (|g)2  /T)  (2) 

The  line  photon  energies  themselves  for  analogous  K-shell  lines  also  scale 

as  Z  .  Assume  for  the  moment  that  a  plasma  of  atomic  number  Z  emits  the 

same  number  of  photons  as  an  argon  plasma  (Z  =  18).  In  this  case  the  yield 

Z  2 

would  be  changed  by  a  factor  (yg)  •  Therefore  the  proportionality 
constant  of  Eq.  (1)  would  be  changed  to 


142  (I8)2 


To  complete  the  Z-scaling  of  Eq.  (1)  we  must  consider  the  conditions 

which  would  result  In  the  same  number  of  photons  being  emitted  from  the 

plasma.  Collisional  excitation  rates  at  the  same  relative  temperature 
2  -3 

(T  —  Z  )  scale  as  Z  .  An  effectively  thin  plasma  -  one  where  the  photons 
mostly  escape  even  after  many  scatterings  -  emits  photons  proportional  to 


the  following  factors:  first  the  product  of  electron  and  ion  density 
2 

N  N  -  N  Z;  second,  the  volume,  V;  and,  third,  the  collisional  excitation 

6  -3 

rates  W  ~  C^Z 

N  .  ~  V  N  2  Z~2  C.  (4) 

(photons)  II 

In  Eq.  (4)  is  a  constant  -  independent  of  Z.  Eq.  (4)  states  that 

scaling  the  ion  density  as  1  would  produce  the  same  number  of  photons 


from  different  elements  in  the  same  plasma  volume, 
terms  of  Eq.  (1)  should  therefore  be  generalized  as 


The  density  function 


4 


The  power  law  of  1.9  In  Eqs.  (1)  and  (5)  accounts  approximately  for 
opacity  effects  on  line  yield  reduction  -  a  perfectly  thin  plasma  would  be 

9 

characterized  by  an  N"  yield  dependence.  Combining  Eqs.  (2),  (3),  and  (5) 
results  in  the  following  analytic  expression  for  yield,  valid  for  any  low 
to  medium  element  of  atomic  number  Z 

Y(kJ)  -  142  (f/  (£5)  (*)  (t^)2( - 5—  —  )X’9  exp  (-2.0  (f/  /T) 

3.7  xl019(^) 

Z  2 

.  min  (1,  2.0  (~)  /T)  (6) 

We  have  applied  Eq.  (6)  to  the  well-diagnosed  shots  which  have  been 
carried  out  with  neon,  aluminum,  and  titanium  [as  well  as  argon  which  was 
done  last  year].  The  table  summarizes  the  results. 


Table  1 


Comparison 

of  Analytic 

Prediction  (Eq.  6)  and  Measured 

Yields 

Neon 

SHOT 

diam. (mm) 

N7flN. 

T(keV)(1) 

Y(kJ, expression) 

Y(actual) 

1472 

1.0 

5.0xl019 

0.35 

46 

45 

Aluminum 

SHOT 

diam. (mm) 

^-rnu 

T(keV)(2) 3 

(Y(kJ, expression) 

Y(measured) 

174 

1.5 

2.2xl019 

0.69 

17.6 

14.6 

177 

1.2 

4.0xl019 

0.53 

22.2 

29.2 

182 

1.2 

3.4xl019 

0.51 

15.1 

20.4 

183 

3.3 

0. 7xl019 

0.79 

11.7 

6.7 

184 

5.5 

0. 13xl019 

0.76 

1.3 

0.6 

Titanium 

SHOT 

diam. (mm) 

Ni.oa 

T(keV)<3) 

Y(kJ, expression) 

Y(actual) 

360 

1.5 

lxio19 

1.33 

1.08 

0.96 

420 

1.5 

lxlO19 

1.48 

1.35 

1.48 

422 

1.5 

lxlO19 

1.43 

1.26 

1.22 

424 

1.5 

lxlO19 

1.53 

1.45 

1.50 

426 

1.5 

lxlO19 

1.78 

1.91 

1.76 

409 

1.5 

lxlO19 

1.65 

1.68 

2.02 

423 

1.5 

lxlO19 

2.00 

2.31 

2.02 

421 

1.5 

lxlO19 

2.13 

2.52 

2.27 

(1) 

Line  ratio  and  continuum 

temperatures 

agree. 

(2)  For  shots  174,  177,  and  182,  T  is  the  average  of  line  ratio  and 
continuum  temperatures,  for  183  and  184  it  is  the  line  ratio  temperature 
(continuum  not  available). 

(3)  Average  of  continuum  and  line  ratio  temperatures. 


The  sequence  of  titanium  shots  presented  involved  changing  the  peak 
diode  power  from  5.2  TVJ  for  shot  360  to  7.5  TW  for  shot  421.  As  evident, 
essentially  no  change  was  observed  in  either  the  size  or  density  of  the 
resultant  pinched  plasmas.  The  increase  in  yield  apparently  reflects  an 
increase  in  the  temperature  with  peak  diode  power,  and  is  tracked  quite 
well  by  the  temperature  function  of  the  scaling  law,  Eq .  (6). 

The  only  serious  discrepancies  present  in  the  comparison  of  the  table 
are  for  the  large  diameter,  diffuse  A l  shots  183  and  184.  The  presumption 
of  plasma  homogeneity  is  doubtless  responsible  for  most  of  this 
difference.  The  difficulty  and  inherent  inaccuracy  in  assigning  a 
"radiating  diameter"  to  the  entire  3  cm  length  of  plasma  could  specifically 
account  for  most  of  the  differences.  The  diagnosed  line  and  continuum 
temperatures  also  diverged  by  20-40%  in  these  two  cases.  Overall,  the 
predictions  for  the  well-characterized  high  yield  shots  agree  with  the  A”, 
measurements  within  an  average  of  23%  and  this  points  to  the  overall 
validity  of  the  general  Z-scaled  yield  expression,  Eq.  (6). 
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Influence  of  Lyman-series  fine-structure  opacity  on  the  K-shell  spectrum 
and  level  populations  of  low-to-medium-Z  plasmas 

J.  P.  Apruzese,  J  Davis,  D.  Duston,  and  R.  W.  Clark 
Plasma  Radiation  Branch,  Plasma  Physics  Division,  Sami  Research  Laboratory, 

Code  47 20,  Washington,  D  C.  20375 
Received  13  July  1983) 

As  laboratory  plasmas  of  increasing  atomic  number,  temperature,  size,  and/or  density  are  pro¬ 
duced,  it  becomes  likely  that  the  details  of  the  doublet  opacity  profiles  of  the  Lyman  senes  will  in¬ 
fluence  the  Jf-shell  level  populations  and  spectrum.  Accordingly,  we  have  analyzed  these  effects  for 
a  range  of  plasma  parameters,  confined  to  densities  low  enough  for  Stark  broadening  to  be  unimpor¬ 
tant.  An  analytic  model  is  developed  which  predicts  line  power  enhancements  and  level-population 
changes  for  A'-shell  plasmas.  This  model  is  based  upon  photon  escape  probability  and  collisional 
quenching  concepts  and  is  valid  for  plasmas  of  atomic  number  -  13- 2b  Additionally,  an  extensive 
set  of  numerical  calculations  of  line  ratios,  line  profiles,  and  level  populations  has  been  earned  out 
for  A'-sheil  argon  plasmas.  Each  computation  was  performed  both  with  detailed  fine-structure  opa¬ 
city  profiles  and  with  a  single- Voigt-profile  approximation.  The  results  of  these  calculations  may  be 
scaled  for  plasmas  of  atomic  number  other  than  18  using  a  simple  set  of  rules  discussed  in  the  text 


I.  INTRODUCTION 

Recent  years  have  witnessed  dramatic  technological  ad¬ 
vances  in  laboratory  production  of  optically  thick  /f-shell 
plasmas  of  low-to-medium-Z  elements.  Consistent  with 
these  developments,  considerable  theoretical  effort  has 
been  devoted  to  elucidating  the  effects  of  opacity  on  the 
spectra  and  intrinsic  properties  of  such  plasmas. 
Typically — although  not  exclusively — the  hydrogenlike 
and  heliumlike  resonance  lines  of  the  plasma  are  quite  op¬ 
tically  thick  but  the  continua  are  optically  thin  in  the  x- 
ray  region.  Thus  it  is  the  line  opacity  which  has  properly- 
attracted  the  most  concern.  As  is  well  known  from  astro- 
physical  literature,  the  precise  form  of  the  line-opacity 
profile — which  determines  the  line-photon  escape 
probability — plays  a  critical  role  in  the  formation  of  the 
spectrum  and  determination  of  the  level  populations.  It  is 
also  true  that — for  quite  specific  physical  reasons — 
regimes  exit  where  simplified  profiles  such  as  pure 
Doppler  may  be  employed  with  acceptable  accuracy.  For 
instance.  Fig.  5(a)  of  Avrett  and  Hummer1  reveals  little 
difference  in  the  steady-state  source  functions  obtained  for 
Voigt  and  Doppler  profiles  up  to  line-center  optical  depths 
of  —50  with  a  Voigt  broadening  parameter  of  10_:.  One 
laboratory  plasma  counterpart  of  this  situation  would  be  a 
1.5-mm  aluminum  plasma  of  ion  density  —  1019  cm-3  and 
600-eV  temperature.  Detailed  calculations2  have  shown 
that  here,  too,  a  Doppler  profile  is  adequate  as  a  represen¬ 
tation  for  line  opacity.  We  will  consider  Z  scaling  of  all 
specifically  quoted  plasma  conditions  below.  At  higher 
densities  and/or  optical  depths  a  Voigt  profile  is  required 
to  adequately  describe  photon  transport  in  the  line  wings. 
At  still  higher  densities,  Stark  effects  must  be  con¬ 
sidered.3-  '  In  this  paper  we  will  delineate  the  plasma 
conditions  where  the  doublet-opacity  profile  due  to  fine 
structure  in  the  Lyman  lines  must  similarly  be  taken  into 
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account.  We  analyze  the  underlying  physics  with  both  an¬ 
alytic  and  numerical  models,  and  obtain  Z-scaled  expres¬ 
sions  governing  the  safe  use  of  simpler  opacity  approxi¬ 
mations.  In  Sec.  II  an  analytic  model  is  developed  which 
serves  two  basic  functions.  First,  it  provides  a  physical 
basis  for  comprehending  the  trends  which  appear  in  the 
numerical  results.  Second,  approximate  but  useful  Z- 
scaled  formulas  expressing  some  of  the  results  are  ob¬ 
tained.  In  Sec.  Ill  specific  numerical  results  for  an  argon 
plasma  (Z  =  18)  are  presented  along  with  a  prescription 
for  scaling  them  to  plasmas  of  different  atomic  numbers. 

II  ANALYTIC  MODEL 
A.  Assumptions  and  restrictions 

Both  the  analytic  and  numerical  results  contained  in 
this  paper  are  restricted  to  plasma  regions  where  Stark 
broadening  has  only  a  small  effect  on  photon  transport 
and  can  be  neglected.  The  calculations  encompass  at  most 
a  few  hundred  optical  depths  in  Ly  a;  hence,  for  only  the 
first  three  lines  of  the  Lyman  senes  will  photon  transport 
in  the  far  line  wings  be  significant.  Stark  profiles  follow 
approximately  an  inverse  2.0— 2.5  power  law3  as  a  func¬ 
tion  of  frequency  separation  from  line  center;  for  Voigt 
profiles,  the  asymptotic  limit  is  an  inverse  power  of  2.0. 
Therefore,  if  the  Stark  width  is  smaller  than  the  Doppler 
width,  the  Stark  wings  will  generally  be  considerably 
smaller  than  the  Lorentz  wings.  This  critenon  is  present¬ 
ly  adopted  in  determining  the  maximum  density  to  which 
our  calculations  can  be  reasonably  extended.  We  obtain  a 
specific  density  for  argon  from  the  tabulated  Stark  pro¬ 
files  of  Kepple  and  Whitney,'  and  scale  the  criterion 
analytically  for  different  atomic  numbers. 

In  examining  the  Stark  profiles  for  Ly  a,  Ly  0,  and  Ly 
y ,  it  is  clear  that  the  upper  density  limit  is  set  by  the 
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broad,  double-peaked  profile  of  Lv  (3.  For  argon,  an  elec¬ 
tron  density  of  -2a  10::  cm-'  is  required  for  the  Stark 
width  to  equal  the  Doppler  width  at  temperatures  of  max¬ 
imum  Ar  Win  abundance.  Accordingly,  the  analytic 
theory  presented  below  is  valid  only  up  to  this  density  for 
argon  and  the  numerical  calculations  described  in  the  next 
section  have  not  been  earned  to  higher  densities. 

This  critenon  for  the  negligibility  of  Stark  effects  is 
readily  Z-scalable.  Gnem,  Blaha,  and  Kepple'  found  that 
the  Stark  widths  expressed  in  A  for  analogous  lines  from 
different  elements  scale  approximately  as  Z"\V,:/3  near 
maximum  ionic  abundance  temperature.  The  Doppler 
width,  also  expressed  in  A,  scales  as  Z"1  !  at  maximum 
abundance  temperature.  Therefore,  the  Stark  width  for 
the  first  three  Lyman  lines  will  be  less  than  the  Doppler 
width  if 


18 

5 

-v,  1:/’ 

18 

Z 

2  x  10::  )  " 

Z 

or 

.V,  <5.1XlOl5Z,:!  .  (2) 

Consequently,  the  present  results  are  not  valid  at  electron 
densities  higher  than  those  indicated  by  Eq.  (2). 

The  relative  populations  of  the  2 p,/;  and  2 p(/;  fine- 
structure  levels,  which  are  likely  to  prevail  under  different 
plasma  conditions,  have  been  a  subject  of  investigation 
and  debate  for  quite  a  few  years.4"14  Irons13  has  summa¬ 
rized  the  status  of  this  work.  Briefly,  at  the  high-  and 
low-density  limits,  e,  the  ratio  of  the  2pwi  to  2 p}/2  sub- 
level  populations,  will  equal  0.5,  the  statistical  weight  ra¬ 
tio  of  the  states.  Rapid  collisional  mixing  ensures  statisti¬ 
cal  equilibrium  at  high  densities,  and  in  the  low-density 
coronal  limit  there  is  little  pumping  of  the  2 p  states  from 
2s;  therefore,  statistical  populations  prevail  because  only 
this  process  can  cause  departures  from  e=0. 5.  In  the  in¬ 
termediate  density  regime,  e  may  reach  values  as  high  as 
0.7-0. 8.  although  the  relative  effects  of  electron-ion  and 
ion-ion  collisions  are  not  yet  precisely  determined10  and 
render  this  conclusion  still  somewhat  uncertain.  As  point¬ 
ed  out  by  Sampson13  and  Irons,1314  substantial  optical 
depth  will  reduce  any  tendency  for  a  departure  from  sta¬ 
tistical  equilibrium.  In  view  of  the  remaining  uncertain¬ 
ties  as  well  as  the  fact  that  the  plasmas  we  analyze  here 
are  generally  optically  thick,  we  have  assumed  statistical 
population  of  the  j  ~  j  and  j—  y  sublevels  <e  =  0.5)  for 
all  the  plasmas  considered  here,  for  each  hydrogenlike  lev- 
el.  We  also  assume  that  the  ion  thermal  velocities  are  suf¬ 
ficiently  greater  than  the  plasma  differential  motion  to  en¬ 
able  the  plasma  to  be  treated  as  stationary.  Irons13,14  has 
explored  some  of  the  consequences  of  substantial  plasma 
differential  motion  for  Ly-a  opacity. 

B.  Analysis  of  opacity  effects 

The  doublet  splitting  of  the  Lyman  series  is  important 
as  an  opacity  determinant  when  the  separation  of  the 
components  is  comparable  to  the  Doppler  half-width 
(Avd)  of  the  line.  For  specificity,  the  Av0  referred  to 
here  is  the  half-width  of  the  Doppler  line  profile  at  e  " 1  of 
the  central  maximum.  As  pointed  out  by  Irons,13  the  ratio 


of  the  fine-structure  splitting  to  the  Doppler  half-width 
varies  as  Z 1  since  the  characteristic  temperature  for 
hydrogenlike-ion  predominance  vanes  as  Z‘.  Specifically, 
for  Ly  a. 


i  Avl7  /( Av)0  =  5.3  (  — 


O' 


For  Ly  /?,  the  splitting  relative  to  the  Doppler  width  is 
about  one-fourth  the  above  value;  for  Ly  •/,  about  one- 
tenth.  Therefore,  fine-structure  opacity  effects  will  have 
by  far  their  greatest  impact  on  Ly  a,  both  because  of  its 
greater  optical  depth  and  because  the  splitting  is  greater 
relative  to  the  Doppler  line  width.  Our  semiquantitative 
analytic  model  consequently  focuses  on  Ly  a. 

Consider  the  n  =2  hydrogenlike  level,  with  the  sublevels 
assumed  statistically  populated.  For  the  moderate  density 
plasmas  considered  here,  this  level  is  populated  primarily 
by  collisional  excitation  from  the  hydrogenlike  ground 
state  and  by  radiative  excitation  due  to  scattered  Ly  a 
photons.  Depopulation  occurs  by  a  variety  of  processes, 
including  spontaneous  decay,  collisional  excitation  and 
de-excitation,  and  (primarily  collisional)  ionization.  We 
will  follow  the  standard  approximate  procedure  in  ac¬ 
counting  for  radiative  excitation  by  diluting  the  spontane¬ 
ous  decay  rate  A  2|  by  the  spatially  averaged  photon-escape 
probability  P,.  This  counts  only  those  decays  not  result¬ 
ing  in  re-excitation  of  the  n—  2  level.  Letting  C;:  stand 
for  the  collisional  excitation  rate  from  n  =  1  in  cm1 
sec-1,  D;  represents  the  sum  of  all  collisional  processes 
depopulating  the  n  =  2  level;  also  in  cm3  sec"1,  a  qualita¬ 
tively  correct  equation  for  the  steady-state  population  of 
the  n  =  2  level  is 


-'  |Ci2-Ve  =  A ;t--I i\P,  +£>;.%* )  • 

In  Eq.  (4),  A",  and  A\  are  the  n  =  I  and  2  level  population 
densities  (cm"3),  St  is  the  electron  density  (cm"3),  and 
the  other  symbols  are  defined  above.  It  is  useful  to  define 
a  quenching  probability  Pq — the  probability  that  the  n  =  2 
level  is  collisionally  depopulated  during  a  Ly-a  scattering 


P  Q  — D-> /l  A  21  -t-Z?2Af  )  . 


1 5 1 


Combining  Eqs.  (4)  and  (5)  leads  to  the  following  expres¬ 
sion  for  the  n  =2  population  density: 


S, 


Pq 


(61 


The  significance  of  Eq.  (6)  becomes  apparent  when  con¬ 
sidering  P,.  The  escape  probability  obtained  will  depend 
on  whether  the  single-Voigt-profile  approximation  or  the 
true  fine-structure  profile  is  used  in  the  calculation.  Let¬ 
ting  P,  now  stand  for  the  escape  probability,  which  would 
be  obtained  from  the  spurious  single-Voigt-profile  approx¬ 
imation.  and  Pff  refer  to  the  true  fine-structure-opacity 
escape  probability,  the  ratio  of  population  densities  ob¬ 
tained  with  the  correct  opacity  to  that  obtained  with  the 
single-profile  opacity  is 

-Vj(  with  fine  structure)  P,{  1  —  Pq  )  +  Pq 

iVjIsingie  Voigt)  ~  P,f{  1  ~pQ)  +  pQ 
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Equation  (7)  is  obtainable  from  Eq,  (6)  with  the  (excellent) 
assumption  that  ,V,  is  not  affected  by  fine-structure  opaci¬ 
ty.  Note  the  limits  implicit  in  Eq.  (7).  At  very  high  den¬ 
sities,  when  Pq— .1,  there  is  no  difference  in  populations. 
Collisional  processes  control  the  level  population  Also, 
when  P'H—Pq  i  and  Ptfi\—  Pq)  are  both  much  smaller 
than  Pq,  fine  structure  again  has  no  effect  on  .V,.  This 
corresponds  to  the  effectively  thick  case,  where  most  pho¬ 
tons  do  not  escape,  even  after  many  scatterings.  The 
number  of  radiative  excitations,  while  possibly  large,  is 
determined  by  collisional  quenching,  not  by  the  detailed 
behavior  of  P,  or  Ptf.  However,  when  Pq  — 0  the  com¬ 
puted  populations  of  n=2  are  inversely  proportional  to 
the  ratio  of  escape  probabilities.  Radiative  excitation 
dominates  the  level  if  P,  (or  P,/ 1  is  substantially  smaller 
than  unity.  For  fine-structure  opacity  to  be  negligible  (de¬ 
fined  as  having  less  than  a  10%  effect  on  population  of 
n  =2)  we  require 
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range 

(0. 18— 0.571Z2  Ry  covers  a  large  range  of  plasmas  where 
hydrogenlike  species  are  significantly  present.  Combining 
Eqs.  ill!  and  (12)  yields  our  final,  Z-dependent  critenon 
for  less  than  10%  level  population  change  due  to  fine- 
structure  opacity 

9. 3  X  10“  16Z  ~  7,V,  >  Pc  .  (13) 

If  the  electron  density  is  sufficiently  high  to  render  the 
inequality  of  Eq.  U3i  valid,  fine-structure  opacity — as 
compared  to  the  single-Voigt-profile  approximation — will 
affect  the  level  populations  by  less  than  10%.  We  again 
stress  that  Eq.  03)  is  to  be  applied  to  optically  thick  plas¬ 
mas  only — when  the  plasma  is  optically  thin  the  line- 
structure  profiles  cannot  affect  populations.  The  Voigt- 


Pt{  1  —  Egl-i-Fg >0.9F,/1 1  —Pq)  +  0.9Pq  .  (8) 


In  general,  Ptf  will  exceed  Pr.  For  Ly  a,  a  single  line  of 
optical  depth  r0  is  in  effect  replaced  by  two  separate  lines 
of  optical  depths  j r0  and  \t0  for  2p\p_  and  2pt/2,  respec¬ 
tively.  The  escape  probability  at  large  optical  depth  for  a 
Voigt  profile  is  proportional  to  r_l/‘;  hence,  if  there  is  no 
overlap  of  the  doublet  wings 
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profile  escape  probability  Pc  has  been  given  in  approxi¬ 
mate  analytic  form  by  Athay.15 

The  effect  on  line  emission  of  the  doublet  profile  is 
perhaps  of  even  more  interest  and  importance  than  the 
change  in  population  densities.  Even  though  the  popula¬ 
tions  may  not  be  affected,  the  different  effective  escape 
probability  can  surely  influence  the  Ly-a  integrated  inten¬ 
sity.  The  ratio  of  computed  emission  with  fine-structure 
opacity  to  that  obtained  from  a  single  Voigt  profile  is 

If  _  AMfine) A2lPef  (J4) 

f0  (V2(singie  Voigt  lAltPt 


a  1.4.  -  (9) 

As  discussed  by  Irons,14  the  assumption  of  no  wing  over¬ 
lap  fails  for  r0  >  102;  this  failure  will  arise  at  lower  r0  for 
small  Z  and  higher  r0  for  large  Z,  according  to  Eq.  (3)  for 
the  doublet  spacing.  However,  since  we  seek  a  criterion 
for  the  onset  of  fine-structure-opacity  effects,  the  ratio  of 
Eq.  (9)  will  be  presently  assumed.  More  exact  results  are 
given  below  in  the  section  on  numerical  modeling.  Com¬ 
bining  Eqs.  (8)  and  (9)  yields  the  following  criterion  for 
the  negligibility  of  fine-structure-opacity  effects  on  popu¬ 
lation  densities: 

0.38  > P,  ,  (10) 

I-Pq 

or  employing  Eq.  (5) 

0.38  >P,  .  (H) 

^21 

To  cast  Eq.  (II)  into  a  more  useful  form,  note  that 
A21  »5x  10U(Z/18)4  sec-1.  Also,  depopulation  of  the 
n  =2  level  is  temperature,  density,  and  Z  dependent.  Fol¬ 
lowing  Refs.  9—11,  we  will  refer  to  the  temperature  in 
units  of  Z1  Ry=13.6Z:  eV.  For  temperatures  of 
(0. 18— 0.57SZ2  Ry  (0.8  to  2.5  keV  for  argon)  the  depopu¬ 
lation  of  it  =2,  summed  over  the  levels,  has  been  calculat¬ 
ed  from  the  Coulomb-Bom  approximation16  and  is  found 
to  vary  less  than  20%  from  its  valu  ■  at  0.3Z2  Ry,  where 


or,  using  Eq.  (7), 

Z/=£fL+_£<L_  j  P  +  .  ,15) 

h  Pe  [  l~PQ  /  1  ~PQ 

The  quantity  Pq  /( 1  —  Pq)  is  obtainable  as  above  for  the 
temperature  range  (0. 18— 0.571Z2  Ry,  from  Eqs.  (5)  and 
(12).  and  the  expression  for.421 


»2.45x  10-15.V,Z-7  .  (16) 

1  ~pQ 
Therefore, 

If  _p,t  P'+  2.45xl0-|XZ-7  (]7) 

I0  P,  Ptf  +  2A5X  10-i5A',Z-7 

Equation  (17)  is  usable  if  Pe  and  P,f  can  be  calculated  or 
estimated.  If,  as  before,  we  assume  P,f  =  1  4F,  an  approx¬ 
imate  criterion  for  10%  enhancement  of  emission  due  to 
fine  structure  is  obtainable. 


5.2x  10~I?.V,Z“'>  P,  . 


The  physical  content  of  Eq.  (18)  is  that,  in  an  optically 
thick  plasma  with  an  electron  density  sufficiently  high  to 
validate  the  inequality,  fine-structure  opacity  will  result  in 
a  Ly-a  intensity  at  least  10%  higher  than  would  be  ob¬ 
tained  from  the  single  Voigt  approximation. 

We  have  demonstrated  that  the  interplay  of  optical 
depth  (as  expressed  by  Pr  >  and  the  collisional  quenching 
rate  (as  expressed  by  .V,),  rather  than  either  factor  alone. 
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determines  the  importance  of  doubiet-opacity  effects  on 
population  densities  and  spectral  emission.  Several  simple 
expressions  obtained  above  express  this  interplay  semi- 
quantitatively  and  convincingly  demonstrate  the  trends  to 
be  anticipated.  We  now  turn  to  specific  numerical  results 
for  further  illumination  of  the  question. 

III.  NUMERICAL  MODEL 
A.  Model  description 

In  order  to  validate  and  clarify  the  analytic  approxima¬ 
tions  presented  above,  as  well  as  provide  more  specific 
spectral  diagnostics,  we  have  performed  an  extensive  series 
of  numerical  calculations  centered  on  one  element  (argon). 
The  calculations  consist  of  a  set  of  simultaneous  solutions 
of  the  coupled  atomic-rate  and  radiative-transfer  equa¬ 
tions  for  argon  plasma  cylinders  where  the  varied  parame¬ 
ters  are  temperature,  density,  and  radius  r,  presumed 
much  smaller  than  the  length.  These  solutions  have  typi¬ 
cally  been  carried  out  with  both  the  more  realistic  doublet 
opacity  and  an  assumed  single  Voigt  profile  for  each  pa¬ 
rameter  set. 

The  atomic  model  consists  of  38  levels  ranging  from 
neutral  Arl  to  the  bare  nucleus  (ArXlX).  Since  our  focus 
is  on  A-shell  radiation,  only  ground  states  are  carried 
through  Be-like  argon  (ArXV).  For  Li-like  Arxvi,  the 
1j22s,  ls22p,  ls23s,  ls23p,  ls23d,  and  Is2,  n=  4  (compos¬ 
ite)  levels  are  modeled.  The  ArXVli  manifold  includes 
Is2 ’5,  ls2s  3S,  U2p  3F,  1s2s  ’S,  ls2p  'P,  the  n  =3  triplets 
and  singlets,  and  n  =4—7  composite  levels.  For  Arxvm, 
n  =  1  — 5  are  included.  The  various  levels  are  connected 
through  electron  collisional  excitation  and  deexcitation, 
collisional  ionization  and  radiative  recombination,  three- 
body  recombination,  and,  where  appropriate,  dielectronic 
recombination.  Photoionization  and  photoexcitation  are 
accounted  for  by  solving  the  radiative  transfer.  All  the 
modeled  ionization  edges  are  included  as  part  of  the 
radiative-transfer  model.  Additionally,  all  the  optically 
thick  resonance  lines  inherent  in  the  level  structure 
described  above  are  also  transported,  along  with  the  2-3 
line  of  Ar Xviii  and  the  \s2p  'P—  ls3d  ’ D  line  of  ArXVli. 
While  some  of  the  higher  levels  have  been  omitted  because 
of  storage  and  efficiency  considerations,  the  model  pro¬ 
vides  a  substantially  complete  description  of  the  processes 
which  form  the  dominant  spectral  lines  of  AT-shell  argon. 
A  discussion  of  the  atomic  rate  calculations  is  contained 
in  Ref.  17.  In  all  results  presented  below,  the  plasma  is 
assumed  to  be  in  collisional-radiative  equilibrium  (CRE); 
namely,  that  each  level-population  density  is  in  a  steady 
state  consistent  with  the  atomic  rates  and  photon  fluxes 
produced  by  the  other  level  populations. 

Some  improvements  have  been  effected  in  our  technique 
for  carrying  out  the  simultaneous,  steady-state  solution  of 
the  radiative-transfer  and  rate  equations. 2-4, 17-22  Refer¬ 
ence  21  describes  a  multicell-coupling  technique  for  radia¬ 
tive  transfer  in  spherical  and  cylindrical  geometries,  based 
on  frequency-integrated  line  escape  probabilities.  This 
method  produces  nearly  exact  solutions  while  employing 
only  one  ray  angle  and  the  computational  equivalent  of 
just  one  frequency  per  line.  To  account  for  doublet  opaci¬ 


ty,  where  escape  probabilities  have  yet  to  be  tabulated  or 
fitted  analytically,  we  have  adopted  a  hybrid  of  the 
method  of  Ref.  21  and  more  conventional  multifrequency 
approaches.  The  single-angle  feature  of  Ref.  21  is  re¬ 
tained  but  each  line  is  divided  into  17—43  frequency 
groups — depending  on  the  line  shape  and  optical  depth. 
The  coupling  constants  are  defined  in  terms  of  individual 
frequencies  (where  the  escape  probability  at  each  frequen¬ 
cy  is  purely  exponential)  rather  than  individual  lines.  This 
also  allows  inclusion  of  continuum  optical  depths  and  line 
and  continuum  cross-pumping.  Given  the  population 
densities  a  set  of  cell-to-cell  photon  coupling  constants  at 
each  of  more  than  350  frequencies  is  thus  obtainable. 

The  iteration  technique  by  which  self-consistency  be¬ 
tween  the  radiative-transfer  and  steady-state  rate  equa¬ 
tions  is  obtained  has  also  been  significantly  improved. 
The  most  straightforward  and  obvious  procedure,  using 
the  radiation  field  from  the  previous  iteration  to  calculate 
populations  which  then  enables  the  radiation  field  to  be 
recomputed  until  consistency  is  achieved,  is  known  as  A 
iteration.  This  technique  is  a  very  poor  choice  for  low- 
density,  high-optical-depth  astrophysical  problems23  since 
the  number  of  iterations  required  is  approximately  equal 
to  the  mean  number  of  photon  scatterings,  which  can  be 
very  large  under  such  circumstances.  However,  the  situa¬ 
tion  is  not  nearly  so  difficult  for  laboratory  plasmas  where 
the  optical  depths  are  smaller  and  the  quenching  probabil¬ 
ities  larger  than  those  which  prevail  in  astrophysical  situa¬ 
tions,  resulting  in  far  fewer  scatterings.  We  employ  the 
mathematical  equivalent  of  Rybicki's  core-saturation 
method24  and  thereby  render  the  A-itcration  technique 
quite  serviceable  for  laboratory  plasma  calculations.  As 
noted  by  Rybicki,24  the  conceptual  and  computational 
simplicity  of  the  technique  is  a  very  desirable  feature, 
especially  for  complex  multilevel  problems. 

To  demonstrate  our  application  of  the  core-saturation 
method  for  A  iteration,  we  consider  one  line  characterized 
by  a  profile-averaged  photon  escape  probability  from  the 
plasma  P,  and  quenching  probability  Pq  Let  S'J 
represent  the  upper  level  density  obtained  on  the  ith  itera¬ 
tion  and  A  the  spontaneous  decay  probability.  The  pho¬ 
toexcitation  rate  for  the  ith  iteration  is  therefore 
jVu'  A  ( 1  —  F, ).  If  C„  and  Du  are  the  collisional  creation 
and  depopulation  rates,  respectively,  for  .V.,  A','""1'  is  ob¬ 
tained  from  the  equation  balancing  gains  and  losses  in  the 
level 

Nl‘  +  '\A+Du)  =  Cu+\u,lAl\-Pt)  ,  (19) 

or 

N,J  +  "=*Nl"  +  \,J\  1  -Pt )( 1  -PQ)  ,  (20) 

where  Pq  =  Du  /( A+D„ )  as  before,  and  the  upper  level 
population  for  an  optically  thin  plasma  is  given  by  the 
first  iterative  solution 


Equation  (20)  yields  a  recursive  sum  for  the  nth  iteration's 
solution 
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n  —  l 

i-0 


(22) 


The  slow  convergence  of  the  sum  of  Eq.  (22)  when 
P„Pq  «  1  precisely  expresses  the  difficulty  of  A  iteration. 
Note,  however,  that  the  sum  is  analytically  evaluable  to 
give 


,V 


(  «  ) 

u 


=  iV 


(1) 

u 


_ 1 _ 

l-(l-pQ)ll-pt) 


(23) 


and  that  this  is  exactly  what  is  obtained  for  A',  on  the  first 
iteration  if,  in  Eq.  (19)  A  is  replaced  by  AP,  instead  of  cal¬ 
culating  excitations  as  Ad  —  P,).  This  is  the  essence  of 
Rybicki’s  technique24  of  eliminating  the  readily  absorbed 
line  core  photons.  In  practice,  we  apply  the  method  by 
calculating,  cell  by  cell,  the  escape  probability  from  the 
plasma  for  each  line  by  integrating  over  the  absorption 
profile — including  the  continuum.  The  spontaneous  emis¬ 
sion  coefficients  are  then  diluted  by  this  factor  on  the  first 
iteration,  which  generally  results  in  an  excellent  approxi¬ 
mation  to  the  final  solution,  especially  near  the  plasma 
center.  On  subsequent  iterations,  the  spontaneous  emis¬ 
sion  coefficients  are  diluted  in  each  cell  by  the  escape 
probability  from  that  cell  and  excitations  are  computed 
explicitly  only  for  the  line  wing  photons  arriving  from 
other  cells.  This  completely  general  procedure  usually 
converges  within  20  iterations  even  for  several  hundred 
Ly-a  optical  depths.  Convergence  has  been  verified  by 
starting  with  entirely  different  initial  iterative  solutions, 
and  by  quadrupling  the  number  of  iterations,  both  of 
which  yield  the  same  final  level  populations. 


B.  Z  scaling  of  numerical  results 

In  order  that  the  results  given  below,  specifically  ob¬ 
tained  for  argon,  may  be  readily  applied  to  other  low-to- 
medium-Z  elements,  we  present  the  methodology  for  Z 
scaling  in  this  section.  Given  a  A-shell  argon  plasma  of 
radius  R0,  ion  density  A"/,  and  temperature  Tt,  the  results 
may  clearly  be  applied  to  a  plasma  of  atomic  number  Z  at 
a  temperature  of  T,(Z/18)2.  However,  scaling  the  density 
and  size  is  not  quite  so  straightforward.  To  maintain  the 
analogous  level  populations  and  spectrum,  it  is  evident 
from  the  discussion  in  Sec.  II  that  both  the  quenching 
probabilities  and  optical  depths  must  remain  the  same. 
The  quenching  probability  is  determined  by  the  ratio  of 
collision  rates — which  scale  as  Z~3 — to  spontaneous  de¬ 
cay  rates — which  scale  as  Z4.  Since  the  quenching  proba¬ 
bility  consequently  scales  as  Z-7,  the  electron  density 
must  scale  as  Z7  to  offset  this.  Therefore,  the  ion  density 
must  increase  as  Z6  to  maintain  the  same  ratio  of  radia¬ 
tive  to  collisional  rates  as  would  prevail  in  an  argon  plas¬ 
ma.  Finally,  to  maintain  the  same  line  optical  depth  rn, 
note  that  it  is  proportional  to  RoN/MM /T,)'n  where  k  is 
the  wavelength,  M  the  ion  mass,  and  Tt  the  ion  tempera¬ 
ture.  Since  Af/~Z4,  k  —  Z~2,  M-Z,  and  Tt~-Z'-. 
r0~R0Z}>.  Therefore,  the  size  must  scale  as  Z-3?  to 
maintain  the  optical  depth.  In  summary,  a  A'-shel!  plasma 
of  atomic  number  Z,  radius  R0d&/Z)3i.  ion  density 
Ar;(Z/18)6,  and  temperature  T,(Z/18)J  will  produce  the 
same  spectrum  as  an  argon  plasma  of  radius  R0 ,  ion  den¬ 


FIG.  1.  Normalized  fine-structure  profiles  of  the  absorption 
coefficients  of  Ar  win  Ly  a,  0.  and  y  are  compared  to  the  sin¬ 
gle  Voigt  profile  which  would  characterize  the  opacity  in  the  ab¬ 
sence  of  fine-siruciure  splitting. 

sity  A’/,  and  electron  temperature  T,.  The  “same"  spec¬ 
trum  is  defined  as  the  same  relative  intensities  of  the  lines. 
Absolute  intensities  will  differ  due  to  the  differing  line 
photon  energies  and  total  numbers  of  emitting  ions.  We 
recommend  Z  scaling  only  for  I3<Z<26  from  the 
Z  =  18  results  presented  below,  since  the  Z'  s  scaling  of 
fine-structure  splitting  [Eq.  (3)]  will  spoil  the  applicability 
of  the  scaling  far  from  Z  =  18. 

C.  Numerical  results 

In  Fig.  I  the  doublet  absorption  profiles  for  Ar  .Will  Ly 
a,  Ly  /?,  and  Ly  y  are  shown  along  with  the  single  Voigt 
profile  which  would  characterize  the  line  opacity  in  the 
absence  of  fine  structure.  The  profiles  are  normalized  for 
case  of  comparison.  Plasma  conditions  of  ion  temperature 
1.4  keV  and  ion  density  3.7X101'*  cm “ '  are  assumed. 
Note  that  the  Ly-a  profile  is  split  into  two  distinct 
features  since  the  doublet  components  are  -  5  Doppler 
widths  apart.  However,  both  Ly  0  and  Ly  y  appear 
quasi-Gaussian,  with  a  skewing  of  the  profile  toward  the 
stronger  of  the  fine-structure  components.  For  Ly  y,  the 
change  from  the  single  Voigt  profile  is  quite  small. 


0  !  - 
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FIG.  2.  Ratio  of  the  n  =2.  3,  and  a  Amin  populations  ob¬ 
tained  with  the  correct  fine-structure  profiles  to  those  obtained 
using  a  single  Voigt  profile  are  plotted  against  ion  density  for 
the  indicated  cylindrical  plasma  temperature  and  diameter 
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Therefore,  rhe  following  results  will  concentrate  on  Ly-a 
and  Ly-f?  power  outputs  and  line  ratios  which  could  re¬ 
flect  those  outputs  as  influenced  by  fine-structure  opacity. 

In  Fig.  2  the  ratio  of  populations  calculated  with  realis¬ 
tic  doublet  opacities  to  those  which  would  prevail  in  the 
absence  of  fine  structure  is  plotted  against  density,  for  an 
argon  plasma  temperature  of  1.4  keV  and  a  diameter  of 
1.34  mm.  The  changes  in  n=  3  and  4  are  15%  or  less, 
whereas,  for  an  ion  density  of  -~3x  10|g  cm'3,  the  n  =2 
population  is  reduced  by  a  factor  of  2.  For  all  three  levels, 
the  population  ratio  approaches  unity  at  both  low  and 
high  densities,  consistent  with  Eq.  (71.  For  the  lsI/;-2pj/: 
component,  the  line-center  optical  depths  are  0.14,  2,  47, 
and  420  at  ion  densities  of  lO1*.  10|q.  I0‘°.  and  10" 1  cm-'. 


respectively.  The  corresponding  optical  depth  of  a  single 
Voigt  profile.  rn,  would  be  50%  greater.  Since  .V,  a  ISA’/, 
tq  =  3.5x  10-;".V,  for  ion  densities  between  1020  and  1021 
cm-3.  The  approximate  escape  probability  for  a  Voigt 
profile13  is  given  as  a  function  of  the  damping  parameter 


a  and  optical  depth  r0  for  r0»  1: 

I  1 1/2 

P,=Z  0.85  —  I 

T0  I 


(24) 


The  damping  parameter  a  for  a  resonance  line  is 

r 

4irAv0 


(25) 


where  the  upper-state  inverse  lifetime  F  for  Ly  a  is  most¬ 
ly  determined  by  radiative  decay  in  this  density  range. 
The  Doppler  width  Av0  at  1.4  keV  is  2.2x  1014  Hz.  In¬ 
cluding  the  collisional  contribution,  a  =2.5 X  10-2  for 
these  densities,  and  Eqs.  (24)  and  (13)  may  be  combined  to 
obtain  the  electron  density  at  which  a  10%  departure  of 
n=  2  population  due  to  fine-structure  opacity  would  be 
expected  according  to  the  analytic  theory 


9.3xlO-16(18r7iV,>(U5 


2.5x  I0~2 
3.5  X  10-20iV, 


1/2 


(26) 


which  yields  an  electron  density  of  6.  IxIO21  cm-3, 
equivalent  to  an  ion  density  of  3.4 x  1020  cm-3  which  is  in 
excellent  agreement,  considering  the  approximations,  with 
the  numerically  calculated  value  of  4.5  x  10:o  (Fig.  2). 

The  most  experimentally  detectable  consequence  of  fine 
structure  is  certainly  spectroscopic  effects.  The  analytic 
results  presented  above  also  contain  [Eq.  (18)]  an 
electron-density  criterion  for  10%  enhancement  of  Ly  a. 
In  Fig.  3  we  present  two  line  ratios  which  are  both  tem¬ 
perature  and  density  sensitive:  He  U2—  ls4p  'P/ Ly  P  and 
He  Is2— ls3p  'P/ Ly  a,  as  a  function  of  density,  with  and 
without  the  correct  doublet  opacity.  The  ratios  decrease 
substantially  with  density  even  with  fixed  temperature  (1.4 
keV!  and  plasma  diameter  (1.34  mm)  since  higher  collision 
rates  and  optical  depths  increase  the  excitation  and  ioniza¬ 
tion  state  of  the  plasma.  However,  at  the  highest  densities 
the  Arxvii  to  Arwin  line  ratios  level  off  and  then  in¬ 
crease.  This  is  primarily  due  to  the  increased  opacity  of 
the  Lyman  lines  coupled  with  the  decreased  opacity  of  the 
heliumlike  lines  as  the  average  charge  state  of  the  plasma 
increases  with  density.  The  double  valuing  of  line  ratios 


FIG.  3.  Line  ratios  Arxvu  |s:-ls3p  'P/Ly  a  and  Arxvii 
Is2— li4p  'P/Ly  f)  are  plotted  against  ion  density  for  a  plasma 
diameter  of  1.34  mm  and  temperature  1.4  keV.  Where  fine- 
structure  opacity  results  in  significant  differences,  it  is  indicated. 
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FIG.  4.  Same  line  ratios  of  Fig.  3  are  p.otted  against  tem¬ 
perature  for  the  same  size  cylindrical  argon  plasma  (1.34  mmi  at 
an  ion  density  of  3.7xI014  cm"',  where  the  presence  of  fine 
structure  makes  no  difference. 
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due  to  opacity  effects  has  been  noted  elsewhere. '“  At  ion 
densities  of  —  1020  and  higher,  the  dichotomy  between 
fine-structure  and  single-profile  calculations  is  noticeable; 
it  reflects  the  higher  power  outputs  of  both  Ly  a  and  Ly  0 
due  to  the  greater  effective  escape  probabilities  of  the 
doublet  profiles.  For  Ly  a,  these  results  provide  another 
point  of  comparison  with  the  theory  of  Sec.  11.  According 
to  Eq.  (18),  the  Ly  -a  line  intensity  should  be  enhanced  by 
10%  when 

5.2  x.  10~  15iV,Z -'>/>,  .  (18) 

Equation  (24),  plus  the  fact  that  r0  =  3.5  x  10_:o.V,  allows 
Eq.  (18)  to  be  solved  for  .V,;  the  result  predicts  a  10% 
power  enhancement  at  .V;  /18=l.!  xI0-’°  cm-3,  in 

excellent  agreement  with  the  numerically  calculated  value 
of  1.3x  lCr°  cm'-  (Fig.  3). 

Figures  4  and  5  present  this  same  ratio  as  a  function  of 
temperature,  again  for  a  plasma  of  fixed  diameter  1.34 
mm.  In  Fig.  4  the  ion  density  is  3.7x10'“  cm'5.  At 
these  relatively  modest  optical  depths  (  —  10  for  Ly  a  I  the 
ratios  are  single  valued  and  a  very  good  temperature  indi¬ 
cator.  Also,  fine  structure  does  not  affect  the  power  out¬ 
put  at  this  density;  therefore,  only  one  curve  appears  for 
each  line  ratio.  This  may  be  contrasted  with  Fig.  5  where 
the  same  information  is  presented  at  the  considerably 
higher  ion  density  of  IO21  cm-'’.  Fine-structure  opacity 
noticeably  depresses  the  line  ratios  for  the  entire  tempera¬ 
ture  range  considered.  Also,  the  much  larger  optical 
depths  virtually  eliminate  the  temperature  sensitivity  of 
the  ratios  above  1.0  keV;  as  the  population  of  Arxvm  in¬ 
creases  with  temperature,  so  does  the  opacity,  which  con¬ 
siderably  reduces  the  tendency  of  the  Lyman  line  to  in¬ 
crease  in  intensity.  At  an  ion  density  of  3.7X  10'“  cm"  ', 
the  Lyman  photons  generally  escape  after  a  few  scatter¬ 
ings;  at  102'  cm"',  they  are  mostly  collisionally  quenched. 

The  fact  that  the  enhancement  of  Ly  /3  is  comparable  to 
that  of  Ly  a  may  seem  puzzling  in  light  of  the  fact  that 
the  Ly-/?  profile  is  not  so  drastically  affected  by  fine 
structure  iFig.  1).  However,  the  Lorentz  wings  of  the  pro¬ 
file  are  substantially  enhanced  by  the  increase  in  damping 
parameter  (i.e.,  reduction  in  level  lifetime)  when  collisions 
transferring  population  among  fine-structure  levels  are 
considered.  At  such  high  opacities,  substantial  numbers 
of  photons  can  escape  only  in  these  far  wings — resulting 
in  the  Ly-/?  power  increase.  For  Ly  a,  as  noted  abov  e,  the 
upper  level  lifetime  is  still  mostly  determined  by  its  very- 
high  radiative  decay  rate,  even  when  collisions  among 
fine-structure  levels  are  considered. 

Finally,  in  Fig.  6  we  present  the  evolution  of  the  emit¬ 
ted  Ly -a  profile  as  a  function  of  density  for  a  temperature 
of  1.4  keV  and  a  plasma  diameter  of  1.34  mm.  At 
A;  =  1018  cm-3,  the  2 p3/2  component  is  twice  the  strength 
of  the  2 pxn  since  the  lines  are  optically  thin  at  this  densi¬ 
ty.  At  6x  1018  cm-3,  r0»  1  and  the  2pi/2  component  be¬ 
gins  to  saturate  to  the  value  of  the  source  function — 
allowing  the  2pi/;  to  begin  to  approach  its  intensity.  At 
3.7x10'“  cm-3,  full  saturation  has  occurred — the  com¬ 
ponents  are  of  equal  intensity  and  a  small  self-reversed 
core  is  seen  in  each.  The  profiles  at  ion  densities  of  1020 
and  3X  1020  cm-3  are  qualitatively  similar  to  that  at 
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FIG  5.  Same  as  Fig.  4,  except  the  density  is  10:’  ions  per 
cm'’.  At  this  density  the  result  depends  upon  the  details  of  the 
fine-structure  profile  and  such  is  indicated. 

3.7x  10‘“  cm-1  except  that  the  higher  optical  depths  pro¬ 
duce  deeper  self-reversed  cores  At  ion  densities  of 
7.5  x  10;o  cm"3  and  I021  cm"  ’,  the  profiles  have  a  pecu¬ 
liar  three-pronged  appearance.  The  three  peaks  appear 
where  the  profile  optical  depth  is  approximately  unity — 
one  peak  on  each  of  the  far  wings  of  the  components,  with 
another  emission  peak  occurring  at  the  absorption 
minimum  between  the  two  components. 

Because  of  this  density  sensitivity,  the  Lv-o  profile  is 
potentially  an  excellent  density  diagnostic.  Since  the  pro- 
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FIG.  6.  Evolution  of  the  emitted  Ly-o  profile  is  presented  3' 
a  function  of  density  for  an  argon  plasma  of  I  4  KeV  and  diam 
eter  1.34  mm.  The  splitting  of  the  iwo  components  visible  at 
the  lowest  density  is  5  4  mA  and  this  horizontal  scale  is  m  un¬ 
tamed  throughout.  The  absolute  intensities  have  been  eorn.i 
ized  for  ease  of  comparison 
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files  depend  on  optical  depth  as  well  as  density,  the  size  of 
the  plasma  must  be  ascertained  to  fully  exploit  this  tech¬ 
nique.  Experimental  methods  such  as  x-ray  pinhole  pho¬ 
tography  would  be  quite  valuable  for  this  purpose.  Also, 
spectral  resolution  of  —1—2  mA  will  be  required,  in¬ 
dependently  of  Z. 
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B.  Ionization  Dynamics  and  Diagnostics 


I.  INTRODUCTION 


In  response  to  the  new  experiments  being  planned  on  the  Gamble  II 
machine  during  FY84,  several  new  atomic  models  have  been  constructed  and 
implemented  in  the  ionization/radiation  models.  With  the  acquisition  of 
the  Hartree-Slater  atomic  structure  code  of  Cowan,  the  compilation  of 
energy  levels  and  radiative  decay  rates  has  been  vastly  improved,  both  in 
the  increased  comprehensiveness  of  the  data  as  well  as  the  enhanced 
accuracy.  This  improvement  is  directly  reflected  in  the  predictive 
capability  of  the  radiation  models. 

II.  MODELS 

Two  new  K-shell  models  were  completed  this  year,  for  silicon  and 
titanium.  Both  models  contain  the  sufficient  level  structure  needed  to 
accurately  diagnose  x-ray  spectra  from  dense  plasmas  of  these  materials. 
In  particular,  this  Includes  up  to  n-values  (principal  quantum  number)  of  5 
in  the  hydrogenlike,  heliumlike,  and  lithiumlike  ions.  The  silicon  models 
expected  to  be  employed  in  imploding  glass  wire  experiments  on  the  pulse 
power  generator  at  NRL.  Conceivably,  glass  could  be  incorporated  as  a  core 
element  or  a  mixture  with  aluminum  in  order  to  test  the  cooperative  photon 
pumping  effects  between  the  two  materials.  An  example  of  a  silicon  K-shell 
spectrum  is  shown  in  Fig.  1. 

The  titanium  K-shell  model  was  likewise  constructed  to  diagnose 
imploding  wire  arrays  of  this  material.  During  previous  experiments  on  the 
BJV  generator,  titanium  received  serious  consideration  as  a  PRS  load 
material.  With  the  upgrades  on  existing  machines,  It  is  possible  that 
titanium  wires  will  be  re-examined  as  a  viable  radiation  source.  A  typical 
titanium  K-shell  spectrum  is  shown  in  Fig.  2. 

It  has  often  been  conjectured  that  it  may  be  far  easier  to  use  the  L- 
shell  x-ray  emission  from  a  high-Z  material  rather  than  the  K-shell  x-rays 
from  a  lower-Z  material  (e.g.,  neon  aluminum,  argon)  to  obtain  radiation 
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emission  of  the  right  color  for  the  DNA-PRS  goals.  To  this  end, 
experimentalists  at  PI  have  been  working  with  a  krypton  gas  puff  during 
recent  months.  Krypton  has  also  been  selected  as  the  primary  load  material 
to  test  the  new  gas  puff  capabilities  presently  being  installed  on  the  NRL 
Gamble  II  generator.  In  order  to  provide  theoretical  support  for  these 
programs,  an  L-shell  model  for  krypton  has  been  constructed.  This  model 
represents  a  Herculean  effort  in  both  data  compilation  and  calculation  and 
is  the  most  advanced  atomic  structure  modeling  tool  to  date.  In  the 
construction  of  this  model,  over  15,000  oscillator  strengths  were 
calculated;  these  were  reduced,  via  statistical  averaging,  to  over  450 
lines  from  over  100  different  excited  levels  in  11  ions  from  Na-like  Kr 
XXVI  to  H-like  Kr  XXXVI.  Nearly  600  electron  impact  excitation  rates,  150 
collisional  ionization  rates  and  photoionization  cross-sections,  and  15 
effective  dielectronic  recombination  rates  were  necessary  to  complete  a 
model  suitable  for  use  in  our  collisional-radiative  equilibrium  (CRE) 
ionization/radiation  code.  The  atomic  energy  levels  used  in  the  krypton 
model  are  listed  in  Table  I  along  with  their  statistical  weights;  as  is 
evident,  significant  level  averaging  has  been  done  to  account  for  the  bulk 
of  the  radiation  while  keeping  only  12  excited  states  per  ionization  stage. 

As  an  illustration  of  how  the  predictions  of  various  approximations 
compare  in  dense  plasmas,  our  CRE  model  has  been  used  to  calculate  the 
values  of  various  Ionization  and  radiation  quantities  in  a  dense  (1020 
ions/cm  )  krypton  plasma.  In  Fig.  3,  the  effective  ionization  state  of  the 
plasma  is  plotted  versus  electron  temperature,  for  three  different  models: 
CRE,  LTE,  and  Corona.  Also  shown  is  the  curve  for  CRE  with  opacity  effects 
included  for  an  optically  thick  cylindrical  plasma  (radius  »  500  urn).  For 
this  density,  the  plasma  is  nearly  in  LTE  at  low  temperatures,  while  the 
corona  model  gives  accurate  results  at  high  temperature.  As  the  density 
increases,  the  CRE  results  will  tend  closer  to  the  LTE  curve.  The 
fractional  ion  abundances  (including  excited  state  densities)  are  shown  in 
Fig.  4  as  a  function  of  temperature.  Only  the  L  and  K-shell  ions  are 
included  as  well  as  the  Mg  and  Na-  like  ions  since  no  excited-level 
structure  was  included  in  the  model  for  Kr  I  -  Kr  XXIV.  Also  shown  is  the 
result  for  the  Ne-like  ion  as  predicted  by  the  corona  model.  Notice  both 
the  peak  abundance  as  well  as  the  temperature  at  which  it  occurs  are  in 
disagreement.  This  is  an  important  result  for  x-ray  laser  studies  of 
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population  inversions  in  Ne-like  systems,  since  the  gain  prediction  for  an 
experimental  plasma  configuration  is  a  strong  function  of  the  ground  state 
density;  an  error  from  1.0  to  1.8  keV  in  the  temperature  corresponding  to 
peak  Ne-like  abundance  will  have  a  significant  effect  on  this  prediction. 

However,  the  peak-abundance  temperature  is  not  necessarily  the  value 
one  wants  to  attain  experimentally.  In  Fig.  5,  the  line  intensities  for 
several  Ne-like  transitions  are  shown  versus  temperature  for  a  krypton 
plasma  at  10  ions/cm  .  .  Note  that  although  the  total  ion  population 
density  peaks  at  about  1.0  keV  (see  Fig.  4),  the  line  intensities  peak  at 
nearly  1.3  keV.  This  is  due  to  the  fact  that  the  line  intensities  are 
nearly  proportional  (exactly  proportional  in  coronal  equilibrium)  to  the 
product  of  the  ground  state  population  and  the  electron  impact  excitation 
rate  of  the  upper  level  of  the  bound-bound  transition.  Thus,  if  a 
population  inversion  can  be  attained  in  the  plasma,  the  electron 

temperature  at  which  maximum  gain  can  be  realized  will  probably  be 

somewhere  between  the  values  predicted  by  Figs.  4  and  5. 

The  gross  radiation  characteristics  of  the  krypton  plasma  are  shown 
,  for  the  various  models  in  Fig.  6  as  a  function  of  temperature.  The  CRE 
"thick"  results  are  also  shown  for  the  same  geometric  conditions  as  stated 
earlier.  Neither  the  LTE  nor  the  coronal  model  predict  the  functional  form 
of  the  radlational  cooling  very  realistically,  until  30-50  keV 
temperatures,  where  bremsstrahlung  is  the  dominant  radiative  mechanism. 
The  LTE  predictions  are  by  far  the  worst,  calculating  an  intense  peak  in 
the  K-shell  radiation  at  the  same  temperatures  that  CRE  results  indicate 
the  L-shell  radiation  is  peaking.  The  corona  model,  on  the  other  hand, 

predicts  an  incorrectly  large  M-shell  peak  at  about  0.6  keV,  since 

collisional  quenching  is  neglected  in  the  coronal  approximation.  Note  that 
opacity  effects  tend  to  reduce  the  radiative  cooling  by  about  a  factor  of 
three  for  these  conditions  at  plasma  temperatures  corresponding  to  peak  L- 
shell  emission  ( ~ 2 . 5  keV) . 

Finally,  a  breakdown  of  the  various  radiating  mechanisms  is  shown  in 
Fig  7,  for  the  same  conditions  as  Fig.  6.  The  contributions  due  to  bound- 
bound  (lines),  bound-free  (recombination)  and  free-free  (bremsstrahlung) 
radiation  are  illustrated,  as  well  as  the  K-line  radiation  as  a  function  of 
electron  temperature.  At  low  temperatures,  line  emission  is  totally 
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dominant,  while  bremsstrahlung  takes  over  at  about  30  keV.  The  radiative 
recombination,  although  clearly  not  a  factor  in  the  total  radiation  at 
1020  ions/cm2,  can  become  dominant  as  the  density  increases.  However,  the 
density  at  which  this  occurs  is  an  increasing  function  of  plasma  material 
Z.  Thus,  in  laser-heated  targets,  for  example,  bound-free  processes  may 
dominate  the  radiational  cooling  for  carbon  and  aluminum  while  materials 
like  molybdenum  and  gold  radiate  mainly  via  line  emission. 

As  an  illustration  of  a  typical  krypton  emission  spectrum,  the 

emission  from  a  puff  plasma  with  a  linear  temperature  gradient  (1. 0-8.0 
keV)  is  shown  in  Figs.  8a-8d. 

Interesting  features  of  the  spectrum  include  line  emission  gap  from 
1.0-1. 6  keV  (Fig.  8b);  the  An=0  L-shell  transitions  lie  below  the  gap  while 
the  An=l,2,...  transitions  all  lie  above  it.  The  M-shell  lines  (not 
included  in  this  model)  all  lie  below  900  eV  and,  hence,  will  not  fill  the 
empty  region  with  line  emission.  Above  2.8  keV,  the  line  emission  ceases 
once  more  (Fig.  8c)  until  finally  the  K-shell  Lyman  series  lines  appear  at 
13.0  keV  (Fig.  8d) . 

This  model  can  be  used  both  as  a  diagnostic  post-processor  (in 

conjunction  with  our  CRE  model)  or  as  an  on-line  radiation  routine,  used 
with  our  l-D  MHD  gas  puff  implosion  code.  Presently,  we  are  attempting  to 
reduce  the  size  of  this  model  significantly  for  more  efficient  and 
economical  implementation  in  the  hydro  codes,  while  yet  maintaining  the 

bulk  radiation  energetic  properties  of  the  model  to  high  accuracy.  It  will 
be  possible  to  insure  this  accuracy  by  detailed  comparison  between  the 
reduced  and  full-blown  versions  of  the  krypton  model. 

III.  DIELECTR0NIC  SATELLITE  LINE  MODEL 

The  main  emphasis  of  the  program  in  plasma  diagnostics  development 
this  year  centered  on  the  construction  of  a  model  for  dielectronic 
satellite  lines.  Satellite  lines  are  emitted  as  a  consequence  of  the 

radiative  decay  of  doubly-excited  levels  and  lie  close  to  main  resonance 
lines  in  the  emission  spectrum.  The  importance  of  these  spectral  features 
to  the  DNA-PRS  program  lies  in  their  value  as  plasma  temperature  and 
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density  indicators  coupled  with  their  prominence  in  the  spectrum  of  gas 
puff  and  wire  array  implosions;  the  lines  of  special  interest  to  the 
program  are  the  lithiumlike  and  heliumlike  satellites. 

One  of  the  main  advantages  of  using  satellite  lines  as  plasma 
diagnostics  is  that  both  they  and  the  resonance  lines  are  populated  by 
electron  collisions  with  the  same  ground  state.  This  removes,  somewhat, 
the  uncertainty  of  the  ionization  balance  from  the  density  and  temperature 
determinations.  The  preliminary  study  of  these  features .describing  the 
construction  of  the  model  and  the  temperature  and  density  variation  of  the 
satellite  line  intensities  have  been  documented  in  a  detailed  publication; 
a  reprint  of  this  paper  can  be  found  in  the  Appendix  B  of  this  report. 
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Table  1  -  Atomic  energy  levels  and 
Kr  I  -  Kr  XXV 

ground  states  only 

Kr  XXVI  (Na-jjOte) 


ls22s22p63s 

2 

3p 

6 

3d 

10 

4s 

2 

4p 

6 

4d 

10 

4f 

14 

51 

50 

6 1 

72 

Kr  XXVII 

(Ne-like) 

statistical  weights  for  krypton 


Kr  XXVIII 

(F-like) 

22s22p5  (2P) 

6 

2s2p6  (2S) 

2 

2s22p43s 

30 

2s22p43p 

90 

2s22p43d 

150 

2s2p53s 

24 

2s2p53p 

72 

2s2p53d 

120 

2s22p'44f. 

480 

2s2p^4?, 

2s22p45;. 

384 

2s2p35(l. 

2s22p46?, 

1350 

2s2p561 

1944 

ls22s22p6  ( 

lS) 

1 

2s22P53s 

(3/2) 

8 

2s22p53s 

(1/2) 

4 

2s22p53p 

(3/2) 

24 

2s22p53p 

(1/2) 

12 

2s22p53d 

(3/2) 

40 

2s22p53d 

(1/2) 

20 

2s2p63s 

4 

2s2p63p 

12 

2s2p63d 

20 

2s22p54.t 

192 

2s22p557. 

300 

2s22p56£ 

432 

Kr  XXIX  (O-like) 


ls22s22p4  (3P2)(1Sq) 

6 

2s22p4  (3P1)(1D2) 

8 

2s 22p4  (3Pq) 

1 

2s2p5  (3p) 

9 

2s2p5  (XPX) 

3 

2p6  (1So) 

1 

2s22p33l 

360 

2s2p33?. 

270 

2s22p34f. 

640 

2s2p44c, 

480 
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Table  1  (Cont'd)  -  Atomic  energy  levels  and  statistical  weights  for  krypton 


Kr  XXX  ( >1-1  ike) 

Kr  XXXII  (B-like) 

ls32s22p3  (4S3/2) 

4 

ls22s22p  (2P1/2) 

2 

2s22p3  (2D)(2Pl/2) 

12 

2s22p  (2P3/2) 

4 

2s22p3  (2P3/2) 

4 

2s2p2  (4P) 

12 

2s2p4  (4P) 

12 

2s2pz  ( “D) ( 4P1/2) 

12 

2s2p4  (2D) 

10 

2s2p2  (2S)(2P3/2) 

6 

2s2p4  (231/2K2P3/2) 

6 

2p3  (4S) 

4 

2s2p4  (2P1/2) 

2 

2pJ  (-D)(-P1/2) 

12 

2p5  (2P) 

6 

2P3  (2P3/2) 

4 

2s22p23.1 

270 

2s231 

18 

2s2p33.2 

720 

2s2p( 1/2)32 

82 

2s22p245. 

480 

2s2p( 3/2)30 

134 

2s2p34?. 

1280 

2s  24?. 

32 

2s2p4!>. 

384 

Kr  XXXI  (C-like) 

Kr  XXXIII  (Be-like) 

ls22s22p2  (3Pq) 

1 

2s22p2  (3P1)(1D2) 

8 

ls22s2  (1Sq) 

1 

2s22p2  (3P2)(1SQ) 

6 

2s2p  (3P) 

9 

2s2p3  (5S2) 

5 

2s2p  ( 1?1) 

3 

2s2p3  (3D) 

15 

2P2  (3P0) 

1 

2s2p3(3P0)1)(1D)(3S) 

12 

2p2  (3P1)(1D2) 

8 

2s2p3  (3P2)(1P1) 

8 

2p2  (3P2)(1S0) 

6 

2P4  <3po,2> 

6 

2s3s/ 3p 

16 

2p4  (3P1)(1n2) 

3 

2s  3d 

20 

2p4  (3S0) 

1 

2p( 1/2)37 

44 

2s22p3C. 

108 

2p( 3/2)3? 

64 

2s2p23d 

540 

2s4f 

36 

2s22p4<> 

192 

2p4?. 

192 

2s2p24f. 

960 

22 
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Table  1  (Cont'd)  -  Atomic  energy  levels  and  statistical  weights  for  krypton 


Kr  XXXIV  (Li-like) 


ls22s 

2 

2p 

6 

3s 

2 

3p 

6 

3d 

10 

4  5 

32 

51 

50 

Kr  XXXV 

(He-like) 

Is2  ( lS)  1 

ls2s  (3S) 

3 

ls2s  (3S) 

1 

ls2p  (3P) 

• 

9 

ls2p  (3P) 

3 

ls3f,  (3L) 

27 

Is  3)1  (lL) 

9 

Is4£ 

64 

ls5A 

100 

Kr  XXXVI 

(H-like) 

Is 

2 

21 

8 

51 

18 

bl 

32 

51 

50 

23 
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1  Theoretical  CRE  silicon  K-shell  spectrum  at  T  =1.0  keV 


Fig.  4  Fractional  ion  abundances  (CRE)  of  an  optically  thin  krypton 

plasma  at  lCru  ions/cra  .  Also  shown  is  the  Ne-like  abundance  as 
predicted  by  a  corona  model. 
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Fig.  7  The  CRE  radiative  cooling  curve  for  krypton,  as  in  Fig.  6,  showing 
the  contributions  by  line,  recombination  and  bremsstrahlung 
emission,  also  shown  in  the  K-line  radiation. 
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C.  Detailed  Modeling  of  Ionization/Radiation  in  Dense  Plasmas 


I.  INTRODUCTION 

In  this  section,  some  of  the  problems  involved  in  modeling  plasma 
radiation  via  the  detailed  configuration  accounting  method  are  discussed. 
This  consists  of  a  collection  of  ideas,  suggestions  and  difficulties 
encountered  over  several  years  of  work  on  the  DNA-PRS  program  by  several 
members  of  the  NRL  Plasma  Radiation  Branch.  The  main  thrust  of  this 
discussion  is  an  illustration  of  why  such  detailed  methods  are  necessarily 
brought  to  bear  on  DN'A  radiation  programs  and  why  more  approximate  modeling 
techniques  yield  unacceptably  incorrect  predictions.  In  order  to  make  this 
point  as  clearly  as  possible,  the  discussion  will  be  maintained  in  a 
totally  qualitative  context. 

As  a  preview  of  how  the  calculations  from  which  the  results  were 
gleaned  are  performed,  a  short  description  of  the  modeling  techniques  is 
required.  Briefly,  a  set  of  atomic  rate  equations  is  solved  for  the 
population  densities  in  a  plasma  as  a  function  of  ion  density,  electron  and 
ion  temperature,  and  geometric  configuration.  Each  equation  describes 
collisional  processes  affecting  a  specific  atomic  level  and  includes  rate 
coefficients  from  all  important  atomic  processes.  In  addition  to  particle 
collisions,  radiation  transport  is  also  included  in  the  model  using  either 
photon  probabilistic  techniques  or  a  multi-frequency  ray-trace  solution  of 
the  equation  of  transfer,  solved  in  variable  1-dimensional  geometry.  The 
basis  for  the  accuracy  of  such  calculations  is  the  quality  of  the  cross- 
sections  employed.  Presently  we  use  a  mixture  cf  state-of-the  art 
calculations  and  simple  hydrogenic  approximations  to  generate  the  atomic 
data;  a  more  in-depth  discussion  of  data  compilation,  however,  will  be 
presented  below. 

The  discussion  has  been  divided  into  three  areas: 

(1)  Atomic  Data  Compilation 

(2)  Model  Construction  and  Numerical  Constraints 

(3)  Model  Implementation  and  Application. 
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ATOMIC  DATA  COMPILATION1 


A.  Level  Structure 

The  first  stage  in  the  radiation  modeling  is  usually  the  calculation  or 
compilation  of  the  necessary  data  base.  Before  this  can  be  accomplished,  a 
suitable  or  appropriate  level  structure  must  be  selected.  The  selection 
involves  a  decision  to  limit  (i)  the  number  of  ions  to  model  in  detail  and 
(ii)  the  number  of  excited  levels  to  include  in  these  ions.  The  former 
depends  on  the  temperature  and  density  regime  to  be  modeled.  To  analyze 
soft  x-ray  spectra  from  a  laser  microballoon  implosion  with  argon-doped 
fuel^,  for  example,  the  two  K-shell  ions  and  the  lithiumlike  system  (which 
can  effect  the  x-ray  diagnostics,  especially  satellite  lines)  are 
sufficient.  A  radiation-hydrodynamic  treatment  of  a  gas-puff  implosion  or 

2 

a  laser-heated  foil  ,  on  the  other  hand,  requires  the  modeling  of  excited 
levels  in  many  ions  to  adequately  treat  the  radiation  over  the  wide  range 
of  electron  temperatures  encountered  in  these  experiments.  Since  this  is  a 
rather  obvious  point,  let  me  proceed  directly  to  the  second  (and  much  more 
subtle)  point.  One  wants  to  limit  the  number  of  excited  levels  per  ion  for 
two  reasons: 

(1)  the  number  of  cross-sections  required  to  be  calculated  increases 
by  almost  n  (where  n  Is  the  principal  quantum  number)  as  more  excited 
states  are  included,  and, 

(2)  the  required  computation  time  for  radiation  modeling  increases  as 
the  number  of  bound-bound  transitions  (hence,  levels)  increases, 
particularly  due  to  radiation  transport. 

Presented  below  are  two  examples  of  how  the  number  of  excited  levels 
included  in  a  model  can  disturb  the  population  densities.  In  Fig.  1,  the 
fractional  ion  abundances  (ground  plus  excited  states)  of  an  aluminum 
plasma  at  10  ions/cm  versus  electron  temperature  is  shown.  The 
calculation  Is  for  a  plasma  in  collisional-rad iat ive  equilibrium  (CRE)  in 
an  "optically  thin"  plasma  environment  (radiation  transport  is  not 
included).  Notice  the  variation  in  the  lithiumlike  A1  XI  abundance  as 
fewer  excited  levels  are  included  in  the  level  structure.  One  is  spared 
the  endless  task  of  including  higher  n-values  only  by  the  lowering  of  the 


ionization  limit,  which,  at  this  density,  lies  just  above  the  n=6  levels 

3  ^23 

(using  the  interparticle-distance  criteria)  .  In  fact,  at  10  ions/cm  , 

Fig.  1  illustrates  one  of  the  worst  cases,  since  the  effect  of  excited 
states  on  the  ionization  balance  decreases  as  the  density  decreases. 

However,  a  careful  determination  of  the  lowering  is  obviously  necessary  in 
order  to  accurately  predict  the  required  number  of  levels  and  the 
ionization  balance. 

While  this  large  effect  on  the  ion  abundances  is  present  only  at 

higher  densities,  we  have  discovered,  very  recently,  that  limiting  the 
level  structure  in  an  ion  can  also  affect  the  diagnostics,  even  at  lower 
densities.  As  a  consequence  of  our  tracer-dot  spectroscopy  program  at 

N'RL,^’^  it  was  found  that  K-shell  x-ray  diagnostics  can  be  misinterpreted 
if  adequate  structure  is  not  included  in  the  model.  As  an  example.  Fig.  2 
contains  a  plot  of  several  hydrogenlike  to  heliumlike  resonance  line 
intensity  ratios  for  a  cylindrical  aluminum  plasma  typical  of  that  produced 
by  laser  heating  of  a  115  um  tracer-dot  implant.  The  ion  density  is  fixed 

at  lO"  ions/cm  and  the  plasma  is  assumed  to  be  in  CRC.  Three 

calculations  were  performed  with  different  excited  level  structure:  (1) 

the  solid  curves  represent  n<  10  for  both  K-shell  ions,  (2)  the  dashed- 

dotted  curves  represent  n<  10  for  the  heliumlike  system  and  n<  5  for  the 

hydrogenic  ion,  and  (3)  the  dashed  curves  are  with  n<  5  for  both  K-shell 
ions.  As  can  be  seen  from  the  curves,  all  ratios  are  affected  by  limited 
level  structure  in  the  heliumlike  system,  while  only  higher  Rydberg  members 
of  the  hydrogenic  ion  are  affected  by  limiting  the  levels  in  that  system. 
The  error  in  the  temperature  determination  is  50-200  eV  at  this  density, 
depending  on  which  line  ratio  is  used.  Of  course,  the  effect  is 

exacerbated  at  higher  ion  densities  but,  again,  is  limited  by  the  cutoff  in 
n  due  to  lowering  of  the  ionization  limit.  The  message  of  Figs.  1  and  2, 
however,  i*»  clear:  care  must  be  taken  in  selecting  level  structure 

adequate  to  describe  the  main  population  densities  in  the  density  regime 
the  model  is  to  be  applied,  particularly  if  one's  intention  is  to  analyze 
experimental  spectra. 
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B.  Large  Data-Base  Generation 


One  major  obstacle  that  must  be  overcome  when  the  modeling  of  a  large 
number  of  ions  is  required  (as  usually  occurs  where  combined 
radiation/hydrodynamic  simulations  are  performed)  is  the  assembly  of  large 
numbers  of  cross-sections  and/or  rate  coefficients.  This  can  be  a  problem 
not  only  in  the  computational  expense  that  must  be  incurred  but  also  simply 
due  to  the  large  amount  of  time  invested  in  amassing  a  large  data  base. 
For  example,  the  kryton  model  just  completed  contains  all  37  ground  states, 
typically  10  excited  levels  in  each  of  11  ions  (Na-like  through  H-like), 
over  500  oscillator  strengths,  nearly  150  collisional  ionization  rates,  150 
photoionization  cross-sections,  12  effective  dielectronic  recombination 

rates  and  over  600  electron  impact  excitation  rates.  To  construct  such 

models  efficiently  requires  rapid,  inexpensive  mathematical  expressions, 

6  7 

such  as:  the  exchange  classical  impact  parameter  (ECIP)  or  Lotz  formula 

for  impact  ionization,  the  hydrogenic  Seaton  formula  for  photo-ionization, 

Q 

the  Bates-Damgaard°  formula  for  oscillator  strengths,  a  semi-classical 

impact  parameter  (SCIP)^  calculation  for  impact  excitation,  and  the 
Burgess-Merts 10  formula  for  dielectronic  recombination.  From  these  rates 
or  cross-sections  one  calculates  collisional  and  radiative  recombination 
and  electron  collisional  de-excitation  via  detail  balancing.  Although 
these  (and  other)  methods  reduce  the  expense  and  time  requirements 

significantly,  the  quality  of  the  coefficients  is  questionable  in  several 
cases.  This  leads  us  to  the  discussion  in  the  next  section. 

C.  Quality  of  the  Rate  Coefficients 

Of  course,  we  are  all  aware  of  the  effect  that  the  fusion  effort  and 
scientific  space  program  have  had  on  the  atomic  physics  community  over  the 
last  two  decades.  The  need  for  comprehensive  sets  of  accurate  collisional 
data  for  higher  Z  materials  and  higher  charge  states  has  stimulated  many 
workers  to  develop  highly  sophisticated  theoretical  techniques  for 
calculating  these  cross-sections.  This  allows  model  builders  to  be  more 
selective  in  making  a  "wish-list”  for  the  type  of  rate  coefficients  they 
would  like  to  include  in  the  models.  For  example,  with  a  Hartree-Slater 
atomic  structure  code  (such  as  those  of  Cowan  or  Scofield)  one  can  obtain 
highly  accurate  bound  and  continuum  wave  functions,  energy  level 
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eigenvalues  and  oscillator  strengths.  With  these  wave  functions  and  the 
proper  potentials,  one  can  also  generate  quality  photoionization  cross- 
sections.  Also,  using  these  wave  functions,  autoionization  rates  are 
obtained  for  use  in  calculating  accurate  ionization  cross-sections  via  a 
Coulomb-Born  Distorted  Wave  calculation  (with  the  effects  of  excitation- 
autoionization  included,  of  course).  Similarly,  excitation  rates  can  be 
computed  via  either  the  distorted  wave  or  close-coupling  techniques, 
including  the  relevant  autoionizing  resonances.  One  now  has  a  highly 
accurate  set  of  computer  codes  for  calculating  everything  needed  for  a 
sophisticated  atomic  model;  however,  the  expense  and  time  investment 
required  to  exercise  them  is  still  prohibitive.  Eventually,  atomic 
physicists  will  no  doubt  be  persuaded  to  use  their  expertise,  coupled  with 
their  sophisticated  structure  and  scattering  codes,  to  parameterize  large 
numbers  of  coefficients,  perhaps  along  iso-electronic  sequences  in  a  way 
that  can  be  easily  adapted  or  implemented  by  the  model  builders.  Some  of 
that  has  already  been  accomplished  by  Sampson^1  and  co-workers  at  Penn. 
State  U.  with  impact  excitation  collision-strengths,  but  more  work  along 
this  line  is  sorely  needed  if  comprehensive  sets  of  accurate  rate 
coefficients  are  to  be  made  easily  accessible  to  the  community  of  workers 
who  need  them. 

III.  MODEL  CONSTRUCTION  AND  NUMERICAL  CONSTRAINTS 

A.  Diagnostics  or  Energetics? 

The  question  posed  by  the  preceding  subtitle  may  be  puzzling  to  some  of 
you,  but  to  radiation  modelers  it's  an  all  too  familiar  one.  Often,  a 
choice  has  to  be  made  during  the  construction  of  a  model  as  to  whether  it 
will  be  used  in  conjunction  with  a  hydrodynamic  model  for  predicting 
radiation  energetics  or  as  simply  a  stand-alone  tool  for  post-processing  or 
diagnosing  experimental  data.  While  these  applications  are  not  necessarily 
at  cross-purposes  with  one  another,  a  model  is  usually  designed  exclusively 
for  only  one  application  for  reasons  discussed  earlier  in  the  section  on 
limiting  the  number  of  levels.  However,  in  many  cases  where  an  energy 
level  is  required  for  both  energetic  and  diagnostic  purposes,  one  can  still 
reduce  computational  time  by  cleverly  averaging  over  a  group  of  levels 
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which  lie  close  in  energy.  As  an  example,  there  are  eight  n?  levels  in 
the  n  *  4  state  of  heliumlike  ions:  the  4s,4p,4d,  and  4f  states  in  the 
single  and  triplet  systems.  In  a  dense  plasma,  however,  one  can  often  make 

1 0 

the  safe  assumption  that  electron-ion  or  ion-ion  collisions  will  maintain 
a  Boltzmann  statistical  equilibrium  among  the  n?.  components.  Thus,  by 
taking  appropriate  weighted  averages  of  the  various  rate  coefficients 
coupling  the  4i  levels  to  surrounding  states,  the  eight  n  -  4  levels  can  be 
reduced  to  one  with  little  loss  of  accuracy  in  calculating  the  radiation 
from  decay  of  the  n  «  4  state.  Of  course,  what  is  sacrificed  is  the  finer 
details  of  the  spectral  emission,  since,  for  example,  a  single  3d-4  line 
will  be  calculated  whereas  both  a  3d-4p  and  3d-4f  will  be  detected  in  the 
experimental  spectrum.  In  this  way,  one  sacrifices  more  detailed  spectral 
calculations  for  reduced  memory  storage  and  computation  time,  both  of  which 
can  impact  the  usefulness  of  large  radiation-hydrodynamic  calculations.  At 
NRL,  the  solution  to  this  problem  is  currently  being  pursued  by  way  of  a 
two-pronged  attack.  The  full  models  with  more-detailed  level  structure  are 
constructed  first.  These  are  then  used  in  experimental  analysis  (e.g., 
tracer-dot  spectroscopy  and  Z-pinch  diagnosis)  or  in  post-processing  hydro 
code  plasma  profiles  (e.g.,  x-ray  laser  studies).  Then,  extensive 
averaging  of  levels  and  rates  is  done  to  construct  a  much-reduced  model 
suitable  for  hydro  code  applications,  with  the  advantage  that  the 
predictive  capabilities  of  the  reduced  model  can  be  checked  and  benchmarked 
with  the  more  complete  version  for  accuracy.  Using  this  technique,  we  are 
currently  reducing  our  500-level  krypton  model  to  a  66-level  version  for 
implementation  in  our  MHD  codes.  The  differences  between  these  models  will 
be  discussed  in  depth  in  a  report  to  appear  later. 

B.  Short-Cuts  to  Population  Densities 

As  was  stated  earlier,  our  CRE  radiation/ionization  models  are  based  on 
the  solution  of  coupled  atomic  rate  equations  to  calculate  the  various 
level  populations.  Because  of  the  complexity  of  the  method  and  the  large 
data  base  requirements,  attempts  are  often  made  to  use  simpler  but  more 
approximate  alternatives  to  determine  the  state  densities.  Considering 
only  schemes  which  still  maintain  detailed  configuration  accounting  (as 
opposed  to  average-atom  models),  the  methods  fall  into  three  basic 
categories:  (1)  local  thermodynamic  equilibrium  (UTE)  models,  which 
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solve  the  Saha  and  Boltzmann  equilibrium  relations  to  get  population 

densities,  and  are  applicable  in  the  high  plasma  density  limit,  (2) 

corona  models  which  solve  a  set  of  equilibrium  balance  equations, 

neglecting  collisions  when  in  competition  with  radiative  processes,  and  are 

valid  in  the  low  density  limit,  (3)  a  group  of  “approximately  CRE"  models 

which,  like  CRE,  are  valid  over  a  wide  range  of  plasma  densities,  but 

employ  some  special  mathematical  techniques  to  reduce  the  necessary  atomic 

data  or  computation  time.  Examples  of  these  would  be  the  method  based  on 
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population  probabilities  by  Salzmann  and  the  "mixed  model"  of 

Busquet.^  In  general,  all  of  these  models  can  give  rather  accurate 
results  for  the  population  densities  when  applied  to  the  plasma  regimes  for 
which  they  were  designed,  and  it  is  not  our  intent  to  detract  from  their 
value  in  determining  the  level  distributions.  However,  we  have  found  that 
using  short-cuts  to  determine  the  population  densities  does  little  to 
reduce  the  labor  or  cost  of  accurately  predicting  the  radiation  output. 
There  are  two  fundamental  reasons  for  this.  First,  one  still  requires  all 
the  oscillator  strengths  to  make  accurate  predictions  of  line  emission,  as 
well  as  all  the  radiative  recombination  rate  coefficients  to  accurately 
calculate  the  bound-free  emission.  Second,  the  major  computational  cost  of 
performing  radiation  calculations  lies  in  the  radiation  transport 
algorithm,  not  in  the  ionization  balance  calculation.  Hence,  one  avoids 
neither  the  toil  of  compiling  the  atomic  data  nor  the  necessary  cost  of 
computation  by  using  more  approximate  models,  if  the  radiation  field  is  a 
required  result.  The  inevitable  conclusion  we  have  arrived  at,  therefore, 
is  that  one  may  as  well  construct  and  implement  the  complete  CRE  model. 

C.  Radiation  Transport 

As  well  as  calculating  the  collisional  effects  on  the  ionization 
balance,  reabsorption  and  transport  of  radiation  must  also  be  taken  into 
account  in  these  dense  plasmas.  Since  the  transport  calculation  usually 
accounts  for  a  major  portion  of  the  computational  time,  it  is  crucial  that 
the  method  be  made  as  efficient  and  economical  as  possible.  Our  group  at 
NRL  uses  two  different  techniques  to  transport  radiation  in 
ionization/radiation  models:  (1)  probabilistic  photon  reabsorption,  and 

(2)  frequency-by-frequency  ray-trace  solution  of  the  transfer  equation. 
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The  details  of  both  of  these  methods  have  been  published  previously"’  ~ 
and  so  I  will  spare  you  complicated  descriptions  and  formulas.  Basically, 
the  probabilistic  technique  treats  individual  spectral  features  (e.g., 
lines  or  recombination  edges)  distinctly,  and  calculates  probabilities  that 
photons  emitted  in  one  plasma  zone  will  be  absorbed  in  another.  By  taking 
into  account  various  line  broadening  mechanisms  and  geometric  factors  in 
summing  and  differencing  the  probabilities,  the  emitted  radiation  intensity 
as  well  as  the  photon  pumping  rates  (the  effect  of  re-absorbed  photons  on 
the  rate  equations)  can  be  calculated.  The  frequency-by-frequency  method, 
on  the  other  hand,  involves  tracing  the  photon  paths  via  rays  through  the 
plasma  at  individual  frequencies,  obtaining  a  solution  for  the  radiative 
intensity,  1^,  by  solving  the  equation  of  radiative  transfer.  The 
difference  between  the  methods  is  that  the  probabilistic  method  is 
essentially  the  ray-trace  technique  averaged  over  path  angle  and  over  the 
frequencies  of  the  spectral  feature.  Both  methods  have  been  carefully 
benchmarked  with  each  other  to  insure  that  the  predicted  spectral 

Intensities  and  energy  transfer  compare  well  over  the  regime  of 
applicability  of  the  models.  What  is  lost  by  the  probabilistic  method  are 
the  finer  details  of  the  line  profile  or  recombination  edge;  these,  on  the' 
other  hand,  are  accurately  calculated  by  the  ray-trace  technique.  Of 
course,  many  more  spectral  features  can  be  treated  by  the  former  (due  to 
decreased  memory  requirements)  as  well  as  being  significantly  more 

economical  In  terras  of  required  computational  time.  Presently,  we  can  use 
these  techniques  interchangeably  or  in  combination  with  each  other.  For 

IQ  'Tt 

example,  in  lower  density  gas-puff  or  wire  implosions  (N  =;  10  ions/cm  ) 
the  probabilistic  technique  is  completely  adequate  for  both  lines  and 

recombination  continua.  However,  at  very  high  densities,  like  laser- 
raicroballoon  target  simulation,  the  multi-frequency  ray-trace  method  is 

needed  due  to  the  dominance  of  Stark  and  opacity  broadening.  At  medium 
densities,  such  as  those  found  in  laser-foil  or  ion  beam-foil  interactions, 
a  "hybrid"  model  is  used,  whereby  lines  are  treated  probabilistically  and 
photoionization  is  calculated  via  the  multi-frequency  model. 
Bremsstrahlung  is  always  transported  frequency-by-frequency  since  no 
probabilistic  technique  has  as  yet  been  developed.  However,  due  to  the 
smooth  nature  of  the  free-free  source  function  with  photon  frequency,  very 
few  frequencies  are  required  for  accurate  transport. 
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The  final  statement  about  radiation  transport,  however,  is  simple  and 
to  the  point:  newer,  faster  transport  techniques  need  to  be  developed! 
The  present  algorithms  are  economical  enough  for  implementation  in  one¬ 
dimensional  hydro  models;  we  demonstrated  this  fact  in  a  paper  earlier  this 

O 

year.  However,  if  ionization/radiation  models  are  to  be  successfully 
integrated  with  2-D  codes,  new,  faster  techniques  are  needed! 

IV.  MODEL  IMPLEMENTATION  AND  APPLICATION 

A.  An  =  0  L-Shell  Transitions 

It  is  an  obvious  fact  to  those  who  work  at  low  densities  (primarily  the 
astrophysics  and  magnetic  fusion  workers)  that  the  bulk  of  the  non-K-shell 
line  emission  results  from  transitions  between  levels  of  similar  quantum 
number  (An=0  transitions),  where  n=2  in  the  L-shell,  n=3  in  the  M-shell, 
etc.  It  is  also  well-known  that  as  the  plasma  density  increases  to  values 
of  interest  to  the  audience  here,  the  An=0  radiation  is  overwhelmed  at  some 
point  by  the  An=l,2,...  emission,  due  to  the  fact  that  at  higher  densities 
the  excited  n=2  states  (in  the  L-shell)  become  collisionally  quenched;  at 
the  same  time,  the  rates  for  populating  the  higher  n  levels  become  great 
enough  that  the  lines  from  these  states  become  more  intense.  The  question 
most  relevant  to  model  builders  regarding  this  topic  is:  "Can  one,  at  some 
point  in  the  density,  neglect  the  excited  n=2  levels  in  the  model  and 
include  levels  charcterized  by  principle  quantum  number  along?" 

To  address  this  question,  the  ratio  of  the  An=0  emission  to  the  total 
line  radiation  for  the  carbonlike  ion  of  krypton  versus  ion  density  is 
shown  in  Fig.  3.  The  electron  temperature  was  fixed  at  3  keV, 
corresponding  to  roughly  the  maximum  of  the  carbonlike  abundance  (as  shown 
in  Fig.  4).  At  low  densities,  the  An=0  transitions  account  for  80%  of  the 
total  ionic  line  emission,  but  this  drops  to  less  than  a  percent  in  the 
high  density  limit.  The  interesting  fact  is  that  the  An=0  emission  is 
strong  up  to  densities  of  several  times  10  ,  indicating  that  including  the 

n»2  excited  states  is  imperative  for  modeling  krypton  Z-pinch  implosions  or 
plasmas  heated  by  x-ray  lasers,  non-coherent  radiation  sources,  or  ion 
beams  (where  low-temperature,  medium-density  plasmas  are  generated);  but, 
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perhaps  they  can  be  neglected  In  laser-plasma  interactions  where  the  low- 
tenperature  L-shell  emitting  region  is  usually  at  high  densities  during 
much  of  the  emission  phase.  Also,  the  density  at  which  the  An=0  radiation 
becomes  negligible  decreases  with  smaller-z  materials,  although  no  attempt 
has  been  made  to  quantify  this  with  our  other  models,  as  yet. 
Nevertheless,  all  indications  are  that  some  accounting  of  An=0  line 
emission  is  necessary  for  accurate  radiation  energetics  in  a  wide  range  of 
plasma  modeling  applications.  Because  of  the  limitations  of  several  other 
techniques,  such  as  the  "average-ion"  method,  no  account  of  An  radiation 
can  be  taken,  making  them  unsuitable  for  modeling  of  PRS  sources. 

B.  Non-Maxwellian  Electron  Distributions 

The  effects  of  non-thermal  electron  distributions  on  various  physical 

processes  in  dense  plasmas  is  a  topic  of  very  recent  interest.  Areas 

presently  under  investigation  include  such  fundamental  processes  as  heat 

Conduction,  laser  absorption  and  bremsstrahlung  emission,  for  example. 

Recently,  we  have  taken  an  interest  in  how  non-Maxwe 1 1 i an  distributions 

affect  the  ’  ionic  abundances  and  subsequent  atomic  radiation.  To  model 

these  effects  requires  complete  re-shaping  of  the  atomic  dat3  base,  since 

rate  coefficients  thermally  averaged  over  Maxwellian  distributions  are  not 

applicable  to  the  problem.  Instead,  tables  of  cross-sections  or  collision 

strengths  are  required,  enabling  one  to  include  the  electron  distribution 

as  an  input  variable,  over  which  the  data  can  be  averaged  to  obtain  rate 

coefficients.  This  problem  was  recently  addressed  in  an  N'RL  paper  which 

discussed  the  possibility  of  obtaining  population  inversion  and  gain  in 
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neon-like  systems  by  pumping  with  a  suprathermal  electron  <iis t ri but  ion . “  ‘ 
A  graph  depicting  the  ion  abundance  of  neonlike  iron  as  a  function  of  the 
hot  electron-to-cold  electron  number  density  is  shown  in  Fig.  4  for  fixed. 
Te(cold)  =  70  eV,  Tg(hot)  =  800  eV,  and  Ng(cold)  =  4x10-®  em~\  For  a 
purely  thermal  plasma  at  70  eV,  the  fraction  of  Fe  XVII  ions  is  below 
20%.  As  the  number  of  suprathermals  is  increased,  the  abundance  increases 
to  over  70%  and  then  to  zero  as  all  the  neonlike  electrons  are  burned 
off.  The  figure  to  note  here  Is  that  only  a  1%  concentration  of  hot 
electrons  can  totally  disturb  the  ionization  balance  calculated  due  to  a 
thermal  bath  alone.  Whereas  many  laser-plasma  interactions  at  high 
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irradiance  generate  this  percentage  of  suprathermals ,  it  is  not  clear  that 
neglecting  this  effect  in  modeling  such  plasmas  will  result  in  an  accurate 
assessment  of  the  population  distributions.  This  may  be  particularly  true 
of  K-shell  populations  where  x-ray  lines  from  K-shell  ions  are  used  to 
determine  the  local  electron  temperature  and  infer  the  extent  of  preheat 
due  to  electron  thermal  conduction,  flux  limits,  etc.  Some  of  the  recent 
unexplained  pre-heat  measurements'^  which  attribute  K-shell  resonance  line 
emission  deep  in  the  plasma  to  thermal  conduction  pre-heat  may,  in  fact,  be 
partly  due  to  enhanced  ionization  by  a  low-level  suprathermal  electron 
population,  as  shown  in  Fig.  4. 

C.  Time-Dependent  Effects 

Another  topic  of  current  interest  is  time-dependent  effects  on  the 
atomic  rate  equations:  the  deviation  from  ionization  equilibrium  due  to 
changes  in  plasma  parameters  occurring  faster  than  atomic  collisional  time 
scales.  The  standard  form  for  the  atomic  rate  equations  is 


Jr1  •  l  “ij  \  \  "ji 


where  is  a  collisional  or  radiative  rate  from  level  i  to  level  j.  For 

a  plasma  in  collisional-radiative  equilibrium,  the  time  derivative  is  set 
to  zero  and  the  set  of  steady-state  equations  is  solved  for  the  population 
densities.  If  the  plasma  is  changing  too  rapidly  to  relax  to  an  equilibrum 
state,  the  time-dependence  must  be  retained.  By  expanding  the  total  time 
derivative,  we  get  two  components: 

dN ,  dN.  3N. 

_1  .  _J.  +  v  — i 

dt  dt  dx 

The  first  term  is  the  explicit  time  dependence  and  comes  into  play  when 
local  plasma  variables  (Nj  and  Tg)  change  very  rapidly  compared  to  atomic 
transition  times.  Thus,  when  the  plasma  is  being  heated  rapidly, 
ionization  lags  behind  the  instantaneous  distribution  of  populations  as 
predicted  by  an  equilibrium  model,  and  a  lower  state  of  ionization  actually 
persists;  the  opposite  is  true  in  a  cooling,  or  recombining  plasma.  As  an 
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illustration  of  this  effect,  a  comparison  between  the  Z  cf  in  a  C07  laser- 
heated  plasma  as  calculated  by  a  CRE  and  time-dependent  CR  model  is  shown 
in  Fig.  5  as  a  function  of  time.  The  time-dependent  temperature  and 
density  profiles  for  one  zone  of  the  aluminum  plasma  were  taken  from  a 
LASNEX  hydrodynamic  simulation  provided  by  Los  Alamos.  The  standard  early- 
time  lag  in  ionization  followed  by  the  late-time  lag  in  recombination  is 
evident.  Since  the  atomic  collisional  time  scales  are  a  strong  function  of 
the  plasma  density,  that  curve  is  also  shown.  For  this  particular  shot, 
the  deviation  from  equilibrium  occurs  as  the  density  falls  through 
approximately  10  u  ions/cm  .  This  value  will,  of  course,  depend  on  the 
laser  pulse  width,  irradiance,  and  Z  of  the  target  material.  These 

effects  are  rather  well-known  and  since  quantitative  studies  have  been 

23  '’A 

published  by  Colombant  et  al.  and  Matzen  and  Pearlnan"  several  years 
ago,  we  will  not  dwell  on  it  further. 

Instead,  we  turn  attention  to  the  second  terra  in  the  equation,  the 
convective  derivative.  This  term  becomes  important  when  the  plasma  flow 
velocity  into  a  local  region  of  interest  is  faster  than  the  atomic 
collisions  governing  relaxation  to  equilibrium.  If  this  is  the  case,  the 
ionization  state  of  this  region  will  reflect  tempertures  and  densities  from 
some  other  location  in  the  plasma,  not  the  local  values.  This  term  is 

usually  neglected  in  most  calculations  since  either  the  flow  velocity  is 
assumed  to  be  small  or  the  ionization  state  is  calculated  on  a  Lagrangian 
hydro  mesh  which  moves  with  the  fluid,  obviating  the  need  for  invoking  the 
convective  derivative  in  the  rate  equations. 

Recently,  however,  we  have  become  aware  of  flow-velocity  effects  in 
conjunction  with  the  analysis  of  x-ray  spectra  from  implanted  tracer-dot 
experiments  at  NRL.  To  provide  a  little  background,  tracer-dot 

spectroscopy  is  a  new  diagnostic  technique  developed  at  NRL,  by  which  small 
"dots"  of  a  higher  Z  material  are  implanted  in  a  lower  Z  laser  target  in 
order  to  localize  in  space  the  emission  from  the  dot  material.  The 
spatially  resolved  K-line  spectra  thus  allow  the  plasma  profiles  to  be 

determined.  This  technique,  as  shown  schematically  in  Fig.  6  for  an 

aluminum  dot  in  a  carbon  target,  has  a  number  of  advantages  over  standard 

A  C  O 

laser-plasma  spectroscopy  and  has  been  discussed  in  detail  elsewhere. 

In  addition,  preliminary  studies  of  the  density  and  temperature  profiles  In 
the  blowoff  region  of  a  long  density  scalelength  laser  plasma  have  recently 


/ 
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been  published"  illustrating  the  valuable  diagnostic  contributions  that 
tracer-dot  spectroscopy  can  make.  However,  the  nature  of  the  technique 
(essentially,  spectral  measurements  made  in  the  lab  frame  of  a 
perpendicularly  flowing  laser-plasma  ablation)  lends  itself  to  being 
susceptible  to  the  flow  velocity  problems  alluded  to  earlier.  Although 
agreement  between  spectroscopy,  interferometry,  and  2-D  hydrocode  modeling 
was  good  in  the  long-scalelength  case,  errors  and  discrepancies  arose  when 
we  attempted  to  analyze  short  density  scalelength  plasmas.  We  surmised 
that  the  spectra  we  were  analyzing  were  in  fact  not  representative  of  the 
plasma  region  localized  by  the  spectrometer  slit,  but  were  colored  by  the 
rapidly  flowing  plasma  upstream  from  the  region  being  detected. 

To  test  this  hypothesis,  hydro  profiles  from  the  Eulerian  laser  code, 
26 

FAST2D,  were  transformed  to  represent  time-dependent  temperatures  and 
densities  which  reflected  the  flow  velocity  of  plasma  from  the  laser 
target.  These  profiles  were  then  post-processed  (assuming  a  steady-state 
flow)  using  a  time-dependent  version  of  our  CR  model  to  predict  K-line 
intensity  ratios  as  a  function  of  distance  form  the  target.  These  ratios, 
therefore,  represented  the  emission  that  would  be  measured  from  a  typical 
laser-plasma  interaction.  A  sample  set  of  line  ratios  versus  distance  from 
the  target,  Z,  is  shown  for  the  H-«/He-m  ratio  in  Fig.  7;  both  a  long- 
scalelength  and  short-scalelength  case  are  shown.  For  comparison,  the  set 
of  ratios  from  a  time-independent  post-process  (reflecting  the  actual 
temperature  and  density  in  the  plasma)  are  also  shown.  As  you  can  see,  the 
line  intensity  ratio  which  would  be  measured  (t-dep)  is  lower  than  those 
which  reflect  true  plasma  conditions  (t-ind)  at  Z  values  close  to  the 
target,  while  at  larger  Z's  where  the  temperature  (and,  hence,  the  line 
ratio)  drops,  plasma  flow  velocity  maintains  the  line  ratio  at  higher 
values.  As  was  expected,  the  faster  the  flow  velocity  (as  typified  by  the 
shorter  scalelength  case),  the  greater  is  the  discrepancy  between  the 
actual  and  flow-affected  spectra. 

How  are  these  effects  manifested  in  actual  spectral  diagnostics? 
Unfortunately,  the  affect  can  be  devastating.  In  Fig.  8,  a  plot  of  the 
plasma  electron  temperature  with  Z  is  shown;  both  the  hydrocode  profile  and 
3  pot  spectroscopy  diagnosis  are  shown  (for  the  long  density  scalelength 
case).  The  characteristic  signature  of  conventive  flow  effects  on  the 
spatially-resolved  spectra  is  evident:  underestimates  of  the  temperature 


near  the  target  (where  plasma  upstream  is  cooler)  and  overestimates  of  the 
temperature  far  from  the  target  (where  the  upstream  plasma  is  hotter).  For 
shorter  scalelength  plasmas  the  errors  are  expected  to  be  significantly 
larger  than  those  shown  in  Fig.  8.  This  is  a  very  disappointing  result 
from  the  standpoint  of  the  spectroscopist  since  it  implies  that  tracer-dot 
analysis  no  longer  depends  on  only  the  experimental  data  and  an 
ionization/radiation  model,  but  the  diagnoses  is  linked  to  a  hydrodynamic 
simulation  in  order  to  predict  the  effect  of  the  plasma  flow.  One  would 
like  to  benchmark  or  verify  the  validity  of  hydro  models  by  using  spot 
spectroscopy  techniques;  now  that  the  results  are  presumably  tied  with  the 
hydrocode  predictions  absolute  verification  is  a  bit  more  questionable  via 
the  technique. 

It  would  be  very  helpful  to  develop  a  method  by  which  the  flow  problem 

is  eliminated  since  it  will  manifest  itself  not  only  with  tracer  dots  but 

in  any  spectroscopic  analysis  where  the  emission  is  detected  perpendicular 

to  the  flow,  whether  it  be  low  concentration  tracer  material  doped  in 

targets  or  just  straightforward  spatially-resolved  full-target 

spectroscopy.  Are  there  any  ways  to  avoid  the  convection  effects?  One 

idea  we  are  presently  pursuing  at  NHL  is  to  investigate  using  only 

diagnostic  indicators  from  the  same  ionization  stage.  This  would,  in 

principle,  eliminate  non-equilibrium  ionization  effects  caused  by  either 

explicit  or  flow-velocity  time  dependence,  since  the  atomic  collisional 

rates  populating  levels  within  a  single  stage  of  ionization  are  usually 

much  faster  than  either  the  rates  between  ionization  stages  or  the 

hydrodynamic  time  scales.  If  this  proves  to  be  a  valid  assumption,  plasma 

temperature  could  be  determined  via  dielectronic  satellite  line  ratios  as 

opposed  to  the  present  method  of  comparing  hydrogenic  and  heliumlike 

2  7 

resonance  lines.  (A  recent  paper  describes  the  numerous  diagnostic 
indicators  possible  using  satellites  in  dense  plasmas.)  Although  the 
plasma  density  is  presently  determined  by  lines  from  the  same  stage  of 
ionization  (the  resonance-to-intercombinat  ion  line  intensity  ratio  in  the 
heliumlike  system),  errors  infiltrate  the  measurement  since  temperature  is 
required  in  order  to  determine  density.  Thus,  if  bogus  temperatures  are 
obtained  from,  for  example,  the  H-n/He-ai  line  ratio,  they  will  disturb 
the  density  mesurement.  If  an  alternative  temperature  diagnostic  is  used 
which  makes  accurate  predictions,  unaffected  by  flow  velocity,  then 


accurate  density  measurements  can  also  be  assured.  These  concepts  are 
presently  under  study  at  NRL  and  the  results  will  be  made  available  in  a 


2  Hydrogenlike-to-heliumlike  1  Ine  intensity  ratios  for  a  10“ 
aluminum  plasma  in  cylindrical  geometry  (diameter  =  115  pm), 
calculated  with  the  CRE  model.  Solid  curves  are  for  n  <  10  for 
both  ions;  dashed  curves  are  for  n  <  5  for  both  ions; 
alternating  dots  and  dashes  are  for  n  <  10  in  AA  XII  and  n  <  5 


log  Nj 


Fig.  3  Radiation  due  to  An  =  0  transitions  for  carbonlike  krypton  as 
compared  to  the  total  Kr  XXXI  line  emission  for  various  ion 
densities;  electron  temperature  is  constant  at  3  keV. 
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Ne-like  IRON 


.  4  The  fractional  ion  abundance  of  neonlike  ion  as  a  function  of  the 
ratio  of  hot  electron  density  to  cold  electron  density  (T  (cold)  =  70 
eV,  T  (hot)  =  800  eV,  N  (cold)  =  4x2020  cni3). 


TIME  (nsec) 

Fig.  5  The  of  a  zone  in  a  CO 2  laser-heated  aluminum  foil  simulation 

as  a  function  of  time,  comparing  equilibrium  CRE  and  time-dependent  CR 
model  predictions,  also  shown  is  the  plasma  ion  density. 


FI  Spectrometer 


6  Schematic  of  the  tracer-dot  spectroscopy  technique  for  obtaining 
laser-plasma  T  and  profiles. 


Fig.  8  The  electron  temperature  in  a  long  scalelength  laser-plasma  as 
predicted  by  tracer-dot  spectroscopy  of  a  115  pm-diameter 
aluminum  implant.  Also  shown  are  those  calculated  by  the  FAST2D 
hydrodynamic  code. 
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Dense  plasma  effects  on  AT-shell  dielectronic  satellite  lines 
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The  dielectronic  satellite  lines  from  the  n  =2  doubly  excited  levels  of  lithiumlike  and  heliumlike 
ions  have  been  studied  theoretically  in  dense-plasma  environments.  A  comprehensive  level  structure 
including  ground,  singly  excited,  and  doubly  excited  levels  was  included  in  a  collisional-radiative 
equilibrium  model  to  calculate  level  populations  and  the  satellite  intensities  self-consistently  for  car¬ 
bon,  aluminum,  and  argon  plasmas.  The  effects  of  inner-shell  excitation  rates  (calculated  via  the 
distorted-wave  approximation  with  configuration  mixing)  and  collistonal  ionization  and  recombina¬ 
tion  of  the  doubly  excited  levels  are  included  in  the  model  and  their  effects  on  the  satellite  intensities 
assessed.  The  well-known  heliumlike  intercombination-to-resonance  line  ratio,  used  in  density 
determinations,  is  shown  to  be  strongly  affected  by  blended  satellite  lines  at  high  density,  causing 
significant  error  in  the  measured  density  when  their  intensities  are  not  accounted  for. 


I.  INTRODUCTION 


Dielectronic  satellite  lines  result  from  the  spontaneous 
radiative  decay  of  doubly  excited  levels  of  atoms  and  ions, 
usually  appearing  in  the  radiation  spectrum  as  less-intense 
features  adjacent  to  resonance  lines  of  the  ion  with  the 
next-highest  charge  state.  In  particular,  a  subset  of  this 
general  family  of  lines,  due  to  transitions  of  the  type 
\s2l-2l2V  and  ls;2/-ls2/2/’,  have  received  much  attention 
due  to  their  value  as  plasma  diagnostic  indicators.  These 
lines  are  often  prominent  features  in  the  /f-shell  spectra  of 
high-temperature  plasmas  and  lie  slightly  to  the  long- 
wavelength  side  of  the  hydrogenlike  Lyman-a  and  heli¬ 
umlike  resonance  lines.  Much  effort  in  plasma  spectros¬ 
copy  has  been  devoted  to  their  detection,  classification, 
measurement,  and  diagnosis. 

The  first  laboratory  observation  of  these  lines  was  by 
Edlen  and  Tyren,1  but  more  recently,  a  number  of  labora¬ 
tory  observations  have  been  reported  on.2-  The  impor¬ 
tance  of  these  satellite  lines  in  astrophysical  plasmas  was 
first  pointed  out  by  Gabriel  and  Jordan/  who  presented  a 
detailed  classification  of  the  lithiumlike  satellites  in  that 
study,  while  the  line  intensities  and  population  mecha¬ 
nisms  were  described  in  more  detail  and  calculations  im¬ 
proved  in  later  papers/-"  Observations  of  both  helium- 
like  and  lithiumlike  satellites  in  solar  spectra  were  report¬ 
ed  by  a  number  of  workers,"- 14  while  further  laboratory 
studies  were  done  with  vacuum  sparks’’  and  laser- 
produced  plasmas. The  effects  of  plasma  microfields 
on  autoionization  rates  and  profiles  of  heliumiike  satellite 
lines  was  first  described  by  Davis  and  Jacobs,"  while 
high-density  collisional  effects  were  first  pointed  out  by 
Weisheit22  and  later  developed  into  a  density  diagnostic 
for  laser-produced  plasmas  by  workers  in  the  Soviet 
Union/1-57  Seely,  Lunr.ey,  and  co-workers^-’1  extended 
the  technique  to  compressed  laser-fusion  microballoon  tar- 
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gets,  while  Duston  and  Davis'2  showed  that  opacity  ef¬ 
fects  could  alter  density  determinations  in  menially  con¬ 
fined  plasmas.  Meanwhile,  lithiumlike  satellites  were  the 
subject  of  continued  studies:  The  effects  of  non- 
Maxwellian  electron  distributions  on  satellite  intensities 
were  investigated  by  Gabriel  and  Philips,"  density  effects 
were  seen  in  vacuum  spark  spectra  by  Kononov  and  co- 
workers,'4  Bely-Dubau  et  al.~ 54  extended  calculations  to 
Is:n/- \s2l'nl— type  configurations  with  n  >  3,  Bitter 
et  a/.,3’-41  the  TFR  group42  and  Kallne  el  a/.4'-45  ob¬ 
served  detailed  satellite  structure  in  tokamaks,  Datla 
el  al.M’  determined  electron  temperature  in  a  neon-filled  0 
pinch,  and  Jacobs  and  Blaha4'  calculated  the  effects  of 
collisions  between  doubly  excited  levels  on  the  satellite  in¬ 
tensities.  Armour  et  al.ir  and  Boiko  et  a/.44  have  report¬ 
ed  on  comprehensive  satellite  line  observation  and  classifi¬ 
cation  in  beam-foil  excited  spectra  and  laser  plasmas, 
respectively. 

The  purpose  of  this  study  is  to  continue  the  work  begun 
by  Jacobs  and  Blaha41  of  describing  satellite  line  behavior 
in  dense  plasmas  by  including  electron  collisional  effects 
explicitly.  However,  we  extend  the  complexity  of  their 
model  by  solving  for  the  doubiy-excited-state  population 
densities  and  satellite  line  intensities  self-consistently  with 
singly-excitcd-stale  densities  and  resonance  and  intercom¬ 
bination  line  intensities.  The  details  of  the  model, 
described  in  Sec.  II.  include  ionization  dynamics,  level 
structure,  and  rate  coefficient  calculations.  The  behavior 
of  the  heliumiike  and  lithiumlike  satellite  lines  as  a  func¬ 
tion  of  plasma  densitv  and  temperature  is  studied  in  Sec. 
III.  including  discussions  of  several  diagnostic  line  intensi¬ 
ty  ratios  In  order  to  allow  comparisons  of  the  satellite 
lines  as  a  function  of  atomic  num'  er  Z.  the  calculations 
are  performed  for  high-donsity  carbon,  aluminum,  and  ar¬ 
gon  plasmas.  Because  of  the  magnitude  of  this  subject,  we 
have  limited  the  scope  of  this  paper  to  include  only  opti¬ 
cally  thin  plasmas.  The  effects  of  photon  self-absorption 
on  the  satellite  line  intensities  will  be  investigated  in  a  fu¬ 
ture  paper. 
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II.  THEORETICAL  MODEL 

The  method  used  to  investigate  the  behavior  of  the 
dielectronic  satellite  lines  is  based  on  a  completely  self- 
consistent  collisional-radiative  equilibrium  (CREr0  calcu¬ 
lation.  Rate  equations  of  the  form 

d.\,  _  _ 

-jr-S^VrS^v,.  (i) 

a  i  i 

are  solved  for  each  atomic  energy  level  included  in  the 
model  for  an  equilibrium  plasma  [i.e.,  the  left-hand  side  of 
Eq.  (1)  is  set  equal  to  zero],  with  Nt  as  the  population  den¬ 
sity  for  level  i  and  representing  the  total  transition 
rate  from  level  t  to  level  j.  The  rate  coefficients,  Wtp  are, 
in  general,  functions  of  the  plasma  density,  temperature, 
and  local  photon  field,  and  are  discussed  in  detail  below. 
The  equations  are  solved  independent  of  a  fluid  descrip¬ 
tion  of  the  plasma  motion;  thus,  the  density,  temperature, 
and  plasma  size  are  the  independent  variables.  Once  the 
set  of  population  densities  has  been  calculated,  the  indivi¬ 
dual  line  and  continuum  intensities  can  be  determined  and 
a  theoretical  radiation  spectrum  calculated. 

A.  Energy  levels 

In  order  to  accurately  assess  the  effect  of  varying  Z  on 
the  satellite  intensities,  the  identical  level  structure  was 
used  for  all  three  materials  studied.  All  ground  states 
were  included  in  addition  to  the  n  =2—4  levels  of  the  hy¬ 
drogenlike  ion,  the  ls2s3S,  ls2p3P,  ls2s'S,  Is  2p  'P, 
n  =  3  triplet,  n=  3  singlet,  n=  4,  n  =5,  2 s2tS,  2s2p3P, 
2 p2  3P,  2p2  lD,  2s  2 p  'P,  and  2 p1  'S  levels  of  the  heliumlike 
ion  and  the  \s22p,  ls23s,  ls23p,  lj23d,  n  =4,  n  =  5, 
ls2s22S,  U2s('P)2p2P,  \s2s(3P)2p  2P,  \s2p22D,  \s2p22P, 
and  ls2p2  2S  levels  of  the  lithiumlike  ion.  At  the  densities 
of  interest  in  this  study,  ion-ion  collisions51,52  should  force 
individual  j  components  of  a  single  nl  level  into  statistical 
equilibrium  with  one  another,  hence,  only  nl  levels  (or,  in 
cases  of  high-rt,  only  n  levels)  are  treated  in  the  model  for 
the  purposes  of  collisional  couplings.  In  order  to  treat  in¬ 
dividual  satellite  line  components  from  j  levels,  however, 
the  radiative  couplings  between  doubly  and  singly  excited 
levels  are  calculated  as  individual  j  components.  While 
the  individual  radiative  decay  rates  are  still  averaged  over 
j  to  obtain  the  nWevel  transition  rate,  this  scheme  allows 
us  to  accurately  represent  the  individual  satellite  line  in¬ 
tensities.  The  energies  of  the  levels  included  in  our  atomic 
model  were  calculated  using  the  relativistic  Hartree-Fock 
atomic  structure  code  provided  by  Cowan.55 

B.  ColliiiOBal  couplings 

The  rate  coefficients  employed  in  this  calculation  were 
obtained  by  a  number  of  methods.  All  levels  of  ion  Z 
were  coupled  to  the  ground  state  of  ion  Z  +  I  by  collision¬ 
al  ionization  (exchange  classical  impact  parameter 
method5*)  photoionization  (hydrogenic  approximation55 
with  Karzas-Latter  Gaunt  factors5,j,  and  collisional  and 
radiative  recombination  (the  detailed  balance  of  the  above 
rate  coefficients!.  Adjacent  ions  were  coupled  by  dielec¬ 
tronic  recombination,  with  rate  coefficients  provided  by 


Jacobs  and  Davis.5  All  dipole-allowed  and  some  impor¬ 
tant  forbidden  radiative  transitions  were  coupled  by  spon¬ 
taneous  radiative  decay.  The  oscillator  strengths  and  Ein¬ 
stein  A  coefficients  were  also  calculated  using  Cowan's 
structure  code,  as  were  the  autoionization  raies  for  the 
doubly  excited  levels.  The  autoiomzatton  and  radiation 
rates  calculated  for  this  study  have  been  compared  with 
those  calculated  by  Vainshtein  and  Safranova.5"  In  fact, 
due  to  the  widespread  use  this  comprehensive  set  of 
dielectronic  satellite  data  has  had,  we  not  only  compared 
the  actual  data,  but,  in  addition,  used  their  data  for  pre¬ 
liminary  satellite  line  intensity  studies  while  the  new  data 
set  was  being  generated  for  the  present  calculations.  Com¬ 
parison  of  the  two  sets  of  values  for  wavelengths  of  the  sa¬ 
tellites  and  the  spontaneous  decay  rates  showed  very  little 
difference,  typically  less  than  0.1%  and  10%,  respectively, 
for  all  three  of  the  materials  investigated.  The  values  gen¬ 
erated  by  Cowan’s  structure  code  for  autoionization,  how¬ 
ever,  differed  by  typically  50—100%  for  most  transitions, 
but  varied  by  as  much  as  a  factor  of  3  in  some  cases.  The 
largest  differences  occurred  in  carbon,  the  smallest  in  ar¬ 
gon.  However,  more  importantly,  we  were  interested  in 
what  effects  these  differences  would  have  on  the  diagnos¬ 
tics.  It  was  found  that,  in  most  cases,  little  or  no  change 
occurred  in  the  satellite  line  intensities  when  one  set  of 
data  was  replaced  by  the  other.  The  largest  differences 
were  on  the  order  of  10—20%,  but  these  discrepancies 
were  manifested  in  only  3  of  the  12  satellite  line  blends 
studied;  all  others  experienced  less  than  10%  changes 
Thus,  differences  in  a  plasma  temperature  or  density 
determination  could  be  expected  to  be  of  the  same  order 
using  the  two  sets  of  data.  Since  the  bound-state  and  con- 
tinuum  wave  functions  calculated  using  Cowan's  model 
represent  a  distinct  improvement  over  the  Coulomb- Born 
wave  functions  used  by  Vainshtein  and  Safranova  [the  ap¬ 
proximate  nature  of  their  autoioinizalion  rates  is.  in  fact, 
acknowledged  and  stated  by  the  authors],  we  have  chosen 
to  use  the  more  recently  calculated  data  in  our  study  The 
reader  is  cautioned,  however,  that  the  differences  quoted 
here  result  from  comparisons  performed  in  a  dense- 
plasma  environment;  the  discrepancies  may  well  be  ampli¬ 
fied  at  lower  densities  where  the  plasma  approaches  a 
more  coronal  state  and  autoionization  plays  a  more  instru¬ 
mental  role  in  determining  doubly  excited  level  popula¬ 
tions. 

The  ground  states  of  the  hydrogenlike  and  heliumlike 
ions  were  coupled  to  their  appropriate  doubly  excited  lev¬ 
els  by  electron  capture;  the  rate  coefficients  for  this  pro¬ 
cess  were  obtained  by  detailed  balance  of  the  autoiomza- 
tion  rates  (the  opposite  process).  The  doubly  excited  levels 
were  coupled  to  the  n  =  2  levels  of  the  next-highest  ion  by 
collisional  ionization  and  recombination  using  the  same 
methods  slated  earlier.  Although  these  processes  are  noi 
thought  to  be  important  until  very  high  densities,  we 
sought  to  quantify  their  effect  and,  therefore,  included 
them  in  the  model  Singly  excited  levels  were  coupled  to 
ground  states  and  also  to  each  other  via  electron  impact 
excitation  [the  Coulomb-Bom  distorted-wave  'DWi'“  code 
of  Blaha  was  used].  However,  for  coupling  the  Is",  1s2j. 
and  Is  2p  levels  of  the  heliumlike  ion.  the  more  accurate 
rates  by  Pradhan  ei  al which  include  the  effects  of  au- 
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toionizing  resonances,  were  used,  Finally,  doubly  excited 
levels  of  a  given  ion  were  coupled  to  each  other  and  to 
singly  excited  levels  by  electron  impact  excitation  and 
deexcitation.  The  calculation  of  these  coefficients  is  treat¬ 
ed  in  detail  in  Sec.  II C. 

C.  Inner-shell  and  cross-coupling  transitions 

In  addition  to  electron  capture  and  collisional  recom¬ 
bination,  the  doubly  excited  levels  can  be  populated  by 
inner-shell  electron  impact  excitation  of  an  n  =  1  electron 
of  the  lithiumlike  and  heliumlike  ions.  For  some  states, 
this  process  is  the  dominant  populating  mechanism.  The 
rate  coefficients  for  these  transitions  were  calculated  in 
the  following  way.  After  the  energy  levels  had  been  deter¬ 
mined,  electron  impact  cross  sections  were  generated  be¬ 
tween  the  various  LS  terms  using  the  DW  code.  In  these 
calculations,  the  relativistic  Hartree-Fock  radial  wave 
functions  generated  by  the  atomic  structure  code  were 
used  to  describe  the  bound  electrons  while  the  continuum 
electronic  states  were  represented  by  distorted  wave  func¬ 
tions  which  take  into  account  the  non-Coulombic  electro¬ 
static  interaction.  Since  the  plasma  electrons  will  have  en¬ 
ergies  of  the  order  of  the  threshold  excitation  energies  of 
these  inner-shell  transitions  for  the  parameters  of  interest 
in  this  study,  electron  exchange  effects  were  included  in 
the  cross-section  calculations. 

Since  the  threshold  transition  energies  between  doubly 
excited  levels  are  typically  2  orders  of  magnitude  smaller 
than  electron  energies  at  the  plasma  temperatures  of  in¬ 
terest  here,  exchange  effects  were  not  included  in  the  cal¬ 
culation  of  the  cross-coupling  cross  sections  between  dou¬ 
bly  excited  levels.  However,  it  was  necessary  to  calculate 
the  cross  sections  at  these  very  high  electron  energies.  In 
general,  it  is  not  practical  to  do  the  summation  over  all 
the  angular  momenta  required  to  evaluate  the  collision 
strengths  n  at  these  high  energies.  Therefore,  the  col¬ 
lision  strength  was  usually  calculated  fully  up  to  some 
large  value  of  electron  energy  E  (typically  20  times  thresh¬ 
old)  and  extrapolated  to  higher  values  of  E.  For  dipole 
transitions  the  Bethe  approximation01  was  used,  while  it 
was  assumed  that  12  varied  as  E~~  for  transitions  involv¬ 
ing  a  change  of  spin.  For  quadrupole  transitions,  fl 
should  approach  a  constant  value  as  E  becomes  very  large, 
however,  this  limit  is  approached  computationally  with 
considerable  effort  and  expense  due  to  the  large  number  of 
angular  momentum  terms  required  for  convergence. 
Hence,  nondipoie,  non-spin-change  transitions  were  plot¬ 
ted  to  examine  high- /T  behavior  and  were  extrapolated  ac¬ 
cordingly. 

Configuration  and  term  mixing  of  the  various  energy 
levels  and  transitions  were  included  in  the  collision 
strength  calculations  using  the  methods  described  by 
Jacobs  and  Blaha/1  The  mixing  coefficients  were  taken 
from  Cowan's  structure  code.  Rate  coefficients  were  then 
calculated  by  thermally  averaging  the  collision  streng'hs 
over  a  Maxwellian  di  unbutton. 

Shown  in  Fig  I  are  the  inner-- neil  collisional  excitation 
rate  coefficients  for  lithiumlike  Al  XI  versa:,  electron  tem¬ 
perature.  while  those  for  heliumlike  Al  XII  are  shown  in 
Fig  2  As  was  expected,  the  Al  XI  results  reproduced  very 


FIG.  !.  Inner-shell  electron  impact  excuation  rale  coeffi¬ 
cients  for  singly  to  doubly  excited  levels  of  the  lithiumlike  Al  XI 
ion  vs  electron  temperature. 


well  the  calculations  of  Jacobs  and  Blaha.  The  difference 
between  the  present  calculation  and  theirs  is  that  they 
used  the  bound-state  wave  functions  of  Clementi  and 
Roetti.0"  However,  these  wave  functions  were  found  to  be 
virtually  identical  to  the  Hartree-Fock  wave  functions 
generated  by  Cowan's  structure  code.  The  lithiumlike  ion 
cross-coupling  excitation  rate  coefficients  for  Al  \i  are 
shown  in  Fig.  3;  the  key  to  the  transitions  is  in  Table  1. 
As  was  stated  earlier,  ihe  plasma  temperatures  of  interest 
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FIG.  3.  Cross-coupling  electron  impact  excitation  rate  coeffi¬ 
cients  between  doubly  excited  levels  of  the  lithiumlike  AlXI  ion 
vs  electron  temperature. 


here  < typically  above  200  eV  for  aluminum,  where  Alxi! 
is  expected  to  be  abundant)  are  well  above  threshold  ener¬ 
gies  for  these  transitions;  this  can  be  seen  clearly  in  Fig.  3. 

The  collision  strengths  for  the  heliumlike  ion  cross¬ 
coupling  transitions  versus  electron  energy  (in  units  of 
E/AE.  where  AE  is  the  threshold  energy)  are  shown  in 
Fig.  4  for  carbon  and  Fig.  5  for  argon;  the  key  to  the  tran¬ 
sitions  is  in  Table  II.  Also  shown  as  dashed  curves  in 
Figs.  4  and  5  are  the  calculations  of  Goett,  Sampson,  and 
Clark,03  and  good  agreement  between  their  results  and  the 
present  calculation  is  evident  for  most  transitions  shown. 
Their  results  for  ft  were  generated  using  hydrogenic  atom¬ 
ic  wave  functions  for  the  bound  electrons  and  Coulomb 
functions  for  the  free  electrons,  with  Z:  scaling  and  in  in¬ 
termediate  coupling.04  In  particular,  the  2s2('S)  and 
2p:(lS)  levels  are  strongly  coupled  in  the  Arxvn  system, 
and  the  mixing  must  be  included  in  the  calculation  of  ft. 
The  effects  of  mixing  as  a  function  of  Z  can  clearly  be 
seen  by  comparing,  for  example,  the  2pz(  'P)-2s  2p\  'P) 
collision  strengths  in  Figs.  4  and  5;  the  behavior  of  ft  vs 


TABLE  I.  Transitions  among  the  lithiumlike  doubly  excited 
states  corresponding  to  the  electron  impact  excitation  rate  coef¬ 
ficients  shown  in  Fig.  3. 

Key  Transition 

1  ls2s2:S-l$2s('Pi2p  ;P 

2  \s2s'--S-\s2s<‘P)2p:P 

3  \s2s::S-is2p::D 

4  Is Zt< 'P)2p  :P  -ls2p2  :D 

5  Is2s('Pi2/>  P- 1  s  2p::P 

6  Is  2s(  'Pi2/>  "P-  Is  2p'  ;S 

’  Is  2si  'P  )2p  :P-]s2p :  D 

8  Is2ri  'P)2p  P-  Is2p:  2P 

4  !s2si'Pi2p;P-ls2p:;5 

10  Is 2p  :D  -  Is  2p‘  ;S 


FIG  4.  Cross-coupling  electron  impact  collision  strengths  be¬ 
tween  doubly  excited  levels  of  the  heliumlike  C  v  ion  vs  incident 
electron  energy  lin  threshold  energy  units’ 


FIG  5  Cross-coupling  electron  impact  collision  strengths  be¬ 
tween  doubly  excited  levels  of  the  heliumlike  \r  XX  II  ion  vs  in 
cident  electron  energy  on  threshold  energy  units 
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TABLE  II  Transitions  among  the  heliumlike  doubly  excited 
levels  corresponding  to  the  electron  impact  collision  strengths 
shown  in  Figs.  4  and  5 


Key 

Transition 

1 

2s:'S-ls2p  ' P 

2 

h'-'S-hlp'P 

3 

hip  'P-lp-'P 

4 

2s Ip  'P-lp'-'D 

5 

hip  'P-hlp  'P 

6 

hip  'P-lp'-'S 

7 

lp  :  P—lp1  ‘D 

8 

Ip- 'P-hlp  'P 

9 

lp:  'P-lp:  '5 

10 

lp'-'D-hlp  P 

11 

Ip-D-lpt'S 

12 

hip  'P-lp- 'S 

E/\E  is  totally  different  for  argon  and  carbon.  No  com¬ 
parison  of  transitions  between  similar  electron  configura¬ 
tions  Q.IV-V21')  is  shown  since  the  Coulomb-Bom 
method  of  Goett  et  at.  did  not  treat  this  type  of  transi¬ 
tion. 

D.  Dielectronic  satellite  lines 

As  discussed  earlier,  the  doubly  excited  levels  for  lithi¬ 
umlike  and  heliumlike  ions  give  rise  to  a  large  number  of 
satellite  lines  via  dipole-allowed  radiative,  transitions  to 
singly  excited  levels.  These  lines  are  represented  schemati¬ 
cally  in  Fig.  6  (lithiumlike  lines)  and  Fig.  7  (heliumlike 
lines).  For  classifying  the  lithiumlike  satellites,  we  use  the 
well-known  letter  designation,  first  introduced  by  Ga¬ 
briel.10  For  the  heliumlike  satellites,  we  use  a  simple  num¬ 
bering  scheme  to  define  the  lines.  A  complete  listing  of 
the  satellite  line  transitions,  designations,  theoretical 
wavelengths,  and  radiative  decay  and  autoionization  rates 
used  in  this  study  is  given  in  Table  III  for  carbon,  alumi- 


FIG.  6.  Schematic  of  the  excited  level  system  of  the  lithium- 
like  ion  which  gives  rise  to  the  n  -l  dielectronic  satellite  lines 
adjacent  to  the  heliumlike  resonance  line. 


FIG.  7.  Schematic  of  the  excited  level  system  of  the  helium¬ 
like  ion  which  gives  nse  to  the  it  =2  dielectronic  satellite  lines 
adjacent  to  the  hydrogenlike  resonance  line. 


num,  and  argon.  In  order  to  allow  visualization  of  typical 
satellite  line  spectra,  theoretical  carbon,  aluminum,  and 
argon  K-shell  spectra  in  the  region  of  the  hydrogenlike 
and  heliumlike  resonance  lines  are  shown  in  Fig,  S.  The 
resonance  lines  have  been  truncated  to  show  the  detailed 
structure  of  the  satellite  lines.  Worth  noting  is  the  relative 
shift  in  several  of  the  lines  with  increased  Z,  such  that 
lines  from  different  electron  configurations  overlap  or  be¬ 
come  blended  as  a  result  of  various  line  or  source  broaden¬ 
ing  mechanisms.  Note  also,  the  proximity  of  several  of 
the  satellites  (m,rt,s,l)  to  the  Is2- Is 2p  }P  intercombination 
(IC)  line,  such  that  they  can  appear  as  almost  a  single 
blended  line  (particularly  at  high  plasma  densities,  where 
lines  are  collisionally  broadened).  It  will  be  shown  later 
that  this  can  have  significant  impact  on  density  diagnos¬ 
tics  in  laser-produced  and  inertially  confined  plasmas. 

III.  RESULTS 

The  main  emphasis  in  this  section  will  be  on  the 
behavior  of  the  n  =2  dielectronic  satellites  as  a  function 
of  plasma  density  and  temperature.  As  stated  earlier,  all 
investigations  presented  here  assume  an  optically  thin 
plasma,  calculated  by  ignoring  the  effects  of  opacity  in 
our  CRE  model  imore  detailed  discussion  of  the  theoreti¬ 
cal  ionization-radiation  model  can  be  found  else¬ 
where65’’6).  A  detailed  study  of  opacity  effects  on  satellite 
lines  will  be  presented  in  another  paper  to  follow. 


A.  Ion  population  distributions 

Prior  to  discussion  of  the  satellite  line  intensities,  it  may 
prove  useful  to  present  the  ion  population  abundances  in 
the  parameter  range  of  interest  to  this  work  This  will  al 
low  later  reference  to  the  relative  populations  of  the  ions 
when  explaining  the  effects  of  individual  collisional  pro- 
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TABLE  HI.  List  of  transitions,  satellite  line  key,  wavelengths,  radiative  decay  rates,  and  autoionization  rates  for  the  satellite  lines 
included  in  this  study.  'Numbers  in  parentheses  denote  power  of  10 


Transition 

Key 

K 

iA> 

Carbon 

A 

•  sec  ~ 1 1 

r 

see  *1 

A 

(A' 

Aluminum 

A 

s  cc~l) 

r 

isee" 

•A» 

Argon 

A 

see  " 1 1 

r 

'sec "  ’  • 

1j;2j  1s2s(  'P  ‘2p  'P;  ; 

r 

41  394 

8  46.11! 

6.00*  12> 

7  8496 

2,56. 13) 

7  "?"1  ■  ]  ^  > 

3  9839 

8  88(  1 3 

1.28,13' 

ls:ls  :S,n-\s2st'P'Zp-Py~ 

4 

41.342 

8  49.11) 

5  8i  );• 

7  8475 

2,61(13. 

4.37' 1*. 

3.4816 

I.05d4, 

2.40.12) 

Is  :2s  :Su:-\s2s<'P\lp:P:  : 

r 

40  860 

9.04.10) 

4  807  3. 

7.8097 

3  66.12' 

7.3  M3) 

3  9692 

2.30U3) 

7,23(13. 

lsJ2s  !SW2~ts2s(>P)2p-P\,! 

5 

40.859 

8.78H0) 

4,81.13. 

7.8085 

2.281 12; 

7.64,13) 

3  9681 

6  29<  12' 

S.2'131 

ls:2p  2P|/j—  Is  2si:Slr. 

p 

43  684 

2 .09<10) 

7.1473/ 

8  0644 

9,54. 11) 

L04. 14) 

4.0660 

1  98(12) 

1 .08' 14) 

ls!2p  lPm.  —  ls2a::S|,j 

0 

43.686 

4,16.10) 

774(13. 

8.0732 

4.77(111 

1.04(14) 

4.0703 

3.13.12) 

1.08(14) 

\s:2p:P,/1-ls2p22Dir. 

k 

41.483 

3,62.1  li 

7  12«13» 

7.8683 

1.37(13) 

115*141 

3.9892 

6  33(13. 

1.16(14' 

ls:2p  2P,/:-lslp2  lDxr. 

1 

41.484 

5  46.10) 

7  1273) 

7  3714 

3.96.11) 

1.15,14. 

3  9932 

5  29. 13' 

1.16(14) 

ls:2p  :P\r.~  ls-p:2D<r. 

J 

41,485 

4.15.1H 

772.13. 

7  j*7|5 

1.34U3. 

1.17(14 

3  9923 

3.5F IP 

1.25.14' 

\s:2p  !Pi/;  — Is  2p:1Pir_ 

a 

41.425 

1.16'12l 

3  3070) 

’.85'} 

3.71.13) 

2.10(12' 

3.4854 

1  48  14, 

1 .00i  1 5 

ls:2 p  '~P\n  —  \s2p~~P \/i 

b 

41  424 

2.13(11) 

3.30701 

7.8538 

4.61(13) 

2  10'12) 

3A814 

9  5402 

1  .Ob'  1 5 

ls:2p  :Pi/j— Is  2p22P.r. 

c 

41.428 

4.56. 11) 

7,8609 

1.32(13' 

1,80.10) 

3.9897 

4  60(13. 

1  56. 1 1 

U-2p:P,/2-ls2p:1Pl.: 

d 

41.426 

9.17,111 

7.85-4 

2.90.13. 

1  SO  1 0 

3.4856 

1  '.9(14' 

1.56.11' 

ls:2p  :Pl/z-U2p::S„’ 

n 

40.736 

1.41(11) 

9.77(12) 

7.79g7 

3.47(12 

1  81.13 

3  963 ! 

8.65:12) 

2.04,13' 

\s2lp1Pi/2-\s2p:2Sl/: 

m 

40.738 

2.88(11) 

9.7772) 

-.8022 

9,65.12. 

l.SL  13 

3  46”, 

4.45.13 

2.04-13 

lslr  ’Si—  2i 2p  ’P0 

i 

34.525 

7.41(11) 

9  12'  12) 

7.248* 

1.73.12 

1  I0-13) 

3  ’652 

6.4  M3’ 

1  16.13 

lsls  ,S>-2s2p,P[ 

2 

34.523 

7.41.11) 

9.1272) 

7.2521 

1.7;.  12' 

1  20<  1 3 

3  *63’ 

6.58  15 

1  24  13' 

U2s  'S,-2s2p  ’/>; 

3 

34.520 

74L1H 

o  1 2<  12' 

7.2537 

l.'2’!2 

1.14  131 

3  *601 

b.44. 13' 

1.16.13 

Is2p  ’Pt-2p:  'P0 

4 

34.594 

1.42(12) 

1.00(9. 

7.2605 

3.34  15- 

1.57.11' 

3  “668 

1.25  14 

9  09;]  ] 

Ulp’Po-lp'-'P, 

5 

34.591 

4.77(11) 

7  2575 

1. 14.13. 

3  *656 

4  25'  13( 

{sip  lP\—2p'  ’F, 

6 

34.592 

3.58(111 

7.2587 

8  44, :  2 

3.7049 

3  14.13 

ls2p  'Pi—lpl,P\ 

7 

34.594 

5.96(11) 

7.2M8 

1.42'  13' 

3.7682 

5  29'  13' 

\s2p  'P,  -2p:  ,P! 

8 

34.589 

3.58(111 

1.1070) 

7.2558 

8.56. 12) 

2/>8'  121 

3  *62 1 

3  35' 13' 

1.64(13' 

Is  2p  iPi-lp22P2 

9 

34.591 

1.07(12) 

1.10(10) 

7,2588 

2  531 3 » 

2.08,12' 

5  "655 

8.98-13' 

1  64.13- 

\sls'So--1slp'Px 

10 

34.283 

7. 50(11) 

9.61(13) 

7.2284 

1  72(13. 

1.51(14. 

3.7536 

6.39'  13 

1  53,14 

Islp  'P\ -Is2  ’So 

11 

35.442 

3.34(11) 

1.65(14) 

7.3635 

7,66.12. 

2.59*  14> 

3  8058 

2.41. 13' 

2.59,14; 

ts2p  'Pl-2p:,D: 

12 

34.649 

1.32(12) 

1,62(14. 

7.2710 

3.24.13' 

2,66.14. 

3  -b94 

1.18'  14. 

2.59:14; 

U2p'Pl-2p1'S0 

13 

33.973 

1.21(12) 

7,20(12. 

7.1049 

2.77(13) 

1.26.13' 

3.7422 

1.06.14 

1.71(13' 

cesses  on  the  various  satellite  it  it  isities.  Shown  in  Fig.  9 
are  the  total  (ground  plus  singly  excited  states!  ton  popula¬ 
tion  distributions  versus  electron  temperature  for  carbon, 
aluminum,  and  argon.  Only  the  populations  for  the  lithi¬ 
umlike,  heliumlike,  hydrogenlil.e,  and  fully-stripped  ions 
are  shown,  since  these  are  the  only  ions  strongly  affecting 
the  satellite  lines.  The  ion  density  was  fixed  at  10:o 
ions/cm'  for  all  these  materials;  at  this  density,  most  of 
the  ionic  population  still  reside;  in  the  ground  state,  with 
the  exception,  of  course,  of  the  ls'2p  level  of  the  lithium- 
like  ion. 

B.  Satellite  intensities  versus  N,  and  T, 

In  order  to  more  clearly  (and  conveniently i  present  the 
satellite  line  dependence  on  electron  temperature  i  T,'  and 
total  ion  density  f.Y/i.  the  satellite  line  intensities  will  be 
presented  normalized  to  their  associated  resonance  line  in¬ 
tensity.  i ,e„  all  heliumlike  satellite  intensities  arc  divided 
by  the  Iv-ls2p7’  ’RES'  intensity  and  all  hydrogenhke 
satellites  are  divided  by  the  ls-2p  ' L.yman-rt>  intensity  In 
Fig  I'-,  the  heliumlike  satellite  intensities  are  presented  a- 
a  ;  m.tion  ,>(  tor  iliimmuit'  plasma  at  1 1  ■ 

ini-  -.  rile  uuliv  ideal  .oinponems  have  Seen 
gp'i. with  f'th-T  lines  tr..m  il.e  same  nl  'L>  leu-1  since 


this  is  how  they  normally  appear  in  a  high-density -plasma 
spectrum.  The  individual  j  components,  however,  can 
often  be  resolved,  al  lower  densities,  e.g.,  in  tokamak  or 
astrophysical  plasmas.  As  ihe  nuclear  Z  is  increased  and 
the  LS  representation  becomes  less  applicable,  j  com¬ 
ponents  begin  to  shift  in  wavelenglh,  and  the  groupings 
shown  in  Fig.  10  are  not  convenient  in  terms  of  measur¬ 
able  spectral  quantities  In  argon,  for  example,  sve  see 
from  Fig.  8  that  the  groupings  (/./>,  tr.a.d i,  it.t.m),  ( q.b\ 
(L,c),  lo.pK  and  i/i'  are  more  appropriate  for  experimental 
purposes.  The  consistent  trend  displayed  by  all  the  satel¬ 
lites  in  Fig.  10  of  falling  off  with  increasing  T,  is  due  to 
different  factors  which  depend  on  the  dominant  collisional 
process  populating  the  doubly  excited  level.  The  lines 
•a.b.c.dK  iq.ri,  and  emanate  from  doubly  excited 

levels  which  are  populated  prunarilx  by  inner-shell  excita¬ 
tion  from  lithiurnlike  levels  As  the  ratio  of  A 1  X l  to 
Al  XII  population  drops  with  7,,  so  do  these  satellite  in¬ 
tensities.  The  i  i.K.i  .  ->.p  ,  and  ‘  >.f  satellites,  on  the  oth¬ 
er  hand,  come  from  levels  populated  he  electron  capture 
from  the  i  .  lev  1  T  :a  capture  rale  coefficient  peaks  .it 
lower  temperature  does  ihe  Iv-lOp  excitation  rate 
coefficien!  due  the  .hllcrence  n  iraiisiinn  energies, 
am!  I li Us,  these  sirc'iOis  .ii,o  t  ill  at:  w t r h  increasing  /, 
From  Fig  111  u  s  clearly  ipparcni  why  ihe*c  line  ratios 
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FIG.  8.  Theoretical  spectra  showing  the  hvdri'genhke  jnd 
heliumlike  satellite  lines  for  carbon,  aluminum,  and  argon  tor 
typical  plasma  conditions. 


T.  eV- 

FKi  *>.  Ionization  abund-inee  curses  »  electron  temperature 
for  carbon.  aluminum,  and  argon  Plasma  dcnsits  has  hern 
fixed  at  Hr  ions/cm 


T.  -«vi 

FIG.  10.  Heliumlike  satellite  line  intensities  normalized  to 
the  resonance  line)  ss  temperature  for  an  aluminum  plasma  ai 
10'  ions  cm 


satellite  to  RES1  are  excellent  indicators  of  temperature  in 
the  plasma.  A  good  temperature  diagnostic,  however.  is 
one  which  is  insensi 1 1 s e  to  densitx  oxer  some  range  of  this 
parameter 

In  Fig.  II.  the  densitv  behavior  of  the  heliumlike  satel¬ 
lite  lines  is  shown  for  an  aluminum  plasma  at  SCO  eV.  At 
lower  densities,  the  double  excited  states  populated  mainly 
by  electron  capture  give  nse  to  satellites  which  arc  some¬ 
what  independent  of  density,  while  the  rt.rn'.  a.r.c.d  . 
and  q.r‘  satellites,  for  which  inner-shell  excitation  is  im¬ 
portant.  show  a  distinct  drop  as  the  density  approaches 
lip"  tons. 'cm-.  It  will  be  shown  later  that  the  inner-shell 
processes  become  less  important  in  determining  satellite 


I  Hi  11  Hehumlike  s.iH-iate  lire  miensines  normalized 
the  resonance  line  \\  nn  densits  tor  an  -ilummum  plasma  ,o 
8(10  eV 
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intensities  as  the  plasma  density  increases.  All  of  the 
satellites  display  a  rise  in  intensity  relative  to  the  reso¬ 
nance  line  as  the  density  increases  from  10*‘;  tons/em  . 
This  is  due  to  the  fact  that  the  doubly  excited  levels  are 
approaching  a  state  of  local  thermodynamic  equilibrium6 
1LTE1  with  the  h:  ground  state,  and  the  intensities  exhib¬ 
it  mild  increases  from  their  coronal  values. 

For  the  hydrogenlike  satellite  lines,  shown  in  Fig.  12,  a 
temperature  behavior  similar  to  that  of  the  heliumlike 
satellites  is  seen  The  values  for  the  satellite  to  Lyman-a 
ratios  shown  are  for  an  aluminum  plasma  at  a  density  of 
10*°  ions/cm3  The  line  emanating  from  the  2p2 1 D  level 
is  the  strongest  emitter  and  can  be  a  valuable  temperature 
diagnostic  in  dense  plasmas.  As  in  the  figures  for  the 
heliumlike  satellites,  we  group  the  individual  lines  accord¬ 
ing  to  their  nlisL  1  upper  levels,  since  an  LS  representation 
is  still  appropriate  for  A1  XII.  For  argon,  however,  this  is 
not  the  case  (refer  to  Fig.  8!,  as  j  components  from  dif¬ 
ferent  nL  levels  are  blended  in  the  spectrum.  The  LS 
grouping  of  the  h 2p  }P—2p~  P  and  \sls 'S—2s2p  }P 
transitions  is  more  easily  discernible  in  the  spectrum  as 
the  1s2j  'S—2s2p  'P:  transition,  with  the  six  components 
of  the  \slp  ' P—2p-  'P  transition  and  the  two  remaining 
components  of  the  \s2s  xS—2s2p 'P  transitions  i'S;- 
P i  —  j5,-3P0)  grouped  as  a  single  spectral  feature. 

The  density  dependence  of  the  aluminum  hydrogenlike 
satellites  is  shown  in  Fig.  13  for  an  electron  temperature 
of  800  eV.  These  satellite  intensities  are  drifting  upward 
toward  LTE  values  as  was  seen  in  Fig.  1 1  for  the  helium- 
like  satellites  but  the  onset  of  the  departure  from  coronal 
values  takes  place  at  slightly  higher  densities.  The  satel¬ 
lite  due  to  the  \s2p'P-2p:'P  transition  shows  very 
unusual  behavior  with  density;  not  coincidentally,  the 


Itf.  v>;  -yv 


Fki  12.  Hvdrogenukc  sairlhie  ii:v-  : --. to*--! i iv -  n'-rraii/..  ;  to 
the  Lyn.in-c  line  vs  temperature  :  ulunumitn  pia-m.  i: 

10‘"  ions/cm  tor  key,  see  Fig 


FIG.  13.  Hydrogenlike  satellite  line  intensities  normalized  :o 
the  Lyman-a  line’  vs  ion  density  !  an  aiuminum  p.asma  ai  spo 
eV  Tor  key.  see  Fig  ’ 

2p~'P  level  is  populated  predominantly  by  inner-shell  ex¬ 
citation  as  opposed  to  electron  capture,  and  the  unusual 
trend  of  this  curve  will  be  discussed  later  in  the  section  on 
inner-shell  effects  iSec.  Ill  D  . 

C.  Z  scaling  of  satellite  intensities 

An  analysis  of  carbon  and  argon  satellite  emission, 
similar  to  that  for  aluminum  as  shown  in  Figs  10- i  3. 
was  also  completed.  However,  after  inspecting  the 
behavior  of  the  emission  with  T,  and  ,V/.  no  new  nr 
unusually  striking  behavior  was  found,  and  it  was  decided 
not  to  include  those  plots  in  this  study.  Instead,  we  defer 
to  a  plot  of  the  satellite  emission  as  a  function  of 
shown  in  Fig.  14  for  plasmas  at  a  density  of  10* 


Z  TOMlC  NUMIBr'-  2 

1  |.  i  14  11.  Jr  vn  ike  and  helnmlike  salcilit,  nerr  aneed 

line  'ruuMUev  a.  .■  run.:  <r:  .'!  /  at  an  ion  density  or  !f 
n  ils  err-  temper  mire  -  trial  a;  which  the  one-  and  two- 
elec  I  roil  ions  h.ase  canal  r>  p.iiaiions  -  see  Fig.  ° 
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ion  s/cm  .  In  comparing  the  satellite  emission  from  dif¬ 
ferent  materials,  a  choice  of  appropriate  electron  tempera¬ 
ture  for  each  material  had  to  be  made.  For  the  compar¬ 
ison  in  Fig.  14,  we  used  the  temperature  at  which  the  heli- 
umhke  and  hydrogenlike  ion  population  densities  were 
equal  at  a  plasma  density  of  10:u  lons/cm3;  from  Fig.  9. 
this  can  be  seen  to  be  50  eV  for  carbon,  650  eV  for  alumi¬ 
num,  and  2200  eV  for  argon.  Once  again,  the  satellite  in¬ 
tensities  in  Fig.  14  are  normalized  to  their  respective  reso¬ 
nance  line  intensity. 

An  explanation  as  to  the  detailed  behavior  of  each  line 
with  Z  is  quite  complex,  since  competing  collisional  and 
radiative  mechanisms  in  both  the  singly  and  doubly  excit¬ 
ed  states  are  strongly  intermingled  effects  controlling  the 
line  intensities.  At  10"°  ions/cm3,  for  example,  many  of 
the  aluminum  and  argon  satellite  intensities  are  still  close 
to  their  coronal  values,  while  for  carbon,  they  have  al¬ 
ready  begun  to  tend  toward  the  LTE  values  at  this  densi¬ 
ty.  This  is  particularly  noticeable  in  the  {a.b.c.d  I  and 
\s2p'P-2pl  P  intensities,  which  exhibit  a  drop  from  C 
to  A1  and  a  subsequent  rise  from  A1  to  Ar.  As  seen  in 
Figs.  11  and  13,  these  two  satellites  are  very  density  sensi¬ 
tive;  for  Z  =6,  the  effect  of  collisional  coupling  to  nearby 
doubly  excited  levels  has  begun  to  populate  the  \s2p::p 
and  2 p2lP  levels  at  10-°  ions/cnr  and  their  intensities 
have  riven  above  coronal  values,  only  to  drop  with  increas¬ 
ing  Z.  as  the  coronal  regime  is  reentered.  Having  only 
three  data  points  to  plot,  it's  not  clear  whether  these  ratios 
continue  to  rise  or  flatten  out  in  the  vicinity  of  Z  =  18: 
another  set  of  values,  say  at  iron  (Z  =  26l,  would  answer 
the  question. 

The  general  trend  in  the  intensity  ratios  from  AI  to  Ar 
can  be  estimated  using  a  simple  coronal  scheme  including 
only  one  autoiomzing  level,  the  next-highest  ground  state, 
and  one  excited  level  to  that  ground  state  (labeled  by  d,  g, 
and  x).  Solving  the  simple  three-level  system  yields  the 
equations 


A  A'j/Ij  X^CAj/fAj  +  T )  CAd 
J^~.XxAx~  NtX  ~  X<Ad  +  D  ’  * 

where  /,  and  IR  are  the  satellite  and  resonance  line  inten¬ 
sities.  Xd  and  .V,  are  the  doubly  and  singly  excited  level 
intensities,  A  is  radiative  decay.  T  is  autoionization,  C  is 
electron  capture,  and  .V  is  the  impact  excitation  rate  be¬ 
tween  g  and  ,x.  The  quantity  C  is  proportional  to  T.  so  the 
expression  becomes 


h  I  \  PA,  j 
Tr~  X  !  •  r  -  I ' 


For  most  of  the  satellite  transitions,  the  term  in  the 
large  parentheses  increases  with  Z  and  the  normalized  in¬ 
tensities  also  exhibit  an  increase  with  Z.  In  some  cases, 
however,  :hc  trend  is  decreasing  <q,r)  or  only  slight  1\  in¬ 
creasing  .j  —  3,131,  and  the  slight  increase  of  .V'  V,  I  with  Z 
iat  the  temperatures  used  in  Fig.  14>  causes  the  satellite  in¬ 
tensities  to  drop. 


D.  Impact  excitation  and  ionization  effects 

One  important  aspect  of  the  detailed  collisional  model 
employed  here  is  that  it  allows  us  to  investigate  the  role  of 
the  various  collisional  processes  in  determining  the  satel¬ 
lite  intensities.  At  high  plasma  density,  Jacobs  and 
Blaha4  showed  that  the  angular-momentum-changing 
collisions  between  doubly  excited  levels  have  a  strong  ef¬ 
fect  on  the  intensities.  Bhalla  et  a!.'1  also  showed  the  im¬ 
portance  of  inner-shell  excitation  on  certain  heliumlike  sa¬ 
tellite  lines  while  Lunnev  and  Seely31  showed  that  inner- 
shell  excitation  of  hydrogenlike  satellites  can  affect  the 
trtplet-to-singlet  density-sensitive  line  ratio.  In  this  study, 
we  show  the  direct  effects  of  omitting  certain  collisional 
processes  on  the  intensities  of  individual  satellites  of  both 
systems,  including  the  neglect  of  collisional  ionization  and 
recombination  to  n  =2  levels  as  well  as  inner-shell  effects. 
In  Fig.  15,  we  have  plotted  the  (normalized)  satellite  line 
intensities  versus  density  for  several  satellite  lines  in  an 
800-eV  aluminum  plasma.  The  solid  curves  are  those  cor¬ 
responding  to  the  complete  collisional  model  and  mimic 
those  shown  earlier  in  Figs.  1 1  and  13.  The  dashed  curves 
are  the  result  of  setting  the  electron  impact  ionization  and 
collisional  recombination  rate  coefficients  to  zero  in  the 
calculation.  Although  only  four  satellite  intensities  are 
shown  in  the  figure,  the  effect  was  identical  for  all  lines 
from  a  given  ion:  Including  the  collisional  couplings  to 
n  =2  levels  of  the  next-highest  charge  state  brings  the  rel¬ 
ative  satellite  intensities  down  at  high  plasma  densities. 


18  13  20  2  :  22  23 
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FIG.  15.  Selected  'normalized)  satellite  line  intensities  for  an 
aluminum  plasma  at  800  eV.  Comparison  shows  the  complete 
collis.inal  model  'solid),  neglecting  collisional  ionization  and 
recombination  idashedi.  and  neglecting  inner-shell  impact  exci¬ 
tation  and  deexcitation  idashed-dotted). 


65 


28 


DENSE  PLASMA  EFFECTS  ON  K-SHELL  DIELECTRONIC 


2977 


For  aluminum  the  onset  of  this  effect  is  between  10‘ '  and 
I022  ions/cm5,  although  it  occurs  at  lower  densities  for 
lower-Z  materials  as  was  encountered  with  the  carbon 
satellites.  It  also  occurs  at  a  lower  density  for  the  helium¬ 
like  satellites  than  for  the  hydrogenlike  satellites,  as  seen 
in  the  figure.  In  general,  the  density  threshold  for  these 
collisional  effects  to  manifest  themselves  in  the  line  ratios 
is  higher  than  the  range  over  which  the  satellite  diagnos¬ 
tics  are  usually  applied.  However,  caution  must  be  taken 
when  dealing  with  superdense  plasmas  such  as  laser- 
imploded  microballoon  targets,  since  these  densities  have 
been  attained 

The  next  effect  of  interest  was  the  impact  of  inner-shell 
electron  collisions  on  the  satellite  intensities.  For  this 
study,  the  inner-shell  excitation  and  deexcitation  rate  coef¬ 
ficients  were  set  to  zero  in  the  model.  The  results  are 
shown  as  dashed-dotted  curves  in  Fig.  15  Only  four 
satellites  are  shown  since  the  remaining  satellite  intensities 
suffered  less  than  a  10 %  perturbation  due  to  the  neglect  of 
the  mnershel!  rates  over  the  density  range  shown.  The  de¬ 
gree  to  which  the  innershell  processes  are  important  can 
be  measured  by  the  difference  between  these  curves  and 
those  of  the  complete  model  (solid  lines).  The  i q.ri  satel¬ 
lite  curves  indicate  that  inner-shell  processes  are  probably 
comparable  to  electron  capture  in  populating  the 
ls2s('P)2p  level  while  inner-shell  excitation  is  dominant 
for  the  Is2p-2S  ( n.m>.  1  s2p~~P  'a.b.c.a  i.  and  2p-'P 
(4—9)  levels.  The  universal  trend  that  is  evident  in  the  fig¬ 
ure  is  that  an  essential  populating  mechanism  for  certain 
levels  is  missing  when  inner-shell  excitation  is  neglected  in 
a  collisional  model  and  the  corresponding  satellite  lines 
are  much  reduced  in  intensity.  Note,  however,  that 
inner-shell  processes  play  a  much-diminished  role  in  the 
ionization  balance  as  the  plasma  density  is  increased.  The 
explanation  for  this  is  as  follows  At  low  densities,  ihe 
doubly  excited  levels  are  collisionally  uncoupled  from  each 
other  and  levels  which  have  small  electron  capture  rates 
depend  on  inner-shell  excitation  to  populate  them  As  the 
density  increases,  however,  the  cross-coupling  rates  come 
into  play,  and  the  levels  which  siepend  strongly  on  the 
inner-shell  processes  are  now  populated  by  electron  col¬ 
lisions  from  the  autoiomzing  states.  This  alternate 
creation  mechanism  reduces  the  role  of  inner-shell  excita¬ 
tion  at  these  higher  densities,  and  the  curve  which  doesn't 
include  inner-shell  effects  approaches  the  solid  curve. 

The  fact  that  these  satellites  have  small  autoionizution 
rates  'and,  hence,  small  electron  capture  rates,  compared 
to  the  inner-shell  excitation  rates'  is  what  makes  them  ex¬ 
cellent  diagnostic  indicators  of  density  .  For  example,  let 
us  examine  the  behavior  of  the  > a.b.c.a '  satellite  lines  in 
Fig.  15,  which  display  an  unusual  double-valued  nature. 
At  10'“  ions/cm',  inner-shell  excitation  via  the 
is~2p  :P-  la  2p~  'P  transition  populates  ihe  upper  level, 
and  the  satellite  lines  are  ten  times  more  intense  than  if 
this  process  was  neglected.  At  about  5  •  1 0 1 "  ions  cm', 
the  lithiumlike  n  =2  state  densities  depart  from  their  co¬ 
ronal  values  and  begin  to  drop  with  density  as  iom/ation 
processes  shift  population  to  higher  ions.  As  the  1  s’2p 
level  depletes,  the  \s2p'~P  level  follows  and  the  satellite 
intensity  drops.  As  the  density  increases  further,  the  col- 
lisiona)  coupling  between  doubly  excited  states  becomes 


significant  (at  about  102°  ions,  cm)1  and  the  1  s2p':P  state 
is  repopulated;  hence,  the  ia.b.c.d)  satellite  intensity  in¬ 
creases.  A  more  detailed  discussion  of  these  satellites  in 
relation  to  density  diagnostics  follows  in  Sec.  Ill  E. 

Two  final  points  should  be  emphasized  regarding 
inner-shell  processes  and  these  satellite  lines.  Inner-shell 
excitation  was  found  to  be  important  in  populating  only  a 
few  of  the  levels  in  an  equilibrium  plasma.  In  a  transient¬ 
ly  ionizing  plasma,  however,  where  the  ionization  tem¬ 
perature  significantly  lags  behind  the  electron  tempera¬ 
ture,  electron  capture  may  be  dominated  by  inner-shell  im¬ 
pact  excitation  for  all  the  doubly  excited  levels,  due  to  the 
reduced  ground-state  population  available  for  the  dielect- 
ronic  recombination  process.-'0  Also,  collisional  deexetta- 
tion  from  doubly  to  singly  excited  levels  does  not  seem  to 
play  a  major  role  in  the  determination  of  the  satellite  line 
intensities  over  the  useful  diagnostic  temperature  range 
until  electron  densities  of  about  lO24  cm--  are  attained. 
For  example,  at  \e  =  102'’  cm'-'  and  T,  =  800  eV,  the 
neglect  of  inner-shell  deexcitation  in  an  aluminum  plasma 
resulted  in  only  a  3. 592  error  in  the  calculat¬ 
ed  t Is  2s  }S—2s2p  P-p  Is 2 p  P-2p::'PV t Is 2 p  'P-2p:  'D 
line  ratio. 

E.  Density  diagnostics 

As  stated  earlier,  several  of  the  satellite  lines  emanating 
from  levels  which  are  not  strongly  autoionizing  are  densi¬ 
ty  sensitive.  In  addition,  the  density  sensitivity  of  certain 
satellites  at  high-density  can  color  other  diagnostics,  in 
particular  the  well-known  intercombination-to-resonance 
line  ratio  {•  ls:-ls2p  ’PYi  ls:-\s2p  This  ratio  is  reg¬ 

ularly  exploited  for  density  determinations  in  laboratory 
plasmas,  despite  the  fact  that  several  heliumlike  satellite 
lines  are  usually  blended  unresolvably  with  the  1C  line. 
From  Fig.  8  one  can  see,  for  example,  that  the  im.nl  satel¬ 
lites  overlap  the  carbon  IC  line,  the  tm.ni  and  :s.t>  satel¬ 
lites  overlap  the  aluminum  1C  line,  and  the  ;m  and  'S.t 
satellites  overlap  the  argon  IC  line  Although  these  satel¬ 
lites  are  not  intense  at  lower  densities,  their  density  depen¬ 
dence  can  influence  this  ratio  as  .V,  increases.  Ilyukhin 
et  at  "  realized  tins  fact  and  attempted  to  correct  for  the 
effect  by  scaling  ihe  blended  sateilile  intensities  with  the 
resolvable  'j.k  line-..  As  an  example  of  this  effect,  the 
IC-to-RES  ratio  :«  shown  in  Fig.  16  for  an  S00-eV  alumi¬ 
num  ano  a  !6U>-eV  argon  plasma.  The  dashed  curve  is 
ihe  IC-to-RES  intensity  ratio  alone,  while  the  solid  curve 
includes  the  1  satellites  for  aluminum  and  the 

s.t. in  satellites  for  argon.  The  largest  contribution  is 
from  the  im  -ateihte.  while  the  >s.n  and  ut‘  are  each 
about  half  the  intensity  of  the  uni  line.  As  the  figure 
shows,  serious  errors  in  the  ion  ‘electron1  density  deter¬ 
mination  can  result  from  ignoring  the  satellite  contribu¬ 
tion  to  the  total  line  intensity  in  ihe  vicinity  of  ihe  IC  line. 
It  will  be  shown  in  a  following  paper  that  self-absorption 
in  an  optically  thick  plasma  further  enhances  the  differ¬ 
ences  seen  in  Fig.  lb. 

Another  density-sensitive  line  among  the  heliumlike 
satellites  is  the  'a.b.c.a  >  biend  as  first  pointed  out  by 
Jacobs  and  Blaha.4  The  (normalized'  intensity  of  this 
feature  is  shown  versus  ion  density  tn  Fig.  17,  for  a  60-eV 
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FIG.  16.  IC-to-RES  line  intensity  ratio  vs  ion  density  for  an 
800-eV  aluminum  plasma  and  a  1600-eV  argon  plasma  (dashed1. 
Also  shown  is  the  IC  +  j,r,m,n-to-RES  ratio  for  aluminum  and 
the  IC  +  s,i, m-to-RES  ratio  for  argon  (solid ). 


carbon,  800-eV  aluminum,  and  1600-eV  argon  plasma  [for 
argon,  the  resolvable  feature  is  the  ( a,d,r )  blend].  Over 
the  10lg— 1023  ion/cm5  range,  the  aluminum  and  argon 
line  ratios  are  double-valued,  as  was  discussed  in  Sec. 
IIID.  Care  must  be  taken,  therefore,  to  ascertain  on 
which  side  of  the  minimum  the  density  lies  when  making 
density  determinations  via  this  diagnostic.  For  carbon, 
however  these  lines  may  be  the  most  valuable  diagnostic 
among  the  heliumlike  lines,  since  the  IC  line  is  usually  too 
weak  to  be  detected  in  most  carbon  spectra.  This  intensity 
ratio  may  also  be  of  special  importance  in  laser-fusion  ex¬ 
periments  which  use  argon  as  a  seed  gas  for  diagnostic 
purposes.  Since  these  heliumlike  satellites  occur  at  lower 
plasma  temperatures  than  do  the  hydrogenlike  satellites, 
they  may  be  more  intense  in  the  dense,  imploded  core 
where  temperatures  are  lower  than  in  the  corona. 

Finally,  the  density-sensitive  3P  hydrogenlike  satellites 
have  recieved  much  attention  from  previous  workers”'1' 
in  recent  years,  particularly  because  of  their  usefulness  in 
laser-target  compression  experiments.'0  In  Fig.  18,  the 
( ls2s  }S  —  2s2p  }P)  +  (\s2p  ~P  —  2p2  }P)/l  Is  2p  ]P-2p:lD\ 
line  intensity  ratio  versus  ion  density  is  shown  for  carbon 
(60  eVl,  aluminum  (800  eV),  and  argon  (1600  eV>.  Also 
shown  are  the  results  of  neglecting  collisional  ionization 
and  recombination  to  the  n  =2  hydrogenic  level  (dashed 


«jg  •«  err' 

I  I(  1  I"  1  ,d  sateliife-M-rcsnnancc  line  intensity  ratio  is 
i*  n  density  loi  a  nit-eV  earhon.  stSt-cV  aluminum,  and  a  It  .0. 
cV  argon  pla.tna  ■the  argon  ratio  is  the  u-d-r-to-RES  lines 


FIG.  18.  Total  "triplet"  to  "stnglet-D"  hvdrogenlike  satellite 
line  intensity  ratio  vs  ion  density  for  the  same  temperatures  as  in 
Fig.  17.  Comparison  shows  the  complete  collisional  model 
(solid),  neglecting  collisional  ionization  and  recombination 
(dashed),  and  neglecting  inner-shell  impact  excitation  and  deex¬ 
citation  (dashed-dotted). 


curvet  and  neglecting  inner-shell  impact  excitation  and 
deexcitation  to  the  Is  2s  and  ls2p  levels  (dashed-dotted 
curves).  The  dashed-dotted  curve  for  carbon  is  not  shown 
since  the  results  vaned  by  less  than  a  few  percent  from  the 
complete  model.  As  seen  from  the  various  curves,  omit¬ 
ting  either  of  the  collisional  processes  results  in  no  more 
than  a  factor  of  2—2  error  in  the  density  determination  in 
the  region  where  the  ratio  is  a  sensitive  diagnostic.  The 
noticeable  bump  in  the  aluminum  and  argon  ratios  is  the 
effect  of  inner-shell  excitation  for  the  \s2p}P—2p'  P 
transition.  It  should  be  noted  once  again,  that  these  are 
results  from  an  equilibrium  model,  in  which  the  atomic 
processes  occur  on  a  rapid  enough  time  scale  that  the  ioni¬ 
zation  balance  reflects  the  local  plasma  electron  tempera¬ 
ture.  The  study  of  satellite  lines  in  a  transient  plasma  will 
be  the  subject  of  a  later  investigation. 


IV.  DISCUSSION 

To  summarize,  a  new  model  has  been  described  for 
characterizing  the  emission  from  the  n—2  dielectronic 
satellite  lines  to  the  heliumlike  and  hydrogenhke  reso¬ 
nance  lines  in  dense  plasmas.  In  this  model,  the  doubly 
excited  level  populations  are  solved  for  self -consistently 
with  a  number  of  singly  excited  and  ground  levels  for  a 
plasma  in  collisional-radiative  equilibrium.  In  addition 
great  care  has  been  taken  to  calculate  rate  coefficients  for 
all  the  relevant  collisionally  induced  and  spontaneous  pro¬ 
cesses  to  high  accuracy  using  sophisticated  atomic  stru,.- 
tur  .Hid  scaiterme  modeis. 

ice  behavior  n!  the  satellite  lines  has  been  studied  fa 
the  ptically  thin  plasma"  approximation  and  their  van- 
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ation  with  temperature  and  density  described.  The  effects 
of  neglecting  inner-shell  processes  and  collisional  ioniza¬ 
tion  were  shown  to  cause  large  errors  in  theoretically  cal¬ 
culated  intensities  in  some  density  regimes  and  to  have  a 
negligible  effect  in  others.  Density  determinations  using 
the  heliumlike  intercombination-to-resonance  line  ratio 
were  found  to  be  strongly  impacted  by  unresolvable  satel¬ 
lites  at  higher  densities.  The  effects  were  shown  to  com¬ 
mence  at  densities  of  about  10;i  e  7cm',  corresponding  to 
the  critical  density  for  1. 06-jam  laser  light,  and  to  become 
more  important  as  the  density  increased,  causing  as  much 
as  2  orders  of  magnitude  error  in  the  measured  plasma 
density.  The  ia—d)  satellite  lines  are  also  shown  to  be  a 
valuable  density  diagnostic  at  temperatures  lower  than 
those  where  the  density-sensitive  hydrogenlike  satellites 
can  be  detectable.  This  diagnostic  is  of  special  importance 
for  low-Z  plasmas  where  the  intercombination  line  is 
often  too  weak  to  be  detected.  Finally,  we  investigated  the 
ratio  of  the  triplet  satellite  lines  to  the  'D  satellite  adjacent 
to  the  Lyman-a  line,  and  found  that  inner-shell  and/or 
ionizing  collisions  can  affect  the  density  determination  by 
factors  of  2  to  3  in  the  useful  density  and  temperature 
range  of  the  diagnostic. 

In  conclusion,  there  are  two  additional  points  which 
should  be  mentioned.  The  first  regards  a  recent  paper  0 
dealing  with  dielectronic  capture  followed  by  ladderlike 
excitation  and  ionization  in  dense  plasmas.  The  con¬ 
clusion  arrived  at  in  this  paper  is  that  the  effective  col¬ 
lisional  ionization  rate  is  significantly  enhanced  when  one 
includes  a  large  number  (in  this  case,  up  to  principal 
quantum  number  201  of  the  higher  doubly  excited  levels', 
thus,  the  actual  ionization  state  in  a  dense  plasma  can  be 


larger  than  one  would  calculate  without  inclusion  of  these 
levels.  That  this  investigation  arrives  at  such  a  result  is 
not  too  surprising  since  excitation  rates  often  exceed  ioni¬ 
zation  rates  by  more  than  an  order  of  magnitude  In  fact, 
the  same  conclusion  was  reached  in  an  earlier  study’’5  re¬ 
garding  ionization  via  excitation  of  singly  excited  levels  in 
attaining  LTE  conditions  However,  no  mention  is  made 
in  their  paper  of  continuum  lowering  'where  pressure  ioni¬ 
zation  may  force  the  outer  electrons  of  many  of  these  dou¬ 
bly  excited  states  into  the  continuum).  The  point  here  is 
that,  in  a  plasma  at  high  electron  densities,  it  may  be  diffi¬ 
cult  to  assess  the  ladderlike  processes  without  self- 
consistently  including  the  effects  of  ionization  lowering 
and  level  shifts,  widths,  and  merging  land,  hence,  Stark 
broadening)  in  the  model. 

Finally,  in  this  study  we  have  neglected  the  plasma  mi¬ 
crofield  effects  on  the  autoiomzation-dielectromc  capture 
processes.  In  addition,  no  attempt  was  made  to  take  ac¬ 
count  of  the  effects  of  strongly  overlapping  lines  1 "  '  or 
of  the  true  Funo-Beutler  profile  4  in  the  calculations 
Both  of  these  effects,  as  well  as  photon  self-absorption 
and  nonequihbnum  plasma  effects,  are  presently  under  in¬ 
vestigation  and  will  be  reported  on  later. 
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D.  Dielectroale  Recombination  of  Krypton  Ions 


I.  INTRODUCTION 

The  rate  coefficients  for  radiative  and  dielectronic  recombination  are 
needed  to  predict  the  relative  abundances  of  various  ions  in  a  high- 
temperature  plasma.  These  relative  abundances  ra ay  be  used  to  calculate  the 
far-ultraviolet  and  X-ray  emission  spectra  produced  by  electron-ion 
collisions.  In  the  corona  model  approximation  (Griera  1964),  which  is  valid 
at  low-densities,  the  recombination  coefficients  are  functions  only  of  the 
local  electron  temperature  and  are  independent  of  the  electron  density. 
This  approximation  is  expected  to  be  valid  for  the  K-shell  ions  of  a  high-Z 
element  over  an  extensive  density  range.  However,  the  coefficients  for  the 
L-  and  M-shell  ions  are  expected  to  be  much  more  sensitive  to  density 
effects . 

Burgess  (1964)  first  demonstrated  that  dielectronic  recombination  is 
often  more  important  than  the  usual  free-bound  radiative  process  in  the 
temperature  region  where  the  ions  have  their  maximum  abundance  in 
ionization  equilibrium.  A  simple  formula  was  derived  by  Burgess  (1965)  for 
estimating  the  dielectronic  recombination  rates  in  the  corona  model 
approximation.  This  simple  formula  requires  knowledge  of  only  the  energy 
levels  and  absorption  oscillator  strengths  of  the  recombining  ion. 
However,  ionization  -  equilibrium  calculations  which  have  been  carried  out 
using  the  Burgess  formula  (Jordan  1969)  do  not  agree  with  more  detailed 
calculations  (Jacobs  et.  al.  1977)  in  which  the  rates  for  all  radiative 
decay  an  autoionization  processes  from  the  doubly-excited  states  are  taken 
into  account.  This  is  attributable  to  the  discovery  that,  for  certain 
ions,  autoionizat ton  into  an  excited  state  of  the  recombining  ion  Is  more 
probable  than  autoionizatlon  to  the  ground  state. 

II.  THEORY 

Dielectronic  recombination  nay  be  described  as  a  two-step  process. 
Radiationless  capture  of  a  plasma  electron  Into  an  nl-state  accompanied  by 
the  excitation  1+j  of  the  recombining  ion  results  in  the  formation  of 
a  doubly-excited  state 
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X+(Z)  (i)  +  e~  >  x+(z-1)  ( j  ,nl) 


Recombination  occurs  if,  instead  of  autoionization,  the  doubly-excited 
state  j,  nl  undergoes  radiative  decay  to  a  sin; ly-excited  state  i,  nl  which 
lies  below  the  ionization  threshold 


X+(Z_1)  (j.nl)  -►  X*2"1*  (i ,nl)  +  hu 


The  initial  state  i  is  assumed  to  be  the  ground  state  in  the  corona  model 
approximation.  In  our  calculations,  the  many-electron  states  1  and  j  will 
be  specified  by  giving  only  the  principal  and  angular  momentum  quantum 
numbers  of  the  active  electron. 

Assuming  that  the  electrons  have  a  maxwellian  velocity  distribution, 
the  dielectronic  recombination  rate  coefficient  in  the  corona  model 
approximation  is  given  by  (Shore  1969) 

,3  3  3/2  ,,  .  3/2 
ad(i)  -  2  aQ  *  (VkBTe) 


S  A  (j,nl+i)  A  (  j  ,nl-i ,nl) 

Y  r-  r  8(  j.nl)  _a _ r _ 

A  3  hi  2g(i)  Aa(J,nl)  +  Ar(j,nl) 


X  exp  {  E(i)kV(J>?1)  1  *  (3) 

B  e 

where  E(i)  and  E(j,nl)  are  the  initial-state  and  doubly-excited  state 
energy  levels,  respectively,  and  g(i)  and  g(j,nl)  are  their  statistical 
weights. 

III.  CALCULATIONS 

For  large  values  of  n  the  autoionization  rates  Aa(j,nl-ki)  can  be 
obtained  from  the  threshold  values  of  the  partial-wave  cross  section  for 
the  electron  Impact  excitation 


X+(Z)  (1)  +  e“  (kjlj)  >  X+(Z)  (j)  +  e“  (k  lj)  ,  (4) 
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by  using  the  quantum-defeat  theory  relationship  derived  by  Seaton  (1969). 
The  distorted-wave  method  (Davis  et.al.  1976)  was  used  to  obtain  the 
excitation  cross  sections.  The  stabilizing  radiative  decay  rates  Ar(j,nl 
>i,nl)  where  approximated  by  the  spontaneous  emission  rate  Af(j>i)  for  the 
recombining  ion  core.  These  approximations  are  expected  to  be  valid  for 
An^  *  0  transitions  of  the  ion  core,  which  occur  predominantly  by  means  of 
large  n-values.  They  are  uncertain  for  An^not  equal  to  0  transitions,  for 
which  small  n-values  play  a  more  important  role,  especially  with  increasing 
Z.  The  total  decay  rates  Aa(j,nl)  and  Ar(j,nl)  include  the  rates  for  all 
allowed  autoionization  and  radiative  decay  processes.  For  some  An^  not 
equal  to  0  transitions,  autoionization  into  an  excited  state  of  the 
recombining  ion  can  occur  and  is  found  to  be  more  probable  than  auto¬ 
ionization  into  the  ground  state.  This  gives  a  substantially  smaller 

recombination  rate  than  predicted  by  the  Burgess  formula. 

The  dielectronic  recombination  rate  coefficients  for  Kr  XXVI  -  XXXVI 
have  been  calculated,  taking  into  account  the  autoionization  processes  and 
stabilizing  radiative  transitions  which  involve  a  single-electron  electric- 
dipole  transition  of  the  recombining  ion  core.  The  transitions  Included 
are  given  in  Table  I.  The  asterisk  has  been  used  to  identify  the 
transitions  whose  contribution  to  is  reduced  by  autoionization  into  an 

excited  state. 

The  total  dielectronic  recombination  coefficient  for  each  ion  is  given 
in  Table  II  for  the  temperature  region  of  experimental  interest.  The  K- 
shell  ions  are  expected  to  have  density-independent  recombination 
coefficients,  even  at  densities  as  high  as  10  *  cm  .  The  results  for  the 
L-  and  M-shell  ions  may  not  be  valid  at  such  high  densities  because  of 
collisional  processes  Involving  the  outer  nl-electron.  To  obtain  the  total 
recombination  coefficient  in  the  corona  model  approximation,  the  direct 
radiative  recombination  coefficients  must  be  added  to  the  results  in  Table 
II.  Estimates  based  on  hydrogenic  approximations  indicate  that  for  these 
Kr  ions  radiative  recombination  is  almost  as  important  as  dielectronic 
recombination.  The  is  to  be  expected  from  Z-scaling  arguments. 


/ 
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Table  1 


Stabilizing  Radiative  Transitions  j  -*•  i 


Recombining 

Single-Electron 

ion 

Transitions  j  -*■  i 

Kr  XXXVI 

2p* Is , 

3p>  Is 

Kr  XXXV 

2p+ls , 

3p+  Is 

Kr  XXXIV 

2p*2s  , 

3p>  2s* 

Kr  XXXIII 

2p*2s , 

3p+2s* 

Kr  XXXII 

2p+2s , 

3s+2p,  3d-*-  2p* 

Kr  XXXI 

2p+2s , 

3s+2p,  3d-*-2p* 

Kr  XXX 

2p-2s, 

3s-*-2p,  3d--2p* 

Kr  XXIX 

2p+2s , 

3s+2p,  3d+2p* 

Kr  XXVIII 

2p+2s , 

3s+2p,  3d+2p* 

Kr  XXVII 

3s-*-2p, 

3d-*-2p* 

Kr  XXVI 

3p>3s , 

4p-*-  3s* 

*Affected  by  autoionization  into  an  excited  state 


Table  II 

Dielectronic  Recombination  Rat 
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E .  CoupLed  Plasma  and  Radiation  Transport  Ion  Beam  Coda 


I.  INTRODUCTION 

In  order  Co  understand  the  nature  of  the  dynamic  material  response  of 
a  planar  target  exposed  to  a  high  energy  ion  beam,  particularly  in  the 
analysis  and  assessment  of  directed  energy  and/or  particle  beam  lethality 
and  vulnerability  experiments,  a  one-dimensional  fully  self-  consistent 
energy  deposition/stopping  power  radiation  hyd rodynamics  model  has  been 
developed . 

In  Section  II  the  physical  processes  included  in  the  model  are 
detailed,  along  with  the  numerical  formulation  employed.  Section  III 
presents  the  basic  results  of  our  treatment  of  the  interaction  of  a  1  MeV 
proton  beam  with  a  planar  aluminum  target.  Section  IV  compares  these 
results  with  those  obtained  when  radiation  transport  is  omitted  from  the 
calculations.  Section  V  compares  our  results  with  a  similar  laser-target 
interaction  study.  Finally,  in  Section  VI,  the  work  is  summarized. 

II.  THEORETICAL  MODEL 

The  interaction  of  the  ion  beam  with  the  target  is  strongly  dependent 
on  the  local  temperature  and  degree  of  ionization  in  the  target.  In 

addition,  optical  pumping  and  energy  transport  by  photons  can  influence  the 
plasma  temperature  and  degree  of  ionization,  and  can  modify  the 

hydrodynamic  response  of  the  target.  Thus,  the  deposition  of  the  beam,  the 
hydrodynamic  evolution  and  atomic  physics  of  the  target,  as  well  as  the 
transport  of  radiation,  must  be  calculated  self-consistentlv . 

Discussion  of  the  theoretical  model  can  be  separated  for  convenience 
as  follows:  (a)  hydrodynamics  and  thermal  conduction  (b)  ionization  and 

atomic  physics  (c)  radiation  emission  and  transport  and  (a)  beam  dynamics 
and  energy  deposition. 

(A)  Hydrodynamics  and  Thermal  Conduction 

The  basic  hydrodynamic  variables  of  mass,  momentum,  and  total  energy 
are  transported  in  one  dimension  using  a  numerical  scheme  with  a  sliding- 
zone  version  of  flux-corrected  transport.^  A  special  gridding  algorithm  is 


used  which  moves  cones  in  a  Lagrangian  fashion  and  adjusts  the  mesh  in 
order  to  resolve  steep  gradients  in  the  flow.  The  hydrodynamic  equations 
solved  are 


Do 

Dt 


(  u-'  ) 


0, 


(n 


D(cu)  _  _  d? 

Dt  ~  dx  ’ 


(ur>)  + 


(2) 


+ 


(3) 


where  o  is  mass  density,  u  is  velocity,  P  is  pressure,  T  is  total  energy 
density  is  the  rate  of  energy  loss  or  gain  due  to  radiation,  ^  is 

che  rate  of  energy  gain  due  to  the  beam  deposition  process,  -  is  the 
thermal  conductivity,  and  N  is  the  Ion  density.  The  thermal  conduction  is 
calculated  implicitly,  using  an  Iterative  Crank-N’lcholson  scheme. 

Since  densities  did  not  much  exceed  solid  density  in  this  study,  a 
simple  equation  of  state  was  assumed  with 


P 


2 


(A) 


where  is  the  ion  potential  energy  due  to  ionization  and  excitation.  A 
single  temperature  model  was  employed, 


kT  =  _ t _ 

(o/nj)  (1+Z) 


(5) 


where  m-j-  is  ion  mass,  and  T  is  temperature.  The  Ionization  energy,  and 
effective  charge,  Z,  are  calculated  from  the  lonl cat  ton-rad ia t ion  equations 
which  are  explained  he  low.  A  single  temperature  assumption  is  valid  in  the 
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cold  dense  target  as  well  as  in  the  beam  deposition  region,  since  thermal 
equilibration  times  ire  short  compared  with  the  time  scales  of  the  ener ?y 
input  and  changes  in  the  hydrodynamic  variables. 


The  local  rate  of  change  of  energy  due  to  radiation  transport 
•  • 

and  that  due  to  the  beam  deposition  £  ’ e p  >  be  discussed  in  the 

appropriate  sections  below. 


(3)  Ionization  and  Atomic  Physics 


The  ionic  populations  in  the  plasma  may  be  characterized  by  a  set  of 
atomic  rate  equations  of  the  form 


dfi 

dt 


jt 


W.  .  f 

i 


( f>) 


where  f^  is  the  fractional  population  of  atomic  level  i,  and  Is  the  net 
reaction  rate  describing  the  transition  from  initial  state  j  to  final  state 
i.  An  equation  of  this  type  is  constructed  for  each  cf  the  atomic  levels 
included  in  the  model. 

For  sufficiently  dense  plasmas  of  the  sort  we  expect  to  model  with  the 
ion-beara  code,  the  effective  populating  and  depopulating  rates  are 
extremely  fast  compared  with  the  hydrodynamic  response.  An  equilibrium 
assumption  can  be  justified,  which  involves  dropping  the  explicit  time 
dependence  in  equation  (6).  The  plasma  is  then  said  to  be  in  collisional- 
radiative  equilibrium  (ORE)  ,  whereby  the  plasma  ionization  state  responds 
ins tantaneous ly  to  changes  in  hydrodynamic  quantities. 

The  rate  coefficients  that  are  used  to  calculate  the  populating  and 
depopulating  rates,  Wji»  are  calculated  using  various  scattering 
techniques.  The  processes  included  in  this  calculation  and  the  methods 
used  in  calculating  the  corresponding  rate  coefficients  are  summarized 
below. 


( 1 )  Collisional  ionization  -  exchange  classical  impact-parameter  (ECIF) 
methods  (the  effect  of  autolonizing  resonances  on  the  cross  sections  has 
been  ignored). 


hydrogenic  approximation  with  Karzas-Lat ter-Gaunt 


( 2 }  Photo  Lonizat  ton  - 

factors  . ^ 

(3)  Dielectronic  recombination  -  the  detailed  calculations  of  Jacobs  et 
al.^  are  used. 

(4)  Collisional  excitation  -  Coulomb-Born  distorted-wave  approximation 
including  exchange  effects,'  or  the  seraiclassical  impact-parameter  (SCI) 

O 

technique . 

(5)  Spontaneous  radiative  decay  -  oscillator  strengths  are  taken  from 

9 

several  calculations  and  measurements. 

(6)  Photoexcitation  -  oscillator  strengths  used  are  those  quoted  above  to 
determine  optical  depths  (see  the  next  section  on  radiation  transport). 

Finally,  collisional  and  radiative  recombination,  collisional  deexcitation, 
and  stimulated  emission  are  all  calculated  as  the  detailed  balance  of  the 
corresponding  opposite  rate  listed  above.  A  simplified  level  diagram  using 
two  ground  states  and  a  single  excited  state  is  shown  in  Fig.  1  to  show 
schematically  the  processes  Included  in  the  model. 

Once  the  set  of  rate  equations  (including  the  radiation  transport)  has 
been  solved  for  the  level  populations,  f^,  the  electron  density  can  be 
calculated , 


:■<  =  z  f  n 

e  1  i  I 


(7) 


where  z^  is  the  ionic  charge  of  level  i  and  Nj  is  the  total  ion  density. 
The  ionization  and  excitation  energy  can  also  be  calculated  by 


=  '  v  f  N 
I  '■  xi  1  I’ 


(8) 


where  is  the  energy  of  level  i,  measured  ^rora  the  ground  state  of  the 
neutral  atom. 


(C)  Radiation  Emission  and  Transport 

Radiation  emission  from  a  plasma  and  its  opacity  are  dependent  on  the 
local  atomic  level  population  densities.  Except  for  optically  thin 
plasmas,  however,  the  level  populations  depend  on  the  radiation  field, 
since  optical  pumping  via  photolonizat ion  and  photoexcitation  can  produce 
significant  population  redistribution.  Thus,  the  ionization  and  radiation 
transport  processes  are  strongly  coupled  and  must  be  solved  self- 
consistently .  In  this  model,  an  iterative  procedure^®  is  used,  where  level 
populations  are  calculated  using  the  radiation  field  from  the  previous 

Iteration,  then  using  these  populations  to  calculate  a  new  radiation  field 
and  recalculating  populations  until  convergence  is  reached.  Even  when  the 
local  Thermodynamic  Equilibrium  (LTE)  populations  are  used,  an  iterative 

procedure  is  required,  since  the  LTE  populations  are  functions  of 

temperature  and  electron  density,  which,  in  turn,  depend  on  the  level 
populations . 

Three  distinct  radiation  transport  schemes  have  been  developed,  and 
can  be  used  interchangeably  in  the  code:  A  probabilistic  scheme,  a 
raultif requency  scheme  and  a  hybrid  method.  Descriptions  of  these  transport 
methods  have  appeared  elsewhere;  in  this  brief  report,  their  general 
properties  will  be  outlined. 

The  probabilistic  model^  forms  local  angle  and  frequency  averaged 

escape  probabilities  for  each  emission  line  and  for  each  bound-free 
process.  Free-free  radiation  is  treated  with  a  multifrequency  formalism. 
The  radiation  transport  and  emission  spectra  are  calculated  from  these 
escape  probabilities.  The  method  is  cost-effective,  can  treat 

comprehensive  atomic  models  and  provides  good  overall  energetics,  but 
cannot  calculate  certain  spectral  details  and  breaks  down  at  very  high 
densities  where  no  systematic  treatment  of  Stark  profile  escape 
probabilities  is  available. 

The  mult i f requency  model  “  solves  the  equation  of  radiation  transport 
at  a  large  number  of  discrete  frequencies,  providing  resolution  of  emission 
lines,  recombination  edges  and  absorption  edges.  It  provides  accurate 
radiation  transport  at  high  density,  and  gives  spectral  details  such,  as 
self-absorption  features.  However,  a  large  number  of  frequencies  is 

requLred  to  provide  adequate  resolution,  and  it  tends  to  be  a  relatively 
costly  method. 
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The  hybrid  model  uses  the  multifrequency  formalism  to  transport  the 
continuum  (bound-free  jnd  free-free)  radiation,  and  frequency  integrated 
escape  probabilities  to  transport  the  lines.  Continuum  opacities  are 
interpolated  from  the  multifrequency  mesh  and  folded  into  the  line 
transport  calculations.  The  line  opacities  are  assumed  to  have  negligible 
impact  on  the  continuum  energetics. 

In  all  of  the  models  outlined  above,  local  ionization  state-dependent 

inner-shell  opacities  are  included,  since  these  processes  are  very 

important  in  the  cool,  dense  plasma  regions.  Inner-shell  photoionization 

cross  sections  for  the  neutral  element  is  taken  from  the  fits  by  Biggs  and 
1 3 

Lighthill  ,  and  the  positions  of  the  ionization-dependent  absorption  edges 
are  taken  from  the  Hartree-Fock  calculations  of  Clementi  and  Roetti^'4. 

The  local  rate  of  energy  change  in  zone  j,  due  to  radiation  transport 
is  given  by 


t.  . 

J 


(9) 


where  Fpk  is  the  rate  of  energy  loss  in  zone  k  due  to  a  discrete  radiative 
process  (or  frequency  group)  P,  and  Cp^j  is  the  radiative  coupling  of  zone 
k  to  zone  j  for  that  process.  The  couplings  are  functions  of  opacity, 
integated  over  process  and  photon  path.  In  the  probabilistic  model,  a 
matrix  of  couplings  must  be  computed  for  each  bound-bound,  bound-free  and 
free-free  process;  for  the  multifrequency  model,  they  must  be  computed  for 
each  discrete  frequency.  In  this  way,  the  net  cooling  and  heating  by 
radiation  emission  and  absorption  among  the  various  zones  of  the  plasma  is 
accurately  accounted  for. 


(D)  Beam  Dynamics  and  Energy  Deposition 


We  assume  that  on  the  timescale  of  the  target  response  (several 
nanoseconds),  the  beam  maintains  uniform  flow,  that  is 


PB  UB  =  m 


J  °E 
e  Z„ 


(10) 


where  pg  and  ug  are  the  local  m~ss  density  and  axial  velocity  of  the 
beam.  The  mass  flux  ra  is  assumed  constant  up  to  the  point  the  beam  is 


stopped.  J  is  the  beam  current  per  unit  area;  mg  and  Zg  are  the  mass  and 
charge  of  a  beam  ion. 

For  a  relativistic  beam  with  negligible  thermal  spread,  the  kinetic 
energy  of  a  beam  ion  is  related  to  the  local  velocity  by 


mBc 


•'(l-Ug  /c) 


2  "  raBC 


(11) 


The  energy  is  also  related  to  the  stopping  power  by 


E  =  E  - 


,dE.  a 
dx^T  dX’ 


(12) 


where  EQ  is  the  initial  energy,  and  with  the  total  stopping  power  given  by 
the  sura  of  the  ion-free  electron,  ion-bound  electron  and  ion-plasma  ion 
interactions , 


/dE.  _  /dE.  /dc,\  j.  /dr.. 
(dJ}T  “  dx'e  +  (Fb  +  (dFi 


(13) 


(1)  Interaction  with  Bound  Electrons 


The  local  oscillator  model  (I.OM)15  representation  of  the  stopping 
cross  section  is  employed,  with  the  spatial  change  of  energy  given  by 


<3§>b  ' 


vv> 


(14) 


with  S.Q  in  ergs-cm' 


4tZ  ^  q*  7 

Sfe(V)  =  - - -  /  4*  r2  o(r)  Kq(t)  dr. 

mV 


(15) 


where  Zg  is  the  effective  charge  of  the  projectile  ion,  e  and  m  are  the 
electron  charge  and  mass,  V  is  the  projectile  velocity,  o(r)  is  the  local 
bound  electron  density  in  the  target  atom,  and  K  (r)  is  a  modified  zeroth 
order  Bessel  function. 


T 


h  clo(r)/raV2  , 


(16) 
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where  h  is  Planck's  constant  divided  by  2-x ,  and  q ( r )  is  the  local  plasma 
frequency  at  radius  r  within  the  atom. 

There  are  three  assumptions  underlying  the  LOM  •  The  first  is  that  a 
loss  function  can  be  defined,  dependent  only  on  the  local  electron  density 
in  the  target  atom.  The  second  is  that  the  longitudinal  dielectric 
response  can  be  represented  by  e(c j)  with  a  single  zero  at  >  ~  ■•)  ,  subject 
to  the  high  frequency  condition 

e(w)  ~  1  -  m2/;.)2  (17) 

o 

appropriate  to  free  electrons.  The  third  assumption  is  that  the  induced 
polarization  charge  is  spread  out  from  the  ion  trajectory  to  some  finite 
radius  of  the  order  of  the  de  Broglie  wavelength  h/mV.  This  last 
assumption  is  justified  by  the  adiabatic  argument  that  energy  will  be 
transferred  only  to  electrons  with  velocity  less  than  V.  A  wave  packet 
with  momentum  of  order  mV  would  have  a  spatial  spread  of  order  h/mV. 

The  electron  density  p(r)  is  taken  from  the  Thomas-Fermi  (TF)  model 
proposed  by  Zink.^ 

In  a  plasma  target,  the  bound  atomic  electrons  can  be  screened  from 

the  projectile  ions  by  the  free  electrons.  Thus,  the  LOM  must  be  modified 

to  take  this  shielding  effect  into  account.  For  an  electron  in  an  isolated 

atom,  the  maximum  impact  parameter  is  V/n,  where  n\  is  a  characteristic 

18 

frequency  of  motion  .  Plasma  screening  limits  this  parameter  to  the  Debye 
length  D.  The  expression  for  t  can  be  rewritten 

t  -  (h/mV)  (<>e(r)/V)  .  (18) 

Whenever  V/oi  (r)  >  D,  x  is  taken  to  be 

x  -  h/raVD  (19) 

Thus,  plasr.a  shielding  of  the  bound  electrons  is  taken  into  account  by 
limiting  the  argument  of  the  Bessel  function. 
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(2)  Interaction  With  Free  Electrons 


For  heated  target  materials,  the  atoms  become  Ionized,  and  stopping 
due  to  free  electrons  must  be  considered.  The  free  electron  stopping  power 
is  calculated  from 


<§>  'Hr1  sf<v>-  (20> 

e  I 

with 

4  2 

2”e  V  Z  D2  2 

S  (V)  -  - t  F(0  I  f.n  (1  +  -j)  +?n(l  +  4c  )  ‘  (21) 

mV  b 

where  Z  Is  the  average  charge  of  the  target  Ions,  D  Is  the  Debye  shielding 
length,  and 


K  -  (mBV2/2kT)1/2, 


(22) 


where  kT  is  the  electron  temperature,  and 


F(S)  =  erf(f)  -  2  £  e  /  A  .  (23) 

The  quantity  b  is  the  minimum  Impact  parameter  for  electron-ion 
scattering  and  is  given  by 

2 


b  -  MAX 


V 

i-V 

raV 


h_ 

2raV 


(24) 


i.e.,  the  maximum  of  either  the  classical  or  quantum-mechanical  impact 

parameter  defined  by  the  uncertainty  principle. 

The  first  term  in  S^(V)  is  the  short  range  ion-electron  binary- 

19 

encounter  scattering  term  and  is  taken  from  the  work  of  Campbell  ,  who 

20 

adapted  it  from  Brueckner  and  Brysk  .  The  second  terra  is  the  polarization 
term  and  is  taken  from  Pines  and  Bohra*"  .  For  distances  larger  than  D,  the 
plasma  acts  as  a  continuous  medium,  and  distant  collisions  cause  loss  of 
energy  by  the  excitation  of  plasma  oscillations,  which  appear  as  an 
oscillating  wake  behind  the  projectile. 
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(3)  Interaction  with  Plasma  Ions 

For  stopping  power  due  to  the  plasma  ions,  an  expression  taken  from 
Campbell1^  and  Mehlhorn^  is  used; 


(25) 


where 


Are4  Z  2  Z 

S  (V)  -  - * -  on  A  F(-:) 

mV 


where  £2  =■  i^V^kT  , 


2  ®I  _r ^ 

F( £ )  =  erf (£)  -  4-  (1  +  — )  ?e 

/Tt  mB 


»“  A  -  (1/2)  *»  (1  +  /  bol2n). 


bfflax  is  set  equal  to  D,  and  bm^n  is  given  by 


b  ,  -  MAX 

min 


1  v2  ’  2umpV 
pm  pV  P 


where  mp  is  the  proton  mass,  and 


mB  ml 


mB  +  mT 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


(4)  Enhanced  Stopping  Power  From  Collective  Effects 

The  stopping  power  of  a  dense  ion  beam  can  differ  from  the  sum  of 
single  particle  stopping  powers.  One  way  in  which  this  can  occur  is 
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through  phase  mixing  of  the  polarization  wakes  produced  in  the  ambient 
medium  by  the  beam  particles.  The  magnitude  of  this  particular  effect  is 
similar  to  that  given  by  McCorkle  and  lafrate  for  electron  beams, 

Sfe  -  Sq  (1  +  nB  a3/3)  (32) 


where  Sb  is  the  enhanced  stopping  power,  SQ  is  the  single  particle  stopping 
power,  ng  is  the  beam  number  density  and  a  is  the  ratio  of  the  beam 
velocity  and  local  plasma  frequency 


a 


V  /■,< 

V3'  pe 


(33) 


For  the  case  considered  below,  a  1  meV 

26  -2  — i 

flux  ng  Vg  ■  10  cm  sec  on  a  solid  density 

correction  to  the  stopping  power  is  quite  small. 


proton  beam  with  a 
aluminum  target,  the 


II.  RESULTS 

The  interaction  of  a  monoenergetic  1  MeV  proton  beam  with  a  planar 
aluminum  slab  of  15  ym  thickness  is  treated  in  this  investigation.  A  slab 
thickness  of  15  ym  is  chosen  because  this  distance  corresponds  to  the 
range  of  a  1  MeV  proton  in  cold,  solid  density  aluminum.'^  Hence,  the  beam 
is  totally  stopped  within  the  slab  and  should  deposit  all  of  its  energy 
within  the  target.  The  beam  is  assumed  to  consist  of  a  square-shaped  pulse 
of  10  nanosecond  duration  with  a  flux  of  10  protons/cm  -sec .  The  beam 

1  “5  O 

intensity  on  target  is  1.6XlCrJ  W/cm  ,  which  is  comparable  to  the 
intensities  available  with  current  devices.  These  calculations  were  made 
on  the  CRAY  computer  at  LANL. 

2  5 

In  an  earlier  report,  the  results  of  a  similar  calculation,  which 
employed  a  probabilistic  radiation  transport  formalism  as  described  In 
Section  I,  were  given.  It  was  postulated  that  the  probabilistic  method  had 
over-estimated  the  bound-free  and  free-free  radiative  cooling  of  the  slab, 
since  it  becomes  difficult  to  calculate  continuum  opacities  accurately  with 
the  probabilistic  model  when  they  become  very  large.  Some  limited 
calculations  carried  out  with  the  more  accurate  hybrid  model,  described 
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above,  confirmed  Chese  concerns.  Thus,  the  present  coupled  hydrodynamic- 
radiation  transport  calculation  was  carried  out  using  the  hybrid  scheme. 
This  is  the  first  time  the  hybrid  method  has  been  employed  in  a  full 
production  calculation.  Previously,  it  was  used  for  post-processing 
profiles  generated  with  the  probabilistic  or  multifrequency  methods,  and 
for  checking  the  accuracy  of  probabilistic  treatments. 

The  cost  of  a  full  multifrequency  production  calculation  would  be 
prohibitive.  The  aluminum  atomic  model  used  in  this  study  contains  about 
140  lines;  properly  resolving  them  with  perhaps  ten  frequencies  per  line 
would  require  1400  frequencies  in  addition  to  those  required  for  treatment 
of  the  continuum.  However,  a  probabilistic  treatment  of  line  radiation 
transport  remains  viable  in  the  present  case.  Although  this  treatment  will 
not  provide  details  of  the  emitted  line  profiles,  it  will  accurately 
describe  radiation  energetics. 

Hence,  the  hybrid  scheme,  while  sacrificing  spectral  detail,  permits  a 
cost  effective  self-consistent  hydrodynamic-radiation  transport  calculation 
to  be  performed  for  the  ion  beam-target  Interaction  study. 

Figure  2  shows  density  profiles  at  various  times  within  the  evolution 
of  the  Interaction.  A  striking  feature  is  the  spreading  of  the  target 
plasma  and  the  correspond  Lag  density  decrease  with  time.  From  an  initial 
width  of  15  -  and  a  density  of  2.7  grams/cm  ,  the  target  has  become  about 
0.3  cm  wide  with  densities  of  the  order  of  10  graras/cm  at  9.8  nsec  just 
before  the  proton  beam  is  turned  off.  The  density  gradient  smooths  out, 
and  the  plasma  becomes  more  uniform  with  time. 

Figure  3  shows  corresponding  profiles  of  temperature.  At  these 
densities  and  temperatures,  the  major  contribution  to  beam  stopping  comes 
from  the  free  electrons.  Indeed,  at  the  times  shown  in  Figs.  2  and  3,  the 
bean  is  entirely  stopped  in  the  low  density,  high  temperature  side  of  the 
plasma.  The  hack  side  of  the  target  is  heated  primarily  through  radiation 
transport,  which  is  substantially  larger  than  the  effect  of  thermal 
conduc  t  Lon . 

Figure  4  shows  profiles  of  density  and  temperature  at  14.8  nsec,  4.8 
nsec  after  the  proton  beam  has  been  shut  off.  The  plasma  is  now  about  0.52 
cm  wide  and  is  becoming  more  uniform.  While  the  temperature  of  most  of  the 
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plasma  has  fallen  considerably  since  beam  cutoff,  the  back  surface 
temperature  has  only  dropped  a  few  eV.  The  density  is  becoming  uniform. 

Figure  5  shows  the  energy  history  of  the  plasma  through  14.8  nsec. 
Note  that  the  total  energy  deposited  becomes  constant  at  10  nsec  when  the 
beam  shuts  off.  At  early  times,  the  energy  in  the  plasma  is  approximately 
equal  to  the  deposited  energy;  at  5  nsec,  the  radiated  energy  amounts  to 
about  20%  of  the  deposited  energy;  beyond  this  point,  the  fraction  of  the 
deposited  energy  retained  by  the  plasma  decreases.  By  10  nsec,  more  than 
40%  of  the  deposited  energy  has  been  radiated  away,  and  by  14.8  nsec,  more 
than  half  the  deposited  energy  has  been  radiated  away.  Figure  5  also  shows 
the  energy  partitioning  in  the  plasma.  Up  until  10  nsec,  all  of  the 
components  are  growing,  but  beyond  this  point,  the  ionization  and  thermal 
energies  decrease  as  the  plasma  cools.  The  kinetic  energy  continues  to 
grow  as  the  plasma  expands.  By  12  nsec,  the  radiated  energy  surpasses  the 
energy  retained  within  the  plasma.  At  times  greater  than  12  nsec,  the 
plasma  shows  a  net  loss  in  energy  due  to  radiative  cooling. 

Calculated  emission  spectra  from  the  front  and  rear  surfaces  of  the 
plasma  at  1.76  nsec  are  shown  in  Figs.  6  and  7.  At  this  early  time, 
radiative  cooling  is  not  a  significant  factor  as  may  be  seen  from  Fig.  5. 
There  is  much  more  energy  in  the  front  spectrum,  where  the  temperature  is 
higher.  Due  to  the  large  optical  depths  at  photon  energies  between  0.01 
keV  and  1.0  keV,  radiation  from  the  hotter,  less  dense  regions  does  not 
reach  the  back  side;  most  of  the  rear  side  emission  arises  from  the  plasn.. 
near  the  rear  surface. 

Figures  8  and  9  show  calculated  emission  spectra  from  the  front  and 
rear  sides  of  the  plasma  at  9.8  nsec.  Again,  there  is  more  radiated  power 
coming  out  of  the  hotter,  less  dense  front  surface  than  the  cooler,  more 
dense  rear  surface.  Free  bound  emission  from  the  L  edges  is  high  at  the 
front  (beam)  edge,  but  due  to  large  opt. cal  depths  in  the  0.3  keV  region, 
this  emission  is  greatly  attenuated  at  the  rear  side.  Around  1.0  keV,  the 
optical  depths  become  of  order  unity,  and  the  emission  rises,  then  falls 
again  due  to  K-edge  absorption.  Emissions  from  several  L-shell  lines  can 
be  seen  in  both  spectra,  as  can  K  line  emissions  in  the  1.5  keV  to  2.0  kev* 
range.  There  are  more  K  lines  appearing  on  the  front  side. 
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Figures  10  and  11  show  front  and  back  emission  spectra  at  14.8  nsec. 
These  spectra  are  more  similar  than  the  front  and  rear  spectra  at  9.8  nsec 
due  to  the  fact  that  the  plasma  is  more  uniform  at  this  time.  From  0.01 
keV  up  to  about  0.25  keV,  the  spectra  are  nearly  identical;  above  this 
frequency,  the  rear  spectrum  is  greatly  reduced  with  respect  to  the  front 
spectrum.  This  is  due  to  the  increased  density  near  the  rear  surface  and 
the  large  optical  depths  at  frequencies  from  about  0.25  keV  to  1.0  keV; 
i.e.,  for  frequencies  about  0.1  keV  and  up  to  about  0.25  keV,  emissions 
from  the  front  and  rear  surfaces  have  significant  contributions  from  deep 
within  the  plasma,  but  for  frequencies  greater  than  0.25  keV,  nearly  all 
the  emission  from  both  sides  originates  near  the  surface  due  to  the  large 
optical  depths.  At  these  frequencies,  the  front  surface  radiates  more 
energy.  Due  to  the  cooler  temperatures,  K-shell  lines  are  missing  from 
these  spectra. 

IV.  COMPARISON  WITH  CALCULATION  WITHOUT  RADIATION  TRANSPORT 

One  might  question  the  need  for  a  fully  self-consistent  hydro-dynaraic- 
radiation  transport  calculation.  Would  it  not  suffice  to  perform  a 
hydrodynamic  calculation,  and  post-process  the  resulting  temperature- 
density  profiles  to  obtain  the  radiation?  Alternatively,  could  not  a 
simple  radiation  transport  model  suffice  in  a  calculation  of  this  sort? 
Certainly  there  are  problems  where  one  of  these  alternatives  would  be 
adequate.  However,  for  the  current  beam-target  Interaction  problem, 
accurate  radiation  transport  is  essential.  To  demonstrate  this,  an 
identical  simulation  wi th  fully  self-consistent  ionization  physics,  but 
without  radiative  cooling  or  transport  was  carried  out. 

Comparisons  of  density  and  temperature  profiles  at  various  times  for 
both  sets  of  calculations  are  given  in  the  following  figures. 

Figure  12  shows  comparison  profiles  at  6  nsec.  On  the  beam  side  of 
the  plasma,  the  "purely  hydrodynamic  plasma"  (PHP)  calculation  yields  much 
higher  temperatures  and  correspondingly  lower  djnsities.  This  clearly 
demonstrates  the  effects  of  radiative  cooling  on  the  hyd rodynamic  response 
The  PH?  is  wider  at  this  time  due  to  the  absence  of  radiative  cooling. 

Figure  13  shows  similar  comparisons  at  10  nsec.  Similar  comments 
regarding  temperature  and  density  can  be  made  here.  In  this  situation,  the 
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back  surface  of  Che  plasma  with  radiation  included  has  expanded  further 
Chan  the  PH?  case  as  opposed  to  the  situation  in  Fig.  12.  This  is  chiefly 
due  to  the  effects  of  radiative  heating  on  the  rear  surface. 

Figure  14  presents  comparison  profiles  at  15  nsec,  5  nsec  after  the 
beam  has  been  shut  off.  Similar  observations  regarding  temperature  and 
density  can  be  made;  the  PHP  still  produces  the  wider  plasma.  A  striking 
feature  of  Fig.  15  is  the  relative  uniformity  of  the  profiles  for  the 
radiation  case  as  compared  to  the  PHP  case.  The  PHP  case  still  retains 
steep  gradients  near  the  rear  surface,  while  the  radiation-hydro  plasma  is 
approaching  an  isothermal  condition  with  fairly  uniform  density. 

Figure  15  gives  the  energy  history  of  the  plasma  along  with  the  energy 
partitioning  for  the  PHP  case.  The  plasma  energy  is  equal  to  the  deposited 
energy  at  all  times.  The  thermal  energy  is  much  higher  than  in  Fig.  5  for 
the  radiation-hydro  case;  this  indicates  that  radiation  cooling  occurs  at 
the  expense  of  the  plasma  thermal  energy.  This  is  also  indicated  by  the 
much  higher  temperatures  in  the  PHP  curves  in  Figs.  12-14.  The  ionisation 
energies  are  not  very  different  in  the  two  calculations.  The  PHP  case  has 
higher  kinetic  energy  at  later  times;  this  accounts  for  the  greater  widths 
of  the  hydrodynamic  plasma. 

These  comparisons  of  the  PHP  and  hydro-radiation  transport  cases 
clearly  establish  the  importance  of  radiation  transport  effects  in  hot 
plasmas.  Any  realistic  treatment  of  hot  plasmas  must  include  self- 
consistent  radiation  effects. 


V.  COMPARISON  WITH  LASER-TARGET  INTERACTION 

It  Is  of  interest  to  make  some  qualitative  comparisons  with  a 
simulation  of  a  laser-aluminum  target  interaction  at  roughly  the  same 

r 

incident  energy  Intensity.  Duston,  et.  al.  0  studied  the  interaction  of  a 

13  2 

1.06  micron  laser  beam  at  an  Intensity  of  10  W/cm '  and  a  long  pulse 
length  (3  nsec  full  width  at  half  maximum)  with  an  3  urn  thick  aluminum 
foil.  Fig.  16  is  taken  from  that  study. 

In  the  laser  case,  after  the  initial  blow  off  plasma  is  formed,  the 
beam  does  not  penetrate  into  the  target  beyond  the  critical  density 
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surface,  whereas  Che  ion  beam  penetrates  deep  into  the  target,  at  least  at 


early  times.  In  both  cases,  the  dense  part  of  the  plasma  at  later  times  is 
heated  predominantly  by  radiation  transport  rather  than  thermal 
conduction.  The  laser  case  shows  much  lower  densities  and  higher 
temperatures  on  the  front  side  of  the  plasma,  whereas  on  the  backside  the 
laser  case  shows  much  higher  density  and  lower  temperatures.  Also,  by 
comparing  Figs.  2  and  16,  it  can  be  seen  that  the  laser  beam  calculation 
produces  much  sharper  density  and  temperature  gradients.  These  results 
substantiate  the  view  that  proton  beams  are  more  efficient  than  laser  beams 
for  heating  up  the  interior  and  backside  of  targets;  this  is  due  to  deeper 
penetration  into  the  target  by  the  ion  beams. 


VI.  SUMMARY 

A  fully  self  consistent,  one-dimensional  treatment  of  the  inter-action 
of  a  one-MeV  proton  beam  with  an  aluminum  slab  has  been  presented.  A  fully 
self-consistent  hydrodynaraic-ionizat ion-radiation  transport  model  was 
used.  The  atomic  populations  as  functions  of  density  and  temperature  were 
obtained  from  a  collisional-radiative-equilibrium  assumption.  Radiation 
transport  was  calculated  using  a  hybrid  scheme. 

Profiles  of  the  plasma  density,  temperature,  and  average  charge 
resulting  from  these  calculations  have  been  presented.  Comparisons  with 
calculations  omitting  radiation  transport  show  that  radiation  can 
significantly  affect  plasma  evolution,  leading  to  substantially  different 
temperature  and  density  profiles,  and  must  be  included  in  any  self- 
consistent  treatment  of  ion  beam-target  interactions. 

Qualitative  comparisons  of  the  ion  beam-slab  interaction  with  a 
treatment^”  of  a  laser-beam  slab  interaction  shows  that  the  laser  produces 
higher  temperatures  and  lower  densities  on  the  front  (beam)  side  of  the 
plasma,  while  the  proton  beam  is  more  efficient  at  heating  up  the  back  side 
due  to  deeper  penetration  into  the  target. 

Since  these  treatments  of  the  beam-target  interaction  are  one¬ 
dimensional,  energy  flow  normal  to  the  direction  of  the  beam  has  been 
neglected;  in  addition,  self-field  effects  have  not  been  considered.  These 


91 


r  - 

J-.-r  -  aaaaa 


effects  could  reduce  the  range  significantly,  and  alter  the  magnitude  and 
shapes  of  the  plasma  profiles  given  here.  A  two-dimensional  hyd rodynami- 
code  is  under  development  and  it  is  planned  to  use  this  model  in  the  beam- 
target  interaction  study  when  it  becomes  available. 

Even  with  these  limitations,  however,  the  above  results  should  he 
useful  in  particle  beam  lethality  and  vulnerable  studies. 
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Fig.  1.  Simple  three-level  diagram  used  to  depict  the  atomic  processes 
Included  In  this  ionization  model.  S  is  collisional  and  photoioniza 
tion,  ct^  is  radiative  recombination  and  three-body  recombination,  is 

dielectronic  recombination,  X  is  electron  collisional  excitation,  Y  is 
collisional  deexcitation,  A  is  spontaneous  decay,  and  0  and  i'j  *  are 

stimulated  absorption  and  emission. 
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Fig.  13.  Comparison  of  temperature  and  density  profiles  at  10  nsec  with 
and  without  radiation  transport. 
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Fig.  14.  Comparison  of  temperature  and  density  profiles  at  15  nsec  with 
and  without  radiation  transport. 
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Fig.  15.  Energy  history  of  the  beam-target  interaction  when  radiation 
transport  is  omitted. 
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Fig.  16.  Temperature,  density,  and  radiative  cooling  profiles  at  about  5 
nsec  in  the  laser  heated  foil  case  (Ref.  26). 
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F. 


Enhanced  Radiation  Sources 


I.  INTRODUCTION 

Exploring  schemes  and  configurations  for  enhanced  plasma  radiation 
will  require  a  more  complete  understanding  of  the  nature  of  radiation 
matter  interactions.  Since  we  are  interested  in  producing  x-rays  for 
simulation  purposes,  it  is  irrelevant  whether  the  x-rays  are  incoherent  or 
coherent.  Naturally,  because  of  the  power  requirements  and  the  lack  of  a 
working  x-ray  laser,  combined  with  the  relative  ease  of  producing 
incoherent  x-ray  radiation  all  of  DNA's  sponsored  research  efforts  in  this 
area  have  potential  application.  However,  it  appears  that  it  may  be 
possible  to  develop  a  coherent  x-ray  source  using  current  gas  puff 
technology.  It  is  in  this  atmosphere  that  we  have  investigated  a  scheme 
that  could  be  implemented  experimentally  using  pulse  power  technology. 
Since  our  initial  goal  is  to  understand  the  physics  issues  and  determine 
the  plasma  conditions  favorable  to  produce  a  coherent  x-ray  source,  we  have 
simplified  things  by  assumming  a  laser  driven  system.  In  future 
investigations  the  laser  will  be  replaced  by  a  PRS. 

For  many  years  now,  there  have  been  a  variety  of  proposed  schemes 
suggesting  ways  to  create  a  population  inversion  and  subsequent  gain  in  the 
x-ray  region*-^.  One  such  scheme  involves  the  radiated  flux  from  one 
plasma  pumping  on  another  plasma  to  create  conditions  which  result  in 
coherent  radiation  in  the  soft  x-ray  region^-^.  There  are  several  cases 
where  the  coincidence  between  the  pumping  and  absorbing  lines  is  close 
enough  to  merit  strong  consideration  as  a  possible  lasing  combination.  Two 
such  schemes  utilizing  a  Si  XIII  -  A1  XII  system  and  a  Na  X  -  Ne  IX  system 
are  modeled  in  detail  in  reference  1  (hereafter,  referred  to  as  ADW) .  The 
Na  -  Ne  scheme  is  particularly  interesting  because  the  lines  match  to 
better  than  one  part  in  10^  at  a  wavelength  of  11*. 

The  proper  conditions  in  the  Na  -  Ne  system  is  created  by  "tuning"  the 
hydrodynamics  and  thermal  properties  of  the  two  component  plasma^.  This 
tuning  acts  both  to  create  a  tenous  neon  blowoff  and  to  provide  an  energy 
source  to  heat  the  sodium.  The  sodium  subsequently  acts  as  a  backlighting 
flashlamp  for  the  neon. 
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For  instance,  absorption  of  a  laser  operating  at  low  intensity  can 
ablate  the  neon  at  a  low  temperature,  typically  about  100  eV.  Second, 

energy  can  be  deposited  into  the  sodium  backlighter  and  locally  heat  this 
region.  This  energy  deposition  can  occur  on  the  opposite  side  relative  to 
the  neon  cooling  to  prevent  over  heating  the  neon  or  it  may  occur  off  to 
the  side  of  the  expanding  neon  plasma  plume.  The  optimum  configuration  has 
not  yet  been  determined.  The  rapid  rise  in  temperature  induced  by  this 
deposition  will  then  conduct  through  the  sodium  and  create  the  proper 
conditions  for  photo-pumping  the  cooler  neon.  Of  course,  within  a  finite 
time  scale  this  thermal  wave  propagates  through  the  neon  and  effectively 
burns  out  any  inverted  population.  The  plasma  system  then  remains  in  an 
isothermal  condition. 

Several  constraints  in  this  problem  are  brought  about  by  considering 
the  plasma  properties  which  are  needed  to  achieve  gain  in  a  resonant  photo- 
pumped  medium1.  First  the  line  absorption  coefficients  are  linear  in  ion 
density.  Therefore,  for  a  fixed  optical  depth,  the  width  d  of  the  ablated 
neon  can  be  scaled  as  d  »  A/N^  where  A  is  a  constant  and  is  the  ion 
density.  The  maximum  optical  depth  before  significant  reduction  in 
population  inversion  occurs  has  been  determined  by  Hagelstein^  to  be 
approximately  5.  The  approximation  given  in  ADW  is  that  this  requires  an 
ablated  blowoff  with  a  scale  length  not  greater  than  10  ym  for  ion 
densities  on  the  order  of  10  cm  .  This  scale  length  can  increase  to 
greater  than  0.1  cm  for  ion  number  densities  on  the  order  of  101®  cm-^. 
For  the  inverted  levels  discussed  in  ADW,  the  gain  coefficient  for  the  Ne 
IX  2p*P  -  3d^D  and  2p*P  -  4d*D  transitions  peak  near  an  ion  density  of  10^° 
cm  J.  The  gain  subsequently  drops  off  dramatically  to  less  than  unity  at 
101®  cm  However,  the  gain  for  the  3d^D  -  4f^F  transition  peaks  near 
1019  cm-^  and  is  still  greater  than  10  at  101®  cm"^. 

A  second  constraint  is  placed  on  the  temperatures  of  the  two  plasmas. 
The  temperature  of  the  lasing  neon  must  be  kept  low  enough  to  prevent 
burnout  of  the  Ne  IX  population.  At  a  density  of  10^  cm~^  the  maximum 
fractional  population  of  Ne  IX  occurs  at  approximately  65  eV  in  the 
presence  of  the  pump  radiation.  On  the  other  hand,  the  pumping  plaraa  (in 
this  case  Na)  is  most  favorable  when  it  is  an  optically  thick  plasma  with  a 
kinetic  temperature  of  about  300  eV.  In  short,  one  would  like  a  hot  dense 
pumping  plasma  irradiating  a  cooler,  optically  thin  lasing  plasma1. 
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As  indicated  above,  the  density  requirements  for  the  3-2  and  4-2 

transitions  differ  from  the  requirements  of  the  4-3  transition.  Therefore, 
we  can  separate  the  following  discussion  into  two  parts.  First,  we  will 
examine  the  3-2  and  4-2  scaling  requirements.  We  consider  cryogenic  neon 
with  density  1.2  gm/cm^  on  the  outside  of  a  sodium  sample.  This  neon  layer 
is  subsequently  ablated  off  the  sodium  to  dynamically  set  up  the  proper 
plasma  conditions  for  the  photo-pumped  system.  In  order  to  prevent  photon 
trapping  and  meet  the  conditions  given  above,  a  layer  of  2.79x10  gm  of  Ne 
is  needed.  This  assumes  a  uniform  blowoff  density  of  10  cm  with  a 

width  of  10  gm.  If  we  assume  that  the  Ne  plasma  is  initially  produced  by 
laser  absorption,  then  the  ablation  velocity  will  be  on  the  order  of  10^ 

cm/sec.  At  these  volocities  expansion  out  to  10  gm  occurs  on  a  time  scale 

of  roughly  100  psec.  Both  the  pumping  material  and  the  lasing  material 
must  reach  the  proper  conditions  for  gain  in  this  time  frame.  The  strong 
dependence  of  the  gain  coefficient^  on  density  results  in  the  existence  of 
a  window  of  approximatly  10  nsec  before  the  3-2  and  4-2  transitions  yield  a 
gain  coefficient  less  than  unity. 

On  the  other  hand,  the  4-3  transition  has  a  gain  coefficient  of  1.0  at 
approximately  10  cm  and  hence,  the  expansion  time  can  be  up  to  1  gsec. 
At  peak  gain  (i.e.,  lO1^  cm-^)  this  allows  1  nsec  for  the  blowoff  to 
develop.  Again  this  assumes  a  uniform  expansion  velocity  of  10^  cm/sec. 

We  now  consider  the  propagation  of  the  thermal  wave  from  the  area  of 
energy  absorption  in  the  sodium  through  the  neon  plasma.  The  electron  heat 
flux  is  given  by 


q  -  -vcfr  .  (1) 

The  conductivity  <  can  be  written  in  the  form 

<  -  CnekXcve  (2) 

where  C  is  a  constant  of  order  unity,  ng  is  the  electron  density,  k  is 
Boltzmann's  constant,  X is  the  electron-ion  collision  mean  free  path,  and 
vg  is  the  electron  thermal  speed.  Equation  (1)  can  be  rewritten  in  the 
form 


113 


(3) 


«  “  CnekTe(  i;  )ve(  “eT  > 


where  L^.  is  the  temperature  gradient  length  and  e^,  is  given  by 
$T/ J  )  .  If  L^,  <  Xc,  the  heat  flux  will  be  carried  in  the  direction 
-e^  at  a  rate  which  is  greater  than  the  electron  thermal  speed.  This  is  an 
unphysical  result  and  indeed,  the  derivation  of  the  conductivity  given  by 
Eq.  (2)  is  predicated  on  the  assumption  that  A^  <  L^.  This  problem  of 
steep  gradients  has  been  studied  extensively  in  the  laser  fusion 
literature®-** .  The  approach  which  has  been  used  most  often  in  connection 
with  this  difficulty  has  been  to  limit  the  flux  q  whenever  Ac  >  fLT  to 


q  -  nekW  ~®T  > 


(4) 


where  f,  the  so-called  "flux  limiter",  is  less  than  one®-*^.  Therefore,  we 
can  base  our  estimate  of  the  propagation  speed  of  the  thermal  wave  on  some 
limited  value  of  vg.  The  ratio  of  the  flux  limited  value  of  the  thermal 
wave  speed  to  the  sound  speed  cg  will  be  given  by 


fv 


e 


c 


s 


(5) 


where  we  have  used  the  atomic  weight  for  neon  and  an  adiabatic  index  of 
5/3.  Recent  theoretical  and  experimental  results  have  put  the  value  of  f 
close  to  0..  for  laser  power  levels  which  are  appropriate  for  our 
conditions  l3-15>  Assuming  Ne  IX  (  i.e.,  Z  ■  8  ) ,  Eq.  (4)  gives  fv£  =  5.2 
c8  and  the  thermal  wave  be  expected  to  overtake  the  expanding  plasma  up  to 
Mach  numbers  of  roughly  5. 

Max,  McKee,  and  Mead*^  have  studied  in  detail  the  properties  of  steady 
state  flow  beyond  the  critcal  surface  in  plasmas  created  by  laser 
absorption  and  the  dependence  of  this  flow  on  the  flux  limiter  in  spherical 
geometry.  While  this  analysis  is  not  appropriate  for  the  conditions  which 
occur  during  the  early  phase  of  rapid  energy  deposition,  it  is  helpful  in 
gaining  insight  into  the  physics  during  the  later  absorption  and  blowoff 
stages.  The  parameters  discussed  above  fall  into  the  "semi-saturated" 
regime.  That  is,  the  heat  flux  will  be  classical  up  to  the  point  M  *  1.76 
where  M  is  the  Mach  number*^.  This  point  will  occur  at  1.38  critical 
surface  radli*^.  Beyond  Chat  point  heat  flow  will  be  flux  limited  and  the 
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Mach  number  will  increase  slowly  thereafter.  Hence,  we  can  expect  that  the 
thermal  front  should  overtake  the  expanding  plasma.  Based  on  these 
results,  we  can  make  an  estimate  of  the  propagation  time  for  the  thermal 

IQ  _  *3 

wave  to  penetrate  the  Ne  plasma.  If  we  assume  an  ion  density  of  10  cm 
and  choose  an  optical  depth  of  5,  then  the  total  distance  to  be  traversed 
in  the  neon  plasma  is  approximately  100  uni.  If  the  value  given  above  for 
fv0  is  chosen  with  cg  =  1.5  x  102  cm/sec,  the  thermal  wave  will  propagate 
halfway  into  the  neon  in  70  psec.  This  puts  a  rough  limit  on  the  time  span 
over  which  effective  gain  will  be  achieved. 

In  order  to  assess  the  validity  of  this  scheme  in  a  more  detailed 
fashion,  we  are  using  a  1-D  Lagrangian  hydrodynamics  code  to  model  the  Na- 
Ne  plasma  combination.  This  code  treats  the  electron  and  ion  temperatures 
separately.  Equation  of  state  models  are  used  for  both  the  ions  and  the 
electrons.  Of  course,  because  of  the  low  densities,  the  plasma  in  the 
subcritical  region  can  be  expected  to  behave  nearly  as  a  perfect  gas.  We 
have  assumed  a  uniform  planar  slab  of  sodium  coated  with  a  thin  layer  of 
neon.  The  Ne  layer  contains  enough  particles  so  that  for  a  uniform 
expansion  of  width  10  yra  the  ion  density  will  be  1020  cm-3.  At  solid  or 
liquid  densities  this  gives  a  submicron  layer  of  uniform  thickness.  The 
initial  zone  sizes  are  varied  in  the  transition  region  so  that  the  zonal 
masses  are  nearly  equal  at  the  Na-Ne  interface.  Velocities  and  positions 
are  defined  on  zone  boundaries.  Density,  temperature,  and  laser  absorption 
are  computed  at  zone  centers.  The  electron  temperature  is  flux  limited 
according  to  Eq.  (4). 

In  a  preliminary  fashion,  we  have  been  abale  to  make  an  estimate  of 
the  gain  produced  by  this  configuration  at  each  of  these  wavelengths  using 
the  line  center  gain  coefficient  curve  for  the  Ne  IX  3d*D  -  4f^F  transition 
given  in  ADW  and  making  the  simplifying  assumption  of  a  65  eV  plasma 
temperature  in  the  neon.  This  also  assumes  that  the  flux  at  line  center 
for  the  11  A  pump  is  6.6xl03  ergs/(cm2-sec-Hz)  as  inferred  from  a 
microballoon  implosion  at  Rochester  .  We  consider  the  case  where  energy 
is  deposited  locally  in  the  sodium  by  a  high  powered  laser  after  ablating 
the  neon  with  a  low  energy  prepulse.  The  sodium  target  is  considered  to  be 
several  microns  thick.  There  is  no  absorption  in  the  neon  of  this  high 
powered  pulse.  This,  might  represent  laser  absorption  from  the  side  or  it 
might  crudely  represent  laser  absorption  occurring  in  the  rear  side  of  the 
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sodium  (i.e.,  the  side  away  from  the  neon).  At  the  present  time,  we  are 
conducting  a  more  detailed  analysis  where  the  high  power  laser  is  absorbed 
non-locally  (i.e.,  inverse  bremstrahlung)  in  the  rear  side  of  a  thin  sodium 
slab.  Results  from  this  study  will  be  given  in  a  future  report.  An 

average  density  can  be  computed  by  taking  the  N’e  plasma  length  at  the  time 

when  the  200  eV  thermal  front  propagates  out  of  the  system  for  10  w/cm“ 
deposition  into  the  sodium. 

We  have  performed  post  processed  CRE  (Collisional  Radiative 

Equilibrium)1^  calculations  of  the  sodium  pump  intensity  which  would  be 

incident  upon  the  neon  lasing  material  for  a  representative  snapshot 

density  and  temperature  profile.  This  has  allowed  us  to  compare  with  the 

3  ergs 

Na  pumping  line  intensity  reported  in  ADW  of  6.6x10  — — -  at  line 

cm  -sec-Hz 

center  obtained  from  the  Rochester  implosion  experiment1  .  The  two 
profiles  are  shown  in  Figure  1.  Self  absorption  is  evident  at  line 
center  but  not  in  the  Rochester  profile  which  showed  no  self  absorption 
within  the  experimental  resolution1  of  1.5  eV.  Except  for  the  self 
absorption,  this  line  profile  closely  resembles  the  microballoon  implosion 
profile  of  the  Rochester  experiment. 

The  results  of  ADW  (where  the  neon  is  assumed  fully  bathed  in  pump 
radiation)  showed  that  the  Rochester  profile  produced  a  peak  gain 

coefficient  of  nearly  200  cm-1.  Proportionality  arguments  would  then  yield 

a  gain  coefficient  of  roughly  2.5  cm-1,  for  the  case  shown  where  the  neon 

is  illuminated  on  only  one  side  by  the  sodium. 

We  must  stress,  that  extensive  modeling  of  the  pulse,  deposition 
configuration,  and  sodium  target  thickness  will  most  likely  produce 
different  line  center  intensities  incident  on  the  neon  plasma  than  the 

crude  results  reported  here.  Of  crucial  importance  will  be  the  correct 
modeling  of  the  energy  deposition  into  the  sodium  backlighter  and  the 
subsequent  interaction  with  the  neon  plasma.  Using  the  best  gain 

coefficient  projected  from  post  processing  the  Na  pump,  it  appears  that 

gain  may  be  measurable  in  an  experiment  of  this  type.  It  may,  however,  be 
more  advantageous  to  employ  a  cylindrical  or  spherical  configuration  in 
order  to  provide  compressional  heating  of  the  pump  to  the  requisite 
temperature  at  high  density. 


cm -gcc-Hz 


Figure  1 
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